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ABSTRACT 


The  diffusion  coefficients  of  Nb*9S  in  PG  were  measured  between  1450 
and  1950^0  in  three  different  types  of  PG  which  were  chosen  to  represent  sub¬ 
strate  nucleated,  combination  substrate  and  regeneratively  nucleated,  and 
regeneratively  nucleated  microstructures.  The  diffusion  coefficients  parallel 
and  perpendicular  to  the  deposition  surface  were  determined  in  as -deposited 
and  partially  graphitized  PG.  Autoradiography  was  used  to  observe  the  distri¬ 
bution  of  diffusing  atoms .  The  structure  of  the  PG  was  characterized  by  the 
application  of  many  techniques  of  metallography,  x-ray  diffraction,  small 
angle  scattering  of  x-rays,  and  density  determinations.  A  mathematical 
analysis  of  diffusion  in  layered  structures  was  made  and  correlations  between 
the  structure  and  diffusion  parameters  noted. 

The  activation  energies  for  Nb  diffusion  in  PG  depended  only  upon  direction 
and  were  independent  of  microstructure  and  state  of  graphitization .  The  activa¬ 
tion  energy  for  Nb  diffusion  parallel  to  the  deposition  plane  was  96.2  kcal/mole 
and  that  perpendicular  to  it  was  73.1  kcal/mole.  The  pre -exponential  factors 
for  a  given  direction  differed  only  slightly  for  the  various  PG's  used  in  this  study; 
thes^  factors  were  altered  by  the  graphitizing  heat-treatment.  Parallel  to 
the  deposition  plane,  the  pre -exponential  factors  were  all  approximately 
3.2  X  10^  cm^sec  for  the  as-deposited  PG’s  and  1.8  x  10^  cm^/sec  for  the 
heat-treated  PG's .  For  diffusion  perpendicular  to  the  deposition  plane,  these 
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factors  were  approximately  1.2  x  10  cm  /sec  for  the  as -deposited  PC's  and 
-3  2 

6.5  X  10  cm  /sec  for  the  heat-treated  material.  Autoradiography  demon¬ 
strated  that  diffusion  perpendicular  to  the  layer  planes  occurs  in  the  primary 
cone  boundaries  while  the  radioactivity  is  uniformly  distributed  in  the  auto¬ 
radiographs  taken  for  diffusion  parallel  to  deposition  plane . 

Diffusion  mechanisms  within  the  crystallites  were  rejected  on  the  bases  of 
valence  and  relative  atom  size.  An  examination  of  mechanisms  external  to  the 
crystallites,  combined  with  an  analysis  of  the  relative  values  of  the  pre- 
exponential  factors,  revealed  that  the  fundamental  jump  distance  for  diffusion 
parallel  to  the  deposition  plane  is  the  crystallite  diameter,  while  the  jump  dis¬ 
tance  perpendicular  to  the  deposition  plane  is  the  layer  plane  separation  distance . 
This  model  was  confirmed  quantitatively  by  applying  the  results  of  an  analysis  of 
diffusion  in  the  primary  cone  boundaries  and  by  utilizing  elementary  interstitial 
diffusion  theory . 

In  another  set  of  experiments  with  Re-PG  diffusion  couples,  a  metastable 
carbide  of  Re  was  discovered.  The  carbide  has  a  simple  cubic  crystal  structure 
with  a  unit  cell  edge  of  2.8+  A  and  reacts  with  Re  to  form  a  eutectic  at  ~2000^C 
(±200°C). 
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Chapter  1 


INTRODUCTION 

A.  Preliminary  Comments 

Massive  pyrolytic  graphite  <PG),  a  recent  development  in  the  field  of 
chemical  vapor  deposition  (1,2)  has  aroused  considerable  interest  because  of 
its  high -temperature  properties.  PG  exhibits  a  very  high  tensile  strength  at 
temperatures  in  excess  of  3000°C  (3,4,5),  is  highly  anisotropic  in  its  electri¬ 
cal  (6)  and  thermal  (3)  conduction  properties,  and  has  nearly  theoretical 
density  (1).  The  major  limitations  of  PG  for  high  temperature  service  are 
its  low  oxidation  resistance  at  high  temperature  and  its  low  shear  stress  parallel 
to  the  plane  of  the  deposit.  Efforts  to  alleviate  these  problems  and  to  produce 
material  with  controllable  electrical  properties  have  been  made  by  co-depositing 
metals  with  PG  (7,8);  to  date,  however,  no  other  methods  have  been  investigated. 
Although  the  diffusion  of  metals  in  PG  could  be  important  relative  to  its  behavior 
in  high  temperature  environments,  no  studies  of  this  kind  had  been  reported  at 
the  time  this  investigation  began . 

In  addition  to  the  above,  PG  posseses  a  unique  and  interesting  anisotropic 
structure,  and  many  techniques  have  been  developed  recently  to  characterize 
it.  Thus,  diffusion  studies  in  PG  arc  of  value  from  both  practical  and  theoreti¬ 
cal  viewpoints .  The  practical  value  of  such  studies  lies  in  the  possibility  of 
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diffusion -alloying  PG  to  improve  its  oxidation  and  mechanical  properties  and 
to  alter  its  electrical  properties.  The  theoretical  value  lies  in  studying  diffu¬ 
sion  in  this  unique  anisotropic  structure,  using  the  recently  developed  techniques 
for  structural  characterization  to  further  the  understanding  of  the  structure  of 
PG  and  to  determine  the  underlying  mechanisms  of  metal  diffusion  in  PG . 

B.  Technological  Aspects  of  PG 

Graphite  has  a  sublimation  temperature  considerably  in  excess  of  3000°C, 
the  triple  point  occurring  at  approximately  3700°C  and  105  Kg/mm^  (9, 10, 11). 
The  high  sublimation  temperature  combined  with  a  very  high  degree  of  covalent 
bonding  within  the  layers  makes  graphite  a  highly  refractory  material,  i.e.,  a 
material  t^at  retains  its  strength  and  is  capable  of  providing  service  at  high 
temperature.  Normal  sjrnthetic  graphites  begin  to  lose  strength  around  2500°C 
(3),  but  PG  actually  begins  to  increase  in  tensile  strength  parallel  to  the  layer 
planes  at  this  temperature  (3, 5, 12).  Furthermore,  PG  is  anisotropic  in  all  of 
its  properties,  being  a  good  conductor  of  electricity  and  heat  parallel  to  its 
layer  planes  and  a  fairly  good  insulator  perpendicular  to  them  (3). 

These  highly  refractory  and  highly  anisotropic  properties  make  PG  a  good 
material  for  many  high  temperature  uses,  including  rocket  nozzles,  and  re-entry 
vehicles.  Nuclear  applications  include  such  items  as  control  rod  housings  (12), 
structural  members  in  high  temperature  nuclear  reactors,  and  coatings  on  small 
particles  of  reactor  fuel  (13, 14). 
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The  diffusion  of  metals  in  PG  can  play  an  imporunt  role  in  many  of  the 
above  applications;  for  example,  it  may  be  possible  to  diffusion -alloy  PG  with 
metals  to  alter  its  mechanical  and  oxidation  properties .  It  has  been  observed 
that  high  temperature  deformation  of  PG  occurs  by  crystallites  sliding  past  each 
other  (5)  and  that  oxidation  of  graphite  occurs  principally  at  the  edges  of  the 
crystallites  (15).  One  of  the  conclusions  of  this  investigation  is  that  the  diffusion 
of  metals  in  PG  takes  place  along  the  crystallite  surfaces,  suggesting  that  the 
high  temperature  properties  of  PG  could  be  improved  by  diffusion  alloying.  A 
proper  selection  of  metal  could  interfere  with  both  crystallite  sliding  and  the 
peripheral  oxidation  mechanisms  to  improve  PG's  performance  at  high  tempera¬ 
ture.  As  noted  above,  alloying  of  PG  by  co-deposition  of  metals  is  presently 
being  accomplished  (7,8). 

C .  Scientific  Aspects  of  PG 

PG  is  a  unique  material.  It  offers  a  high  degree  of  anisotropy  in  both  defect 
structure  and  properties  and  nearly  theoretical  density  with  very  low  porosity 
and  permeability.  Furthermore,  it  is  the  closest  approximation  presently  avail¬ 
able  to  large  single  crystals  of  graphite.  The  layer  structure  of  graphite  and  the 
method  of  production  make  for  an  unusual  anisotropic  defect  structure  which  will 
be  important  in  mass  transfer  as  well  as  many  other  properties. 

This  structure,  which  will  be  discussed  fully  in  Chapter  III,  is  unique  to 
PG.  However,  some  of  the  conclusions  reached  in  this  study  will  also  apply  to 
diffusion  in  other  oriented,  layer  structures  such  as  BN  and  Bi2Te2.  Diffusion 
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in  ^  has  not  been  investigated,  but  its  structure  is  very  similar  to  that  of 
graphite  (9).  The  diffusion  of  copper  in  single  crystal  Bi2Tc^  was  found  to 
be  highly  anisotropic  (16),  and  impurity  diffusion  in  oriented,  polycrystalline 
material  might  well  behave  in  the  same  fashion  as  that  in  PG. 

D.  Brief  Review  of  Diffusion  in  Graphites 

Diffusion  in  graphites  will  be  reviewed  in  more  detail  and  referenced  in 
Chapter  II;  however,  a  brief  summary  will  be  given  here  to  acquaint  the  reader 
with  the  literature .  Carbon  self -diffusion  in  an  ideal  graphite  lattice  has  been 
treated  theoretically  and  C-14  tracer  diffusion  experiments  have  been  carried 
out,  with  relatively  little  success,  on  natural  single  crystals  and  bulk  synthetic 
graphite.  The  migration  of  carbon  atoms  displaced  by  radiation  damage  has 
been  studied  extensively  and  some  data  on  carbon  diffusion  have  been  derived 
from  electron  microscopy  studies  on  quenched-in  dislocation  loops  in  graphite. 

There  have  been  several  investigations  of  the  migration  of  fission  products 
from  individual  graphite  fuel  cells,  but  these  experiments  do  not  generally  lend 
themselves  to  unequivocal  interpretation  because  the  graphites  consist  of  vary¬ 
ing  amounts  df  "amorphous"  and  "crystalline"  carbon  with  widely  varying  orien¬ 
tations  and  also  arc  subject  to  radiation  damage.  These  factors,  combined  with 
the  presence  of  a  large  number  of  elements  in  differing  amounts,  make  a  funda¬ 
mental  interpretation  regarding  the  diffusion  mechanisms  impossible. 

Under  more  controlled  conditions,  the  diffusion  rates  of  U,  B,  and  Th  in 
various  forms  of  carbon  have  been  studied  and,  recently,  an  extensive  study  of 
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the  diffusion  of  Ag,  Ni,  U,  Th,  and  Ra  in  PG  and  other  forms  of  carbon  was 
made.  To  date,  this  latter  study  is  the  only  investigation  of  metal  diffusion  in 
graphite  using  standard  radioactive  tracer  techniques  and  resulting  in  reason¬ 
ably  clear  interpretations.  Although  this  study  showed  how  the  diffusion  kinetics 
varied  for  different  atoms,  the  role  of  the  structure  of  the  PG  and  giaphitcs  in 
the  diffusion  process  was  not  investigated  except  that  it  was  shown  that  the  dif¬ 
fusion  rate  of  Ni-63  in  extremely  well  annealed  PG  is  quite  small  compared  to 
that  in  as -deposited  PG. 

E .  Structure  of  PG 

The  structure  of  PG  has  recently  been  investigated  in  some  detail,  as  will 
be  discussed  in  Chapter  III.  X-ray  diffraction  yields  information  concerning 
the  degree  of  perfection  within  the  crystallites,  their  size  in  the  <a>  and  <c> 
directions,  the  degree  of  preferred  orientation,  and  the  variation  of  crystallite 
size  and  perfection  with  orientation  relative  to  the  deposition  plane  normal . 

Small  angle  scattering  (SAS)  provides  information  concerning  the  void  size  and 
distribution  and  transmission  electron  microscopy  can  be  used  to  observe  small 
aggregates  of  crystallites.  Optical  microscopy  provides  information  concerning 
the  manner  in  which  the  PG  has  nucleated  and  grown  during  its  manufacture  and, 
because  of  the  optical  anisotropy  of  PG,  the  relative  orientation  of  major  struc¬ 
tural  features  can  be  determined.  Density  measurements  and  externally  con¬ 
nected  porosity  measurements  permit  a  comparison  with  the  theoretical  x-ray 
density  and  thereby  provide  a  link  between  the  x-ray  diffraction  and  SAS  informa¬ 
tion.  Techniques  have  been  developed  for  all  of  the  above  measurements  and 
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their  combined  application  to  a  particular  PG  results  in  a  very  detailed  descrip¬ 
tion  of  the  structure  that  can  be  correlated  with  many  properties . 

F.  Description  of  Research  Problem 

This  program  was  initiated  to  examine  the  possibility  of  diffusion -alloying 
PG  and  to  determine  the  role  of  structure  in  metal  diffusion  in  PG 

The  purpose  of  the  investigation  was  to  determine  how  the  structure  of 
PG  affects  the  diffusion  of  metal  atoms  in  PG.  To  do  this,  three  different  kinds 
of  PG  were  selected  and  the  diffusion  coefficient  of  Nb-9S  was  measured  perpen¬ 
dicular  and  parallel  to  the  deposition  plane,  in  the  as -deposited  and  in  one  heat- 
treated  condition  for  each  of  the  PG's.  The  structure  of  the  PG's  was  character¬ 
ized  by  analysis  of  x-ray  diffraction  and  small  angle  scattering  data,  by  optical 
and  electron  microscopy,  and  by  density  determinations.  Correlations  with  the 
defect  structure  of  the  PG  deposits  and  the  changes  of  this  structure  with  heat- 
treatment  were  made  with  the  observed  diffusion  coefficients  in  the  two  principal 


directions . 


Chapter  11 

REVIEW  OF  DIFFUSION  IN  GRAPHITES 

A .  Introduction 

Carbon  self-diffusion  has  been  studied  both  theoretically  and  experimen¬ 
tally,  while  high  temperature  impurity  diffusion  in  graphites  has  been  examined 
only  experimentally .  The  different  theoretical  treatments  of  carbon  self¬ 
diffusion  arc  in  fair  agreement,  but  the  experimental  determinations  of  self¬ 
diffusion  arc  open  to  question.  In  addition,  most  of  the  high  temperature 
impurity  diffusion  studies  have  very  questionable  value.  Many  of  these  investi¬ 
gations  have  been  concerned  with  the  migration  of  fission  fragments  in  poorly 
characterized  graphites  under  very  complex  conditions;  the  few  less  compli¬ 
cated  experiments  are  marked  by  a  considerable  lack  of  data  or  rather  question¬ 
able  interpretations.  Only  the  recent  data  of  Wolfe,  McKenzie,  and  Borg  (17) 
and  possibly  that  of  Hennig  (18)  warrant  serious  discussion  regarding  the  mechan¬ 
isms  of  diffusion  in  graphites;  the  other  data,  however,  will  be  reviewed  here 
for  general  reference. 

B.  Theoretical  Analyses  of  Self-Diffusion  in  Graphite 

Carton  self-diffusion  in  graphite  has  been  theoretically  analyzed  by 
Dienes  (19).  Using  arguments  based  on  the  variation  of  bond  strength  with 
distance,  he  calculated  the  activation  energies  of  self-diffusion  for  the  layer 
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vacancy  mechanism,  the  direct  interchange  of  carbon  atoms  within  a  layer 
plane,  and  the  interlayer  or  interstitial  migration  mechanism.  Dienes  defines 
interstitial  migration  as  the  motion  of  atoms  between  the  layers  of  the  graphite 
lattice  and  not  as  passage  through  the  hexagon  rings .  Although  physically  there 
is  a  hole  in  the  center  of  each  hexagon,  the  electronic  interactions  probably 
make  this  a  very  high  energy  path.  A  better  term  for  the  motion  of  atoms 
between  the  layers  would  be  ’’interlayer  migration”  which  will  be  used  here  to 
describe  any  atoms  residing  between  the  layers.  The  calculated  activation 
energy  for  the  interlayer  diffusion  mechanism  included  a  term  for  removing 
an  atom  from  the  edge  of  the  layer  and  a  term  involving  repulsive  force  inter¬ 
action  of  the  interlayer  atom  with  the  layers.  The  repulsive  force  was  repre¬ 
sented  as  an  exponential  function  of  distance  and  the  parameters  were  evaluated 
using  experimentally  determined  elastic  constants . 

All  of  Dienes'  calculations  were  based  on  early  data  by  Pauling  (20)  con¬ 
cerning  the  dependence  of  the  bond  strength  on  the  carbon-carbon  equilibrium 
distance  for  singly  bonded,  resonance  bonded,  and  doubly  and  triply  bonded 
carbon  atoms.  Kanter  (21, 22)  pointed  out  that  more  recent  data  reported  by 
Coulson  (23)  on  the  bond  strength  and  carbon-carbon  distance  raise  Dienes’ 
value  of  79.6  kcal/mole  for  the  bond  energy  in  graphite  (C-C  distance  of  1.42A) 
to  113  kcal/mole.  Table  1  contains  the  results  of  Dienes’  calculations,  cor¬ 
rected  by  Kanter  (21).  The  very  high  value  for  the  formation  of  an  interlayer 
carbon  atom  should  be  noted.  More  recently,  Dienes  (24)  has  estimated  the 
energy  for  motion  of  an  interlayer  carbon  atom  to  be  about  2.6  kcal/mole. 
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Table  1 

CALCULATED  ACTIVATION  ENERGIES  FOR  CARBON 
SELF-DIFFUSION  IN  GRAPHITE(») 


Process 

Energv  (kcal/mole) 

Formation  of  vacancy 
within  hexagon  layer 

170 

Motion  of  vacancy  within 
layer 

93 

Total  energy  for  vacancy 
mechanism 

263 

Direct  interchange  within 
a  layer 

113 

FormatiOL.  of  an  interlayer 
Interstitial  carbon  atom 

467 

Motion  of  an  interlayer 
carbon  atom 

not  discussed 

(a)  Originally  calculated  by  Dienes  (19)  uid  corrected  by 
Ranter  (21,22) 


The  energies  of  formation  and  motion  of  interlayer  C  and  Xe  in  the 
graphite  lattice  have  been  calculated  by  Iwata  et  al .  (25)  using,  as  a  model  of 
the  graphite  layers,  thin  elastic  sheets  subjected  to  bending  due  to  the  presence 
of  the  interlayer  atom.  Their  results  are  as  follows: 


Process 


Activation  Energy— kcal/mole 


Formation  of  interlayer  57.6 

carbon  atom 

Formation  of  interlayer  346 

xenon  atom 

Motion  of  interlayer  0 . 368 

carbon  atom 

Motion  of  interlayer  0.691 

xenon  atom 
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The  value  obtained  by  Twata  et  al .  for  forming  the  carbon  interlayer  atom 
differs  considerably  from  that  by  Dienes.  This  difference'  is  partly  due  to 
the  fact  that  Iwata  did  not  include  tite  energy  to  remove  an  atom  from  the 
surface  and  partly  because  Iwata 's  calculations  allow  for  relaxation  of  the 
layer  sheets  about  the  interlayer  atom  while  Dienes'  calculations  do  not . 
Experimental  determination  of  these  energies  should  be  possible  by  observ¬ 
ing  the  annealing  behavior  of  neutron-irradiated  graphite;  however,  this 
behavior  is  quite  complex  (26,27).  Many  different  processes  occur  during 
the  annealing  of  radiation-damaged  graphite  and  interpretation  of  the  change 
of  a  property  v/ith  time  at  temperature  is  not  simple  (27).  In  any  case,  the 
annealing  of  relatively  light  radiation  damage  can  be  accomplished  below  about 
200°C,  indicating  that  the  simplest  types  of  radiation  damage  have  rather  low 
migration  energies.  Using  a  model  based  on  the  migration  of  single  interlayer 
carbon  atoms  produced  by  neutron  irradiation,  Schweitzer  (28)  determined  the 
activation  energy  for  motion  of  these  atoms  to  be  1.6  ^  0.1  kcal/mole  between 
77  and  490°K. 

The  energies  for  the  formation  and  motion  of  layer  vacancies  have  been 
estimated  by  Baker  and  Kelly  (29)  using  observations  with  the  electron  micro¬ 
scope  of  dislocation  loops  formed  in  graphite  single  crystals.  The  loops  were 
formed  by  quenching  the  crystals  from  3000®C  and  then  annealing  between  1200 
and  1400°C  at  which  temperatures  the  excess  vacancies  should  coalesce  and  con 
dense  into  basal  plane  dislocation  loops.  With  arguments  based  on  the  distance 
vacancies  must  travel  to  coalesce  into  dislocation  loops.  Baker  and  Kelly 
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estimated  the  vacancy  formation  energy  to  be  about  78  kcal/mole  and  the 
motion  energy  to  be  69  kcal/mole,  giving  an  activation  energy  for  self¬ 
diffusion  of  150  kcal/mole.  Wooley  (30)  has  taken  the  ideas  of  Baker  and 
Kelly  somewhat  further  and  calculated  the  energy  to  form  a  vacancy -interlayer 
atom  pair  in  graphite  as  138  ±55  kcal/mole.  Several  configurations  such  as 
interlayer  clusters  were  also  examined.  In  addition,  Wooley  cites  a  value  for 
the  motion  energy  of  an  interlayer  carbon  atom  of  less  than  9  kcal/mole. 

Hennig  (31)  was  unable  to  reproduce  the  results  given  by  Baker  and  Kelly 
and  came  to  the  conclusion  that  the  lower  limit  for  the  formation  of  a  vacancy 
should  be  about  140  kcal/mole.  Hennig's  observations  were  supported  some¬ 
what  by  theoretical  calculations  made  by  Coulson  and  co-workers  (32,  33)  on 
the  energy  to  form  layer  vacancies.  Their  calculations,  based  on  electron 
configurational  energies  in  a  carbon  macro-molecular  layer  with  and  without 
vacancies,  gave  246  kcal/mole  for  the  energy  to  form  a  single  vacancy  (32). 

They  further  showed  that  a  divacancy  is  more  stable  than  two  single  vacancies 
by  110  kcal/mole  of  divacancy  (33). 

To  summarize,  calculations  and  estimates  of  the  energy  to  form  a  single 
layer  vacancy  in  graphite  vary  from  78  to  246  kcal/mole  with  the  more  likely 
value  being  in  the  middle  of  the  range.  The  energy  to  move  a  vacancy  is 
estimated  to  be  between  69  and  93  kcal/mole.  Estimates  of  the  energy  to  form 
an  interstitial  carbon  atom  range  from  58  to  467  kcal/mole,  but  the  lower  value 
is  more  probable  because  it  takes  relaxation  into  account  and  does  not  include  a 
term  for  removing  an  atom  from  one  of  the  layers.  The  energy  of  an  atom  between 

the  edges  of  two  crystallites  is  probably  rather  high,  and  the  energy  to  remove 
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it  from  its  environment  would  not  be  as  great  as  removing  an  atom  from  the 
edge  of  a  layer.  The  formation  of  an  interlayer  foreign  atom,  such  as  Xe, 
requires  about  350  kcal/mole,  and  the  energy  for  motion  of  an  interlayer 
atom,  C  or  Xe,  is  very  small,  being  less  than  10  kcal/mole  and  more  likely 
less  than  1  kcal/mole. 

There  are  no  specific  theoretical  treatments  for  the  diffusion  of  atoms 
other  than  carbon  in  the  graphite  lattice  except  that  given  for  Xe .  Many  metals 
can  enter  into  directional  bonds  with  carbon  and  other  atoms,  as  in  the  case  of 
the  metal -organic  compounds  (34);  however,  the  large  size  of  nearly  all  metal 
atoms  would  pr;5vent  diffusion  by  direct  substitution  and  direct  vacancy-atom 
exchange  in  the  graphite  lattice.  Boron  is  the  only  ’’metallic"  atom  that  might 
be  expected  to  enter  the  graphite  lattice  substitutionally . 

C .  Carbon  Self  Diffusion  -  Experimental  Data 

TVo  studies  of  tracer  self -diffusion  in  graphite  have  been  reported. 

Feldman  et  al.  (35)  studied  tracer  self  diffusion  in  National  Carbon  AUF  grade 
graphite  by  depositing  a  slurry  of  carbon  and  water,  containing  C-14,  on  the 
ends  of  small  graphite  rods  and  annealing  for  various  times  near  2000°C.  The 
cylindrical  surfaces  of  the  rods  were  machined  before  and  after  annealing  and 
sections  along  the  rod  axis  were  cut  with  a  microtome.  It  was  not  clear  whether 
the  data  were  best  interpreted  as  log  concentration  vs.  distance  to  the  first  or 
second  power.  After  some  justifications  regarding  activation  energies,  the 
data  were  interpreted  on  the  basis  of  an  analysis  of  simultaneous  grain  boundary 
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and  volume  diffusion  proposed  by  Fisher  (36).  Assuming  a  value  for  the 
activation  energy  for  volume  diffusion  of  90  kcal/mole,  Feldman  et  al. 
obtained  75.4  kcal/mole  as  the  activation  energy  for  grain  boundary  diffusion. 

Kanter  (21,22)  studied  tracer  self  diffusion  in  natural  graphite  crystals 
by  heating  them  in  an  atmosphere  of  radioactive  carbon  monoxide  and  then 
measuring  the  concentration  profile  by  a  controlled  burning  process .  The 
profiles  could  not  be  measured  conveniently  in  any  other  way  due  to  the  small 
size  of  the  natural  graphite  crystals  (  ~  1mm  diameter).  Experiments  utiliz¬ 
ing  conditions  that  approximated  the  instantaneous  source  solution  to  the  dif¬ 
fusion  equations  did  not  yield  meaningful  information,  but  experiments  approxi¬ 
mating  a  constant  surface  composition  produced  seemingly  better  data:  an 
activation  energy  of  163  kcal/mole  was  determined  for  layer  plane  diffusion  in 
graphite,  based  on  the  assumption  that  carbon  self -diffusion  occurs  only  per¬ 
pendicular  to  the  c-axis  of  the  crystal. 

D.  Impurity  Diffusion 

Investigation  of  impurity  diffusion  in  graphites  includes  studies  of  the  forma¬ 
tion  of  lamellar  compounds,  migration  of  fission  fragments  in  graphite,  and  direct 
measurement  of  the  diffusivity,  using  different  methods  for  measuring  the  con¬ 
centration  profiles  produced  by  diffusion. 

A  phenomenon  related  to  impurity  diffusion  is  the  formation  of  lamellar  com¬ 
pounds  with  graphite  in  which  foreign  atoms  enter  into  the  graphite  lattice  at  cer¬ 
tain  position  between  the  layers  (9, 37).  The  process  of  entering  the  lattice  is 
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termed  intercalation  and  such  lamellar  compounds  are  also  termed  intercala¬ 
tion  compounds.  Formation  of  these  compounds  does  not  have  a  direct  bearing  on 
high  temperature  diffusion  in  graphites  but  these  compounds  will  be  discussed 
briefly  for  the  sake  of  completeness.  A  great  many  elements  and  compounds  will 
intercalate  the  layer  structure  of  graphite  and  form  lamellar  or  intercalation  com  - 
pounds.  The  formation  of  the  compounds  requires  the  migration  or  diffusion  of  the 
element  or  compound  into  the  graphite  lattice;  however,  the  process  is 
complicated  by  the  formation  of  the  new  compound  and  therefore  kinetic 
studies  of  the  formation  of  these  compounds  would  not  be  directly  applicable 
to  the  study  of  metal  diffusion  in  graphites.  Lamellar  compounds  can  be 
formed  directly  with  the  alkali  metals  when  they  are  molten,  and  with  the 
halogens  when  heated  to  400  to  500°C  (37).  In  addition,  aqueous  solutions 
of  many  metal  chlorides,  oxides,  and  sulfides  (38)  will  intercalate  the 
graphite  structure.  Although  the  process  of  lamellar  compound  formation 
is  not  well  understood,  one  theory  is  that  the  metal  atom  must  be  sufficiently 
ionized  to  accept  a  tt  electron  from  the  graphite  before  it  can  enter  the  region 
between  the  layers  (38).  Two  phenomena  occur  with  the  formation  of  lamellar 
compounds .  The  first  is  a  large  expansion  of  the  lattice  (9)  perpendicular  to 
the  layer  planes  which  is  accompanied  with  some  visible  swelling  of  the  graphite 
body.  The  second  is  the  marked  tumescence  (9)  produced  if  the  intercalated 
graphite  structure  is  heated  to  a  few  hundred  degrees  C .  Expansions  of  ten 
fold  are  not  uncommon. 
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Studies  of  the  migration  of  fission  fragments  in  graphite  are  all  based  upon 
measuring  the  fraction  of  activity  remaining  for  any  or  all  the  fission  fragments. 
These  studies  can  be  divided  into  two  groups,  those  in  which  only  the  concentra¬ 
tion  vs.  time  at  temperature  data  are  reported  and  those  in  which  the  authois 
thought  the  data  were  of  sufficient  quality  to  warrant  computing  activation  ener¬ 
gies.  In  the  first  group  (39,40,41)  the  loss  of  Sr,  Y,  Zr,  Nb,  Mo,  Ru,  Rh,  Te, 

I,  Cs,  Ba,  Ce,  and  most  of  the  actinides  and  lanthanides  was  measured  for 
various  times  at  temperatures  between  900  and  2600°C .  The  general  conclusion 
was  that  the  more  volatile,  lower  melting  elements  migrated  out  of  the  graphite 
rather  rapidly  and  the  more  refractory  metals  Zr,  Nb,  Mo,  Ru,  Rh,  and  U 
demonstrated  very  little  tendency  to  leave  graphite  when  heated  as  long  as 
10  minutes  at  2400°C. 

Table  2  contains  a  summary  of  the  studies  of  the  migration  of  fission 
f  gments  in  graphite  in  which  the  authors  felt  they  could  determine  an  activa¬ 
tion  energy.  Rates  for  the  Arrhenius  plots  were  determined  by  a  variety  of 
methods,  the  most  common  one  being  to  plot  the  log  of  the  time  required  to  reach 
a  certain  relative  composition  vs.  reciprocal  temperature.  Although  the  numbers 
given  in  these  investigations  have  some  engineering  utility,  they  are  of  little  use 
in  determining  the  fundamental  mechanisms  of  diffusion  in  graphite  because  of 
the  complicated  and  varied  structure  of  the  graphite  used  and  the  presence  of 
many  differant  chemical  species. 
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Table  2 


SUMMARY  OF  STUDIES  OF  THE  MIGRATION  OF 
nSSION  FRAGMENTS  IN  GRAPHITE 


Diffusing 

Species 

Ai^rent  Activation 
Energy,  Q 
(kcal/mole) 

Reference 

Remarks 

Total  Activity 

20.7 

Q  varied  with  tern- 

Ba-140 

68.2 

42 

perature,  being  lower 

at  higher  temperatures 

Sr-189 

26.4 

As-77 

32 

Br-83 

27 

Pd-109,  112 

9 

Cd-115 

30 

^43 

Sb-127 

50 

Te-18 

18 

1-131,  133 

42 

Xe-135 

49 

44 

D^  =  3  X  10"^ 

0 

1-131 

46  -  49 

45 

Cs 

18 

D^  ~  lO"^ 

Sr 

10 

46 

®  -4 

D  ~  10  ^ 

0  -4 

Ba 

9 

D^  ~  10  ^ 

Xe-135 

23 

Natural  Graphite  1 

Xe-135 

23 

Natural  Graphite  2 

Xe-135 

27 

►  47 

Natural  Graphite  3 

Xe-135 

23 

Synthetic  Graphite 

Xe-135 

49 

Amorphous  Carbon 

A 

» 
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The  remaining  impurity  diffusion  studies  discussed  in  this  section  were 
made  without  the  additional  complications  associated  with  the  presence  of 
many  different  elements  but,  for  various  reasons,  the  results  are  subject 
to  question . 

3 

The  diffusion  of  U  from  an  U-impregnated  graphite  (0.25  gm  U/cm  ) 
into  an  unimpregnated  graphite  (National  Carbon  grade  AGOT)  has  been 
studied  by  Loftness  (48).  The  diffusion  coefficient  was  measured  at  three 
temperatures  and  the  data  represented  by 

=  1.28  X  lO'^exp  -53,000/RT 


The  data  were  inadequate  because  only  2  or  3  points  were  determined  in  each 
composition  profile.  Diffusion  of  U  into  graphite  from  a  UCij  deposit  on  the 
surface  was  examined  by  Loch  et  al.  (49)  who  postulated  two  mechanisms  for 
the  process.  His  concentration  profiles  were  divided  into  two  regions:  a  low- 
concentration  high -penetration  tail,  and  a  high -concentration  low -penetration 
beginning.  The  large  penetration  portion  was  attributed  to  migration  in  the 
graphite  pores  and  the  steeper  profile  to  volume  diffusion.  The  volume  diffu¬ 
sion  portion  was  analyzed  and  the  data  represented  by 


7.6  X  10"^  exp  -87,000/RT 


between  1800  and  2300”C.  Above  2300°C,  the  apparent  diffusivity  was 
observed  to  increase  sharply  and  no  good  explanation  was  given. 
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Kislyi  and  Samsanov  (50)  measured  the  diffusion  of  B  in  synthetic  graphite, 
using  standard  sectioning  techniques  combined  with  chemical  analysis,  and 
found  the  following  relationship 

^B(G)  ^  3.02  exp  -57,250/RT 

to  hold.  They  claimed  extensive  substitutional  solubility  (up  to  26  wt  %)  of  B 
in  graphite  based  on  a  smooth  increase  in  the  (0006)  x-ray  diffraction  peak  of 
graphite  with  increasing  B  content.  However,  no  other  x-ray  evidence  was 
cited,  and  a  shift  in  an  (OOOf)  line  would  more  likely  signify  interlayer  pene¬ 
tration.  Aside  from  this,  the  diffusion  data  seem  to  be  reasonable  but  very 
little  information  was  given. 

Truitt  et  al.  (51)  measured  the  diffusion  of  Th  in  synthetic  and  pyrolytic 
graphite.  Specimens  were  coated  with  Th-232  by  ionically  bombarding  them 
with  singly  charged  ions,  and  following  the  diffusion  anneals,  the  specimens 
were  exposed  to  a  thermal  neutron  flux  to  produce  Th-233  which  decayed 
rapidly  to  Pa-233.  The  specimens  were  then  sectioned,  and  the  residual  Th 
content  was  taken  to  be  proportional  to  the  Pa-233  assay.  The  authors  state 
that  the  data  "suggest  60  and  70  kcal/mole,  respectively,  as  the  energy  of 
activation  for  diffusion  of  Th  in  porous  and  pyrolytic  graphite,  c -direction." 
They  further  note  the  similarity  of  predicted  diffusivities  for  U  in  graphite 
(48,49)  and  carbon  in  graphite  (35).  The  authors  do  not  mention  that  the  latter 
data  for  carbon  in  graphite  were  interpreted  on  the  basis  of  grain  boundary  dif¬ 
fusion  while  their  Th  data  were  interpreted  as  volume  diffusion.  In  a  note 
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added  in  proof,  the  diffusion  of  Th  perpendicular  to  the  c-axis  in  PG  was 
given  as  2.5  orders  of  magnitude  greater  than  that  parallel  to  the  c-axis. 

E.  Contributions  of  Wolfe,  McKenzie,  and  Borg 

Wolfe,  McKenzie  and  Borg  (WMB)  (17)  recently  completed  the  first  investi¬ 
gation  of  metal  diffusion  in  graphites  under  carefully  controlled  conditions . 

They  measured  the  diffusion  coefficients  of  Ag,  Ni,  U,  Th,  and  Ra  in  PG 
(perpendicular  and  parallel  to  the  c-axis),  National  Carbon  grade  ZTA  graphite 
(a  highdensity  synthetic  graphite),  ar..l  Speers  grade  S-700  (a  porous,  low  density 
graphite).  Diffusion  specimens  were  prepared  by  evaporating  an  aqueous  solu¬ 
tion  of  a  salt  of  the  tracer,  by  flashing  from  a  heated  filament,  and  by  vacuum 
evaporation.  The  technique  of  deposition  by  evaporation  of  an  aqueous  salt 
solution  was  not  suitable  for  the  more  porous  synthetic  graphites  due  to  capil¬ 
lary  action  but  functioned  quite  well  for  the  PG.  The  difficulty  with  the  salt 
de]X)sit  is  the  lack  of  knowledge  regarding  the  chemical  state  of  the  diffusing 
species.  WMB  did  not  find  any  anamolies  caused  by  this  deposition  technique, 
which  is  to  be  expected  when  one  considers  the  high  temperature  instability  of 
many  metal  salts  and  the  high  affinity  of  most  metals  for  carbon. 

Diffusion  anneals  were  made  in  vacuum  and  sectioning  accomplished  by 
attaching  the  sample  to  a  piston-in-cylinder  arrangement  and  abrading.  The 
amount  removed  and  penetration  distance  were  determined  by  the  loss  of  weight 
of  the  sample  and  piston.  Penetration  profiles  were  established  by  counting 
the  activity  remaining  after  each  sectioning.  In  addition,  autoradiographs  were 
taken  to  demonstrate  any  non-uniform  distribution  of  the  diffusing  atoms. 
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A  summary  of  WMB's  data  is  given  in  Table  3.  Besides  these  summviry 

data,  WMB  also  noted  that,  in  PG,  there  is  a  slight  increase  in  the  c-axis 

diffusion  coefficient  perpendicular  to  the  deposition  plane  near  the  substrate 

surface  (see  Chapter  III)  as  compared  with  that  near  the  terminal  deposition 

surface.  In  autoradiographs  of  Ag,  Ni,  and  U  in  PG,  WMB  observed  that, 

when  diffusion  occurred  perpendicular  to  the  deposition  plane,  the  diffusing 

specie  was  concentrated  at  the  cone  boundaries  (see  Chapter  III  and  IX);  while 

parallel  to  the  deposition  plane,  the  diffusing  Atoms  are  uniformly  distributed. 

In  the  synthetic  graphites,  the  autoradiographs  showed  inhomogenieties  on  a 

very  small  scale  (~lmm)  but  this  did  not  seem  to  affect  the  results.  WMB 

also  found  that  the  diffusion  coefficient  of  Ni  in  very  well  annealed  PG  (3000°C, 

.4 

15  minutes)  was  about  3  x  10  times  that  of  as -deposited  PG. 

Some  general  conclusions  concerning  the  data  of  WMB  are  as  follows: 

•  Lower  melting  metals  (Ag,  Ni)  have  lower  activation  energies 
(Q)  and  much  higher  diffusivities  than  do  the  higher  melting 
metals  at  a  given  temperature. 

•  For  Ni  and  U,  >  Q  but  for  Th,  <  Q  . 

^a  ^c  a  c 

•  For  Ni,  U  and  Th,  but  by  varying  amounts. 

•  For  diffusion  perpendicular  to  the  deposition  plane,  the  diffusing 
species  tends  to  concentrate  at  the  cone  boundaries  of  the  PG  but 
for  diffusion  parallel  to  the  deposition  plane,  it  seems  to  be 
uniformly  distributed. 


♦ 
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Table  3 

SUMMARY  OF  THE  DATA  OF  WOLFE.  MCKENZIE  AND  BORG<*) 

D  =  D 
o 


Element 

Graphite 

Direction 

Temp. 

Range 

CO 

“o 

cm^/sec 

Q 

kcal/mole 

Avg.  Deviation 
From  Predicted 
Diffusivity 
(%) 

Ag-110 

PG 

<a> 

600-  800 

9. 28  X  10^ 

64.  3  ±2.6 

6.4 

Ni-63 

PG 

<a> 

550-1000 

1.02  X  10^ 

47.2  ±  2.2 

9.3 

Ni-63 

PG 

<c> 

750-1050 

2.2 

53. 3  ±4. 8 

14.8 

U-232^ 

PG 

<a> 

1400-  2200 

6.76  X  10® 

115.0  ±  2 

6.0 

U-232 

PG 

<c> 

1800-2200 

3. 85  X  10® 

129.5  ±  5 

25.0 

Th-228 

PG 

<a  > 

1650-2150 

1.33  X  10® 

145.4  ±  3 

10.5 

Th-228 

PG 

<c> 

1800-2300 

2.48 

114.7  ±  3 

5.1 

Ra-224^®^ 

PG 

<a> 

1650-2150 

1.18  X  10^ 

99.4  ±  3 

5.2 

U-232 

ZTA 

— 

1600-2100 

6.22  X  10® 

133.1  ±  3 

7.9 

U-232 

S-700 

— 

1600-2100 

9.29  X  10® 

136.6  ±  3 

16.0 

Th-228 

ZTA 

i.  and  II 
to  extru¬ 
sion 

direction 

1889 

D  /D  =2. 
y  X 

68 

(a)  Ref.  17 

(b)  U-232  freshly  prepared  to  avoid  complications  caused  by  daughter  products 

(c)  Impurity  present  in  Th-228  tracer 


•  The  D^'s  for  U  in  synthetic  graphite  are  appreciably  higher  than 
either  in  PG  but  the  activation  energies  in  synthetic  graphite 
are  comparable  to  that  parallel  to  the  deposition  plane  in  PG . 

•  Diffusion  of  Ni  parallel  to  the  deposition  plane  in  very  well 
annealed  PG  is  considerably  less  than  in  as-deposited  PG. 

F.  Contribution  of  Hennig 

Hennig  (18)  has  recently  reported  data  on  the  diffusion  of  B  in  graphite 
single  crystals.  The  measurements  are  based  on  the  influence  B  has  on  the 
chemical  activity  of  graphite  and  the  technique  is  essentially  one  of  etch  pit 
decoration  and  quantitative  electron -metallography .  Several  arguments  were 
put  forth  that  B  enters  the  graphite  lattice  substitutionally.  The  diffusion 
coefficients  of  B  in  graphite  were  determined  between  1760  and  2290°C  as 

D„  =  6320  exp  (-157, 000/RT) 

lie 

=  7.1  exp  (- 153,000,  RT) 

The  agreement  with  Kanter's  carbon  self-diffusion  data  (21,22)  (Q  =  163,000 
kcal/mole)  was  noted,  but  Hennig  did  not  think  that  B  Jiffuses  by  a  vacancy 
mechanism.  The  question  of  mechanism  was  deferred  until  his  work  on 
lattice  vacancy  mobilities  is  complete. 
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G .  Summary 

Theoretical  calculations  for  self -diffusion  in  graphite  are  far  from 
complete  and  considerable  differences  exist  between  the  various  authors. 
Furthermore,  even  though  graphite  single  crystals  larger  than  a  few  milli¬ 
meters  are  not  available,  there  have  been  no  theoretical  treatments  or  even 
a  listing  of  the  many  possible  planar  defect  diffusion  mechanisms  such  as  along 
the  layer  edges,  down  the  edges,  across  the  surfaces  of  crystallites,  and  adja¬ 
cent  to  the  the  crystallites  but  in  the  region  where  the  bonding  differs  from  that 
of  nearly  ideal  hexagonal  graphite.  Some  other  mechanisms  for  diffusion  within 
the  crystallites  which  have  not  been  examined  are: 

•  Interstitial  diffusion  along  the  c-axis  by  going  through  the 
center  of  the  benzene  rings  which  compose  the  layers  of  graphite 

•  Vacancy -atom  interchange  between  layers,  especially  at  internal 
defects 

•  Atomic  transfer  from  one  layer  to  the  next  at  the  edges  of  internal 
defects  consisting  of  many  missing  atoms  arranged  as  "vacancy" 
benzene  rings  (52). 

Table  4  lists  the  probable  values  of  the  activation  energies  of  the  processes 
which  have  been  treated  theoretically . 

Tracer  self-diffusion  experiments  have  not  been  very  successful.  The  only 
activation  energy  with  any  degree  of  reliability  is  that  due  to  Kanter  (21,22)  of 
163  kcal/nrole,  which  was  subject  to  many  assumptions.  Hennig's  activation 
energies  for  B  in  graphite  of  155  kcal/mole  lend  considerable  credence  to 
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Table  4 


ESTIMATES  OF  ACTIVATION  ENERGIES  FOR 
VARIOUS  DIFFUSION  PROCESSES  IN  GRAPHITE 


Process 


Activation  Energy 
(kcal/mole) 


Formation  of  layer  vacancy  140  ±  60 

Motion  of  layer  vacancy  80  ±  20 

Formation  of  divacancy  ~  150 

Direct  interchange  within  layer  ~  100 

Formation  of  interlayer 

carbon  atom  >  60 

Formation  of  interlayer  xenon 

atom  (4. 4  A  diameter)  ~350 

Motion  of  interlayer  atom, 

carbon  or  xenon  ~  1 


Kanter's  data,  but  a  careful  tracer  experiment  with  large  single  crystals, 
j/crhaps  produced  pyrolytically  at  very  high  temperature,  is  still  needed. 

The  understanding  of  metal  diffusion  in  graphites  is  even  less  than  that 
for  self-diffusion.  The  many  studies  of  the  migration  of  fission  fragments 
have  been  hopelessly  complicated  by  radiation  damage,  the  complex  nature 
of  synthetic  graphites,  and  the  presence  of  many  chemical  species.  With 
the  exception  of  WMB's  data,  studies  of  metal  diffusion  in  graphite  are  marked 
by  an  absence  of  careful  characterization  of  the  graphite,  a  lack  of  use  of  the 
full  spectrum  of  analytical  tools,  such  as  autoradiography,  and  generally  poor 
data.  Even  if  attention  had  been  given  to  these  details,  the  complexity  of  syn¬ 
thetic  grapltite  would  have  made  fundamental  interpretations  difficult.  Although 
WMB's  investigation  was  generally  of  high  quality,  it  did  not  delve  deeply  into 
the  structural  dependence  of  metal  diffusion  in  graphites  but  did  demonstrate 
that  PG  is  the  least  complicated  graphite  available  for  study. 


} 
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Chapter  111 


STRUCTURE  OF  PYROLYTIC  GRAPHITE 


A .  Introduction 

PG  is  one  of  many  forms  of  pure  carbon.  Thus,  the  structural  chemistry 
of  carbon  must  be  understood  before  examining  PG  whose  structure  is  directly 
related  to  the  manner  in  which  it  is  made,  and  whose  properties  depend  upon 
that  structure.  This  structure  can  be  characterized  using  metallographic  and 
x-ray  techniques,  with  additional  information  available  through  density  and 
porosity  measurements.  Studies  of  the  structure-sensitive  properties  of  PG 
and  their  change  with  heat -treatment  also  provide  information  concerning  the 
nature  of  the  PG  deposit. 

In  this  chapter,  the  structural  chemistry  of  carbon  will  be  discussed,  as 
will  the  techniques  for  manufacturing  PG.  Following  this,  the  details  of  the 
structural  characterizations  of  PG  using  metallography,  x-ray  diffraction,  and 
other  techniques  will  be  outlined,  and  a  short  discussion  of  graphitization  in 
PG  will  be  given. 

B.  Structural  Chemistry  of  Carbon 

The  element  carbon  is  capable  of  participating  in  a  wide  variety  of  struc¬ 
tures.  In  inorganic  compounds,  carbon  reacts  with  metals  and  semi-metals 
to  form  solid  solutions  and  simple  compounds  (53, 54  ,  55).  Carbon  also  provides 
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Chapter  III 


STRUCTURE  OF  PYROLYTIC  GRAPHITE 


A.  Introduction 

PG  is  one  of  many  forms  of  pure  carbon.  Thus,  the  structural  chemistry 
of  carbon  must  be  understood  before  examining  PG  whose  structure  is  directly 
related  to  the  manner  in  which  it  is  made,  and  whose  properties  depend  upon 
that  structure.  This  structure  can  be  characterized  using  metallographic  and 
x-ray  tecliniques,  with  additional  information  available  through  density  and 
porosity  measurements.  Studies  of  the  structure -sensitive  properties  of  PG 
and  their  change  with  heat -treatment  also  provide  information  concerning  the 
nature  of  the  PG  deposit. 

In  this  chapter,  the  structural  chemistry  of  carlwn  will  be  discussed,  as 
will  the  tecimiques  for  manufacturing  PG.  Following  this,  the  details  of  the 
structural  characterizations  of  PG  using  metallography,  x-ray  diffraction,  and 
other  techniques  will  be  outlined,  and  a  short  discussion  of  graphitization  in 
PG  will  be  given. 

B.  Structural  Ctiemisiry  of  Carlxin 

The  element  carbon  is  capable  of  participating  in  a  wide  variety  of  struc¬ 
tures.  In  inorganic  compounds,  carbon  reacts  with  metals  and  semi-metals 
to  form  solid  solutions  and  simple  compounds  (53,  54,  55).  Carbon  also  provides 
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the  foundation  for  organic  chemistry  by  reacting  with  itself,  hydrogen,  and, 
to  a  lesser  extent,  with  other  elements  to  produce  a  wealth  of  organic  com¬ 
pounds.  These  compounds  are  divided  into  four  general  types:  aliphatic, 
aromatic,  alicyclic,  and  hetero-cyclic  (56),  each  of  which  contains  many  dif¬ 
ferent  structures.  The  aliphatic  compounds  are  based  upon  carbon -bonded 
chains  and  the  aromatics  are  based  upon  six-membered  benzene  rings.  Ali¬ 
cyclic  compounds  are  composed  of  different -sized  carbon -bonded  ring  struc¬ 
tures  with  alicyclic  properties;  hetero-cyclic  compounds  are  ring  compounds 
which  incorporate  other  atoms  —  usually  nitrogen,  oxygen,  and  sulfer  —  into  the 
rings  and  behave  generally  as  either  aromatic  or  alicyclic  compounds.  Carbon 
can  further  combine  with  itself  to  produce  macro-molecules  or  networks  of 
carbon  atoms  which  are  the  structural  units  of  the  solid  carbons .  Solid  carbon 
products  are  divided  into  two  general  categories  (57);  graphitic  carbons  and 
non-graphitic  or  "amorphous"  carbons.  The  graphitic  carbons  consist  of 
approximately  planar  nets  of  carbon  hexagon  rings  which  have  a  tendency, 
upon  heating,  to  produce  the  ideal  hexagonal  graphite  lattice.  Non-graphitic 
or  "amorphous"  carbons  result  when  these  macro-molecules  include  many  dif¬ 
ferent  kinds  of  carbon  bonds.  Tney  are  termed  amorphous  because,  although 
there  is  some  tendency  for  the  macro-molecules  to  lie  parallel  to  each  other, 
there  is  a  very  high  degree  of  disorder  which  produces  a  very  broad,  diffuse 
x-ray  diffraction  pattern. 

The  ideal,  hexagonal  graphite  structure,  first  proposed  by  Bernal  (58)  and 
shown  in  Fig.  1,  consists  of  parallel  sheets  of  carbon  hexagons  arranged  in 
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=  2. 4612  A 

=  6.  7076  A 

1  =  3. 3538  A 

carbon-carbon  separation  =  1.415  A 

4  atoms  per  unit  cell 
X-ray  density  =  2. 267  gm/cm^ 


Fig.  1  Crystal  Structure  of  Ideal  Hexagonal  Graphite 
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an  abab. . .  stacking  sequence.  In  this  structure,  the  interlayer  separation  is 
3.3538  A,  and  the  carbon -carbon  distance  within  each  layer  is  1 .415  A  (9.  59). 
The  carbon-carbon  distance  within  the  sheet  is  virtually  independent  of  inter¬ 
actions  with  other  layers  (60)  because  of  the  bond  structure.  Carbon  has  four 
valence  electrons  per  atom;  three  of  which  are  used  to  form  regulat  covalent 
bonds  with  each  of  the  three  nearest  neighbors  in  the  hexagonal  network,  while 
the  fourth  electron  (called  the  electron)  is  free  to  form  resonance  bonds 
over  the  entire  layer  (32, 59).  These  n  electrons  are  responsible  for  the 
electrical,  thermal,  and  magnetic  properties  of  graphite  and  for  the  interlayer 
bonding  which  permits  the  graphite  structure  to  exist  (9). 

A  slightly  less  stable  form  is  rhombohedral  graphite  which  has  an  abcabc . . . 
stacking  sequence .  The  heat  of  transformation  of  rhombohedral  to  hexagonal 
graphite  has  been  estimated  to  be  -144  ±41  cal/mole  (61).  The  hexagonal 
form  is,  by  far,  the  most  prevalent;  however,  hexagonal  graphite  can  be  par¬ 
tially  converted  into  rhombohedral  graphite  by  grinding  (62).  The  macro- 
molecular  layers  are  considereu  to  be  held  together  by  van  der  Waal's  forces; 
however,  interlayer  charge  transfer  bonds  have  also  been  postulated  (9).  The 
interlayer  bond  energy  for  either  structure  is  estimated  to  be  between  8  and  20 
kcal/moie  depending  upon  the  method  used  to  compute  it  (9).  The  small  dif¬ 
ference  in  energies  for  the  two  stacking  sequences,  the  large  interlayer  separa¬ 
tion,  and  the  rather  low  interlayer  bonding  energy  combine  to  permit  the  partial 
or  complete  lack  of  a  layer  stacking  sequence.  This  lack  of  order  is  known  as 
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"turbostratic."  disorder  (63).  Its  principal  effect,  to  be  discussed  later,  is 
to  increase  the  interlayer  spacing.  When  the  spacing  is  greater  than  the  ideal 
3.3538  A,  the  structure  is  referred  to  as  nearly  ideal  hexagonal  graphite  (9). 

Carbon  can  exist  in  still  another  form,  the  diamond  structure,  in  which 
the  atoms  are  situated  symmetrically  with  four  equal  tetrahedral  bonds  to 
nearest  neighbors.  The  diamond  structure  is  unstable  relative  to  the  hexagonal 
graphite  structure,  with  the  heat  of  transformation  at  298  K  estimated  to  be 
453  ±20  cal/mole  (9). 

Graphite  products  are  available  in  many  forms.  Their  differences  result 
from  the  relative  proportions  of  "amorphous"  carbon  and  nearly  ideal  hexagonal 
graphite.  Natural  graphite  consists  principally  of  large  crystallites  with  per¬ 
fectly  stacked  layers,  while  carbon  blacks,  soots,  and  hard-baked  carbons  con¬ 
sist  of  randomly  bonded  carbon  macro-molecules,  with  some  regions  of  small 
stacks  of  nearly  parallel  hexagon  sheets.  For  reviews  and  listings  of  some  of 
the  many  different  forms  of  carbon  ‘.ee  references  27,  59,  64,  65,  66. 

Synthetic  graphites  of  widely  variable  properties  are  produced  by  heating 
organic  cokes  with  suitable  binders  to  high  temperature  in  the  absence  of  air 
to  drive  off  the  non-carbon  atoms  and  produce  the  various  macro-molecules  of 
carbon  described  above.  The  resulting  body  contains  crystallites  of  approxi¬ 
mately  ideal  graphite  surrounded  by  somewhat  less  organized  carbon.  This 
material  is  relatively  porous  and  soft,  rubbing  off  onto  other  objects  easily 
(the  reason  why  graphite-clay  mixtures  are  used  in  pencils).  The  less 
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graphitized  carbon  products  do  not  mark  other  objects  as  readily,  because  they 
contain  smaller  volume  fractions  of  the  soft,  nearly  ideal,  graphite  particles. 

Synthetic  carbon  products  can  also  be  produced  by  the  pyrolytic  decompo¬ 
sition  of  simple  organic  compounds  like  methane  under  controlled  pressure  and 
temperature  conditions  (1,67,68,69).  The  resulting  bodies  differ  from  the 
other  graphites  m  that  there  is  very  little  truly  disorganized  carbon  and  the 
crystallites  tend  to  lie  parallel  to  the  deposition  substrate,  producing  a  high 
degree  of  anisotropy  in  many  properties.  The  density  of  pyrolytic  graphite 
is  a  function  of  the  deposition  temperature,  being  higher  the  higher  the  tem¬ 
perature.  Deposits  made  above  about  2100°C  (1, 67)  approach  the  theoretical 
x-ray  density. 

All  tlie  as-deposited  PC's  have  a  high  degree  of  turbostratic  disorder, 
leading  some  authors  to  term  this  material  pyrolytic  carbon  (PC).  There  is 
no  definite  distinction  between  PG  and  PC,  but  the  less  ordered  material  with 
very  small  crystallite  size  produced  by  deposition  below  about  1900°C  should  be 
called  PC  and  the  material  with  the  larger  crystallite  size  and  nearly  theoreti¬ 
cal  density  deposited  above  1900®C  sh'-Jld  be  termed  PG. 

In  contrast  with  the  other  forms  of  synthetic  graphite,  PG  is  hard  and 
dense,  and  pieces  of  PG  a  few  inches  square  produce  a  "metallic  ring"  when 
dropped  on  a  hard  surface. 
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C.  Manufacture  of  PG 


Methods  of  making  PG  were  described  as  early  as  1880  by  Sawyer  (70) 
and  1883  by  Edison  (12).  The  technique  involves  thermally  decomposing  a 
carbon -bearing  gas  at  or  near  a  substrate,  causing  carbon  to  be  deposited 
on  the  substrate.  Until  the  mid-1950’s  the  size  of  deposits  was  limited, 
and  the  only  use  for  carbon  made  by  this  process  was  the  carbon  electrical 
resistor  (71).  In  1958,  Brown  and  Watt  (1)  gave  the  first  report  on  the  pro¬ 
duction  of  massive  PG.  Their  work  was  followed  closely  by  Diefendorf  (67, 
68,69)  of  General  Electric,  and  Pappis  and  Blum  (3)  of  Raytheon  Company. 

PG  technology  advanced  very  rapidly,  and  by  1963  at  least  four  companies 
were  manufacturing  many  different  forms  (66). 

PG  is  deposited  at  low  gas  pressure  with  substrate  temperatures  between 
800  and  28(X)°C  (12, 67, 68, 69).  For  structural  use,  it  is  generally  deposited 
above  2100°C  to  provide  high  density  and  good  mechanical  properties.  TVo 
methods  of  heating  the  substrate  are  used:  self-resistance  heating  and  induc¬ 
tion  heating  (67,  72).  Two  other  variations  in  furnace  design  are  the  use  of 
cold  wall  and  hot  wall  furnaces  (67).  The  cold  wall  furnace  generally  produces 
non-homogeneous  material  because  of  large  temperature  gradients  through  the 
thickness  of  the  PG  and  large  gradients  in  composition  and  temperature  in  the 
gas  (67).  As  a  result,  most  commercial  facilities  use  hot  wall  furnaces. 

The  important  parameters  controlling  the  structure  of  the  deposit  are  the 
deposition  temperature  and  gas  pressure.  Both  of  these  control  the  rates  at 
which  nucleation  and  growth  occur  at  the  deposition  surface,  and  these  rates 
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control  the  degree  of  crystalline  perfection  (72),  the  degree  of  preferred 
orientation  (73),  the  crystallite  size  (73),  and  many  other  structural  parameters 
(67). 

D .  Wetallographic  Structure  of  PG 

In  physical  appearance,  PG  is  dark  gray  and  can  be  polished  to  a  dark 
metallic  luster.  'I’he  surface  away  from  the  substrate,  the  terminal  deposition 
surface,  is  covered  with  slightly  rounded  nodules  of  a  size  characteristic  of 
the  deposition  process.  The  substrate  surface  (the  surface  at  which  nuclea- 
tion  first  started)  is  covered  with  much  smaller  circular  dimples  which  reveal, 
approximately,  the  original  nucleation  sites.  The  structure  of  PG  is  formed 
by  the  competitive  nucleation  and  growth  of  cones  or  paraboloids  with  their 
major  axes  peipendicular  to  the  deposition  surface.  Two  major  types  of  cone 
growth  are  recognized  (12),  substrate  nucleated  and  regeneratively  nucleated. 
Substrate  nucleated  PG  consists  principally  of  cones  that  were  nucleated  at  the 
substrate  and  continued  to  grow  throughout  the  thickness  of  the  deposit,  while 
regeneratively  nucleated  PG  is  characterized  by  a  recurring  nucleation  of  cones 
through  the  deposit.  Material  representative  of  both  structures  was  used  in 
this  study,  and  photomicrographs  of  the  different  PC’s  are  shov/n  and  discussed 
in  Chapter  VII.  In  the  substrate  nucleated  material,  other  cones  or  paraboloids 
will  occasionally  nucieate  in  the  middle  of  the  deposit  and  grow  in  a  fashion 
similar  to  the  substrate  nucleated  cones .  The  presence  of  soot  or  other  foreign 
matter  on  the  advancing  deposition  surfaces  causes  this  nucleation  (67)  which 
results  in  the  extra  cones  and  in  the  regeneratively  nucleated  structure. 


Actually,  the  structure  of  PG  deposits  ranges  from  that  in  which  the  substrate 
nucleated  cones  are  dominant  to  that  dominated  by  the  regeneratively  nucleated 
cones. 

PG  microsiructure  is  further  complicated  by  the  presence  of  a  fine, 
secondary  substructure  that  is  presumably  associated  with  many  inefficient 
nucleation  sites  on  the  advancing  deposition  surface.  Pappis  (74)  has  charac¬ 
terized  the  PG  microstructure  according  to  the  following  scheme: 

Primary  Structure  Dominant  (substrate  nucleated) 

(A)  Straight  sided  cones —Al,  A2,  A3* 

(B)  Curved  sided  cones  -  Bl,  B2,  B3 

Secondary  Structure  Dominant  (combination  substrate  and  regenera¬ 
tive  nucleated) 

(C)  Straight  sided  cones  -C2,  C3  (no  Cl  observed) 

(D)  Curved  sided  cones —Dl,  D2,  D3 

Intra-layer  Structure  (regeneratively  nucleated) 

(E)  Straight  sided  cones  —  El,  E2,  E3 

(F)  Curved  sided  cones  —  FI,  F2,  F3 

and  this  scheme  was  employed  to  characterize  the  material  used  in  this  investi¬ 
gation.  The  anisotropy  of  PG  is  reflected  in  its  optical  behavior  under  cioss- 
polarized  light;  photomicrographs  are  generally  taken  at  15  degrees  from 


*1,2,3  correspond  to  coarse,  medium,  and  fine  cone  sizes,  respectively. 
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extinction  to  provide  maximum  resolution  of  the  structure.  The  fact  that  the 

cone  structure  is  clearly  observed  under  polarized  light  is  proof  that  the 

boundaries  seen  in  the  microscope  represent  the  boundaries  between  groups 

of  crystallites  having  different  orientations;  however,  within  the  major  cones 

and  within  the  secondary  cones,  there  is  a  uniform  variation  in  the  optical 

> 

activity  with  sample  rotation,  demonstrating  that  the  crystallographic  orienta¬ 
tion  within  even  the  minor  cones  varies  from  one  side  of  the  cone  to  the  other. 

High  magnification  transmission  electron  microscopy  studies  (75,76,77) 
have  shown  that  PG  consists  of  stacks  of  overlapping  wrinkled  sheets  a  few 
hundred  A  thick,  with  facets  on  the  wrinkled  sheets  about  200  A  wide.  Dis¬ 
location  structures  such  as  the  networks  observed  in  graphite  single  crystals 
(78)  and  synthetic  graphite  (79)  have  not  been  observed  in  PG.  As  will  be  dis¬ 
cussed  in  the  next  section,  the  crystallite  diameter  of  PG  is  about  200  A  while 
that  of  other  graphites  is  about  3000  A.  Because  it  is  difficult  to  re^solve  struc¬ 
tural  details  within  a  200  A  diameter  region,  no  dislocation  structure  has  been 
observed  in  PG. 

E  .  X-ray  Structure  of  PG 

In  this  section,  the  structure  of  PG,  as  determined  by  x-ray  diffraction, 
will  be  discussed,  using  the  interpretations  of  diffraction  patterns  developed 
for  carbon  blacks  and  graphitic  carbons.  First,  the  layer  structure  and  layer 
stacking  will  be  discussed,  followed  by  discussions  of  crystallite  size  and  pre¬ 
ferred  orientation.  Next,  the  variation  of  structure  with  orientation  will  be 
noted  and  a  discussion  of  small  angle  scattering  in  PG  will  be  given. 
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Brown  and  Watt  (1)  first  demonstrated  that  PG  has  an  x-ray  structure 
similar  to  that  of  graphitic  carbon  so  the  detailed  theories  developed  for  these 
carbon  materials  can  be  applied  to  PG.  For  instance,  Warren  (80)  showed  that 
graphitic  carbons  are  made  up  of  small  regions  of  approximately  parallel 
layer  planes  which,  in  turn,  are  composed  of  hexagonal  benzene  rings,  just 
like  the  graphite  structure.  Warren  (81)  further  demonstrated  that  such 
carbons  show  only  two  kinds  of  diffraction  peaks:  0001  reflections  from  the 
parallel  layer  planes,  and  two-dimensional  1^  bands  which  are  the  result  of 
a  relaxation  of  the  crystalline  order  in  the  c -direction.  This  concept  was 
further  developed  by  Biscoe  and  Warren  (63)  to  prove  that  carbon  blacks  con¬ 
sist  of  graphite-like  layers  stacked  in  parallel  groups  but  otherwise  without 
mutual  orientation.  The  concept  of  a  stacking  order  in  carbon  black  has  no 
meaning  because  the  layers  are  randomly  rotated  about  their  mutual  layer 
plane  normal.  As  a  result  of  this  structure,  the  OOOi  diffraction  peaks  are 
symmetric,  while  the  hk  bands  are  asymmetric  with  the  intensity  falling  off 
very  steeply  on  the  low  26  side  and  decreasing  very  gradually  on  the  high 
26  side. 

The  concept  of  a  crystallite  size  in  carbon  blacks  is  still  valid,  being  a 
measure  of  the  range  over  which  the  layer  planes  are  parallel  (L^)  and  the 
average  diameter  of  the  layers  in  the  stack  (L^).  This  will  be  discussed  in 
more  detail  shortly. 

Carbon  blacks  are  grouped  in  two  general  classes  (57, 82),  those  which 
approach  the  graphite  layer  stacking  order  upon  heat  treatment  and  those  that 
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do  not.  Carbons  in  the  first  class  are  called  graphitizing  and  those  in  the 
second,  ncn -graphitizing.  When  graphitizing  carbons  have  become  partially 
graphitized,  the  1^  bands  show  modulations  which  indicate  the  beginning  of 
diffuse  1^  reflections.  At  the  same  time,  the  interlayer  spacing  computed 
from  the  OOOf  diffraction  lines  is  found  to  decrease,  following  the  decrease 
in  the  degree  of  disorder  in  the  blacks.  This  observation  has  been  treated 
by  Franklin  (82),  Bacon  (83),  Bacon  and  Franklin  (84)  and  Bacon  alone  (60) 
in  some  detail  with  the  conclusion  that  four  different  interlayer  spacings  exist 
in  a  stack  of  graphitic  layers.  A  description  of  these  layers  and  their  spacings 
is  given  below: 

1.  Between  two  oriented  layers;  e.g.,  with  the  corners  of  the 

hexagons  in  one  layer  directly  above  the  centers  of  the  hexagons 
of  the  other  (d  =  3.354  A) 

2.  Between  two  disoriented  layers,  each  of  which  is  disoriented  with 

respect  to  its  neighbor  (d  =  3.440  A) 

3.  Between  two  disoriented  layers,  one  of  which  is  oriented  with 

respect  to  its  neighbor  (d  -  3.408  A) 

4.  Between  two  disoriented  layers,  both  of  which  are  oriented 

with  respect  to  their  neighbors  (d  =  3.376  A) 

Franklin  (82)  and  Bacon  (60)  were  able  to  show  that  the  probability  of  a  stacking 
error,  p  ,  as  determined  from  analyses  of  the  diffraction  breadths  and  shapes 
of  hkl  diffraction  peaks,  could  be  related  to  the  observed  interlayer  spacing. 
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^0002  ’  relationship  is  shown  in  Fig.  2.  Various  theoretical  expressions 

were  derived  to  fit  the  observed  curve  (60)  and  they  agree  fairly  well  over  cer¬ 
tain  ranges  of  p  .  Although  a  direct  comparison  of  p  versus  dQQ2  PG  has  not 
been  made,  it  is  reasonable  to  assume  that  the  same  relations  hold,  provided 
accurate  determinations  of  c1qqq2  are  made. 

The  a  parameter  of  PG  can  be  obtained  from  the  1£  and  JLl  bands,  using 
the  hexagonal  d-spacing  relationship 


l/d^  «  4/3a2  (h^  +  Ic2  +  hk)  -f  l^/c^ 


Warren  (81)  has  derived  a  correction  factor  for  the  peak  shift  of  an  hk  band 
from  the  theoretical  position  of  the  corresponding  hko  peak.  The  shift  is  to 
higher  angles  and  of  magnitude 

A  sin  e  =  0.16  X/L_ 

d 

For  PG,  is  typically  greater  than  100  A,  making  this  shift  negligibly  small. 
A  precision  extrapolation  is  possible  if  several  hk  bands  are  observed  but  the 
data  invariably  yield  2,46  A,  regardless  of  the  condition  of  the  PG,  indicating 
the  relative  stability  of  the  intra -planar  bonds  over  the  inter -planar  bondk. 

Estimates  of  the  crystallite  sizes,  and  ,  of  graphitic  carbons  .  ng 
the  Scherrer  equation,  with  appropriate  modification  (60,  85, 86)  or  using  a 
Fourier  analysis  of  the  diffraction  peaks  (72, 86, 87, 88)  have  been  made.  The 
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STACKING  ERROR  PROBABIUTY  l(p) 


Fig.  2  Interlayer  Spacing  vs.  Stacking  Error  Probability  for  Graphitic  Carbons 


difficulty  of  estimating  crystallite  thickness  lies  in  establishing  the  definition 
of  a  crystallite  in  graphitic  carbon  and  interpreting  the  effect  of  the  four  inter^ 
planar  spacings  postulated  by  Bacon  and  Franklin  on  the  line  breadth  and  shape . 
Houska  and  Warren  (88)  have  attempted  to  account  for  the  effect  of  two  layer 
spacings  in  carbon  black  and  Guentert  and  Cvikevich  (72,86)  have  applied  the 
rigorous  Warren -Averbach  method  (89)  to  PG.  Guentert  and  Cvikevich  (GC) 
(86)  reported  their  data  with  sufficient  detail  to  make  possible  a  comparison  of 
their  Fourier  results  with  those  from  the  more  convenient  Scherrer  equation. 
Table  5  shows  the  results  of  this  comparison.  For  crystallite  thicknesses 
below  about  700  A,  the  0002  line  breadths  with  the  Warren  correction  (90) 


Table  5 

COMPARISON  OF  METHODS  OF  COMPUTING 
CRYSTALUTE  THICKNESS  IN  PG 


Specimen  ' 
Number 

Crystallite 
Thickness -L^  „  . 
(Fourier  Analysis)*'^' 

Crystallite  Thickness  -  l^, 
(Warren  Correction)  (^) 

0002  0004  (c)  0006  (c) 

AD 

200  A 

190  A 

100  A 

70  A 

25 

310 

310 

180 

120 

27 

380 

380 

230 

170 

28 

700 

600 

350 

300 

30 

850 

1120 

680 

630 

(a)  Specimen  designation  and  Fourier  analysis  by 
Guentert  and  Cvikevich  (Ref.  86) 

(b)  =  b2  -  b^  ,  Lc  =  A//3  cos  9 

(c)  Different  orders  from  (0001)  planes 
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give  fairly  good  agreement  with  the  more  rigorous  Fourier  treatment.  GC  ^6) 
further  showed  that  the  line  broadening  of  the  bands  in  PG  is  due  only  to 
particle  size  so  the  Warren  modification  (81)  of  the  Scherrer  equation 

I  =  1-84X 

a  jS  cos  6 

is  directly  applicable:  however,  the  peak  should  be  corrected  (86),  for  instru¬ 
mental  broadening  (90),  ^1^2  broadening  (911,  preferred  orientation 

effects  (86,92),  and  absorption  broadening  (93).  Although  a  correction  factor 
for  the  hk  bands  for  preferred  orientation  effects  has  been  derived  by  GC  (86, 
92),  they  did  not  apply  it  to  their  earlier  (86)  data  because  they  estimated  the 
effect  to  be  less  than  20%.  The  effect  of  the  preferred  orientation  correction 
IS  to  increase  the  line  width  while  the  absorption  correction  tends  to  decrease 
it;  thus,  if  the  corrections  are  ignored,  their  effects  tend  to  be  cancelled. 
Partial  ordering  leads  to  modulations  in  the  hk  bands  of  PG  and  carbons  (72, 
86,87),  and  line  broadening  theories  for  this  condition  have  not  been  fully 
developed  (86) . 

The  layer  structure  of  graphite  causes  it  to  exhibit  a  high  degree  of  -re¬ 
ferred  orientation  when  pyrolytically  deposited  (67).  The  preferred  orienta¬ 
tion  in  PG  has  been  studied  by  GC  (72,  86,87),  Bragg  with  Hammond  and 
co-workers  (5,  94),  and  Noda,  et  al.  (73).  Because  they  find  it  convenient 
to  use  in  their  calculations,  Noda  and  Guentert  represent  the  preferred  orienta¬ 
tion  with  a  cos'^o;  function,  where  a  is  the  inclination  angle  of  the  crystallite 
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layer  planes  relative  to  the  deposition  surface.  On  the  other  hand,  Bragg  and 
Packer  (94)  found  that  the  preferred  orientation  had  a  Gaussian  distribution 
over  two  orders  of  magnitude  in  pole  density,  and  they  represent  their  data 
as  an  exponential  function  of  a  .  Both  functions  are  approximations, 
with  the  cos*^Q  function  showing  deviations  at  a  layer  plane  pole  Jensity 
of  about  10%  (86,92)  while  the  Gaussian  function  behaves  well  out  to  1%  pole 
density  in  as-deposited  material  (94)  and  out  to  about  5%  in  heat  treated  mate¬ 
rial  (5).  The  effects  of  heat  treatment  will  be  discussed  in  more  detail  later. 

Bragg  and  Packer  (94)  also  demonstrated  that  the  structure  of  PG  is  orien¬ 
tation  dependent,  with  the  degree  of  perfection  generally  decreasing  with  increas¬ 
ing  inclination  to  the  deposition  surface  and  with  the  disparity  in  perfection  in¬ 
creasing  markedly  with  heat  treatment.  This  change  in  degree  of  perfection  re¬ 
quires  the  use  of  integrated  intensities  for  all  preferred  orientation  studies  in  PG . 

Small  angle  scattering  by  PG  has  been  studied  by  Bragg,  with  Hammond  and 
others  (95).  They  showed  that  PG  contains  disk-like  voids,  22  A  thick  and  44  A 
in  diameter,  that  lie  in  the  deposition  plane  and  have  a  spacing  comparable  to 
that  of  the  crystallites .  Such  voids  in  PG  could  very  easily  be  associated  with 
the  wrinkled  sheets  found  in  Stover's  electron  microscopy  studies  (75). 

F.  Other  Structural  Parameters 

In  addition  to  the  parameters  discussed  above,  the  PG  structure  can  be 
further  characterized  by  density  and  pore  spectrum  determinations.  It  has 
been  noted  that  PG  is  impermeable  to  inert  gases  (67)  and  liquids  (96),  in 
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agreement  with  the  lack  of  visible  porosity  in  the  microstructure.  This  means 
that  the  pore  spectrum  cannot  be  determined  by  gas  or  liquid  pwrosimetry,  as 
is  done  with  other  synthetic  graphites  (27),  but  must  be  measured  using  small 
angle  scattering  of  x-rays.  The  displaced-volumo  density  is  slightly  less  than 
tnat  computed  from  the  x-ray  layer  spacing,  indicating  a  small  sub-microscopic 
porosity.  The  x-ray  density  of  PG  is  a  function  of  the  unit  cell  dimensions  and, 
since  the  a  parameter  is  essentially  constant  regardless  of  the  condition  of  the 
graphite,  the  change  in  x-ray  density  is  related  only  to  the  changes  in  the  layer 
spacing.  With  an  a  parameter  of  2.461  A,  the  x-ray  density  of  nearly  ideal 
hexagonal  graphite  can  be  expressed  as 


15.207  ^  7.603 
-  ^^0002 


3 

where  is  the  density  in  gm/cm  ,  c  the  unit  cell  height  in  A,  and  dQQQ2 
the  average  layer  spacing  in  A.  Table  6  lists  the  graphite  x-ray  densities  for 
various  values  of  dQQQ2  and  gives  some  experimental  values  of  dQQQ2 
for  comparison.  The  measured  density  is  about  0.02  gm/cm  or  1%  below  the 
theoretical  value  and  this  difference  must  be  associated  with  the  porosity  indi¬ 
cated  by  the  small  angle  scattering. 

Information  regarding  the  structure  of  PG  can  also  be  obtained  by  measur¬ 
ing  structure-sensitive  properties,  such  as  transport  and  mechanical  properties. 
The  thermal  (97)  and  electrical  (6, 98)  conduction  properties  have  been  used  to 
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Table  6 

X-RAY  DENSITY  AS  A  FUNCTION  OF  INTERLAYER  SPACING 
AND  COMPARISON  WITH  EXPERIMENT 


/q  \ 

Computed  Observed'  ’ 


*^0002 

3 

gm/cm 

*^0002 

^  3 

gm/cm 

3.440  A 

2.210 

3.420 

2.223 

3.420 

2.200 

3.400 

2.236 

3.380 

2.249 

3.370 

2.256 

3.372 

2.225 

3.360 

2.263 

3.366 

2.238 

3.354 

2.267 

3.360 

2.254 

(a)  Observed  data  from  Pappis  (Ref.  74  ) 


demonstrate  the  anisotropy  and  preferred  orientation  of  PG  and  to  gain  an  under¬ 
standing  of  its  band  structure.  In  addition,  mechanical  property  measurements 
at  high  temperature  (5, 77)  have  been  used  to  determine  the  nature  of  the  bonding 
between  the  crystallites  in  PG.  One  of  the  purposes  of  this  investigation  is  to 
employ  the  mass  transport  properties  of  PG  to  increase  our  understanding  of 
PG  structure . 

G.  Graphitization  in  PG 

Graphitization  is  the  general  process  by  which,  on  heat-treatment,  carbon 
products  approach  the  ideal  graphite  structure.  The  intermediate  stages  of 
this  process,  in  the  case  of  PG.  arc  indicated  by  the  following  changes  (74,86,99): 

•  An  elimination  of  the  secondarj'  microslructure 
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•  A  decrease  in  the  interlayer  spacing 

•  A  sharpened  layer  plane  texture 

•  An  increase  in  crystallite  size  ^ 

•  A  shrinkage  perpendicular  to  the  deposition  plane 

•  A  slight  elongation  parallel  to  the  deposition  plane 

•  An  increase  in  density 

For  a  given  thermal  treatment,  those  crystallites  oriented  nearly  parallel  to 
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the  deposition  surface  graphitize  more  fully  than  those  inclined  to  the  surface 
(5, 94).  Furthermore,  the  anisotropy  ratios  of  transport  properties,  i.e.,  the 
electrical  and  thermal  properties,  are  increased  by  graphitization  (74).  The 
terminal  stages  of  the  process  are  noted  by  a  nearly  complete  elimination  of 
the  microstructural  detail  aitd  a  closer  approach  to  the  theoretical  values  of 
the  various  x-ray  parameters.  Perfect  single  crystals  have  not  been  produced 
by  graphitization,  but  material  approaching  a  mosaic  single  crystal,  as  deter¬ 
mined  by  the  interlayer  spacing,  crystallite  size,  and  texture,  have  been 
produced  (86).  One  problem  with  such  material  is  the  presence  of  delamina¬ 
tions  parallel  to  the  deposition  surface  which  are  introduced  when  the  texture 
is  sharpened. 

The  rate  of  graphitization  in  PG  becomes  appreciable  around  2500°C  (5) 
and  reaches  terminal  stages  in  reasonable  times  somewhat  above  3000^0  (99). 
Concurrent  stress,  either  in  tension  parallel  to  the  deposition  plane  (5,77)  or 
in  compression  (lOO)perpendiculartoit,  considerably  enhances  the  graphitiza¬ 
tion  process,  lowering  the  graphitization  temperature  50  to  100°C  and  greatly 


greater  perfection. 


V 

Kinetics  of  the  above  changes  have  not  been  studied  in  detail;  instead, 
the  subject  has  been  treated  much  as  the  early  studies  of  recrystallization  in 
metals.  In  these  studies,  changes  in  various  properties  were  observed  follow¬ 
ing  fixed-time  anneals  at  different  temperatures.  From  a  plot  of  property 
versus  temperature  of  anneal,  a  recrystallization  or  graphitization  tempera¬ 


ture  can  be  defined  as  the  temperature  at  which  a  large  change  in  property 
occurs.  One  study  of  the  kinetics  of  graphitization  has  been  made  by  Fischbach 


(101)  who  used  magnetic  susceptibility  and  d 
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measurements  to  characterize 


the  PG  structure.  The  rates  of  change  of  these  properties  with  time  at  tem¬ 


perature  were  analyzed  using  sequential  first  order  reaction  rate  theory,  but 


no  justification  of  this  procedure  was  given.  It  is  by  no  means  clear  that 


these  properties  bear  a  homogeneous  relationship  to  the  amount  of  transfor¬ 
mation  as  is  required  to  apply  reaction  rate  theories  based  on  the  Law  of 
Mass  Action  (102). 


H.  Summary 

Studying  the  effects  of  PG  structure  on  a  property  requires  a  complete 
characterization  of  the  structure.  This  characterization  involves  the  simul¬ 
taneous  application  of  several  techniques  of  metallography  and  x-ray  diffrac¬ 
tion,  together  with  accurate  density  measurements.  These  techniques  have 
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been  Applied  to  the  8tructur«l  aepects  of  mechanical  and  electrical  properties 

- 

and  will  be  required  to  examine  the  role  of  structure  on  metal  diffusion  in  PG . 
Some  details  of  these  techniques  and  their  general  results  have  been  given  in 
the  discussion  above . 


K 
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Chapter  IV 

EXPERIMENTAL  APPROACH 

A .  Statement  of  Problem 
It  is  evident  from  the  discussion  in  Chapters  I  and  II  that  a  study  of  the  role 

of  structure  in  metal  diffusion  in  graphites  would  be  both  timely  and  useful .  PG 
is  the  least  complex  massive  graphite  available,  and  techniques  for  characteriz¬ 
ing  its  structure  have  been  developed  to  the  point  where  comparisons  of  structure 
and  property  are  possible.  A  study  of  diffusion  in  PG  requires  developing  the 
necessary  experimental  techniques  and  analytical  methods  to  treat  the  problem 
and  relating  the  results  of  the  diffusion  studies  to  the  structure  of  the  material 
being  examined . 

The  purposes  of  this  investigation  are  to  determine  the  fundamental  mechan¬ 
isms  of  metal  diffusion  in  PG,  to  increase  the  understanding  of  the  structure  of 
PG  through  these  studies,  and  to  investigate  the  possibility  of  diffusion -alloying  PG 

B.  Program  Summary 

Because  of  the  expected  low  solubility  of  metals  in  graphite,  high  sensitivity 
radioactive  tracer  methods  are  ideally  suited  for  studying  metal  diffusion  in  PG . 
Precision  lapping  techniques  are  also  well  suited  to  this  study  because  they  are 
sensitive  and  readily  adapted  to  radioactive  tracer  methods. 


rm 
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To  accomplish  the  objectives  mentioned  above,  three  different  types  of  PG 
were  chosen  and  specimens  from  each  were  given  a  single  heat -treatment  so  that 
both  the  as-deposited  and  partially  graphitized  states  could  be  studied.  Radioiso¬ 
topes  of  Re  and  Nb  were  selected  as  the  diffusing  species,  with  the  emphasis 
placed  on  the  Nb.  Diffusion  reactions  of  Re  with  PG  were  examined,  and,  with 
a  precision  lapping  device,  the  diffusion  of  Nb  in  all  three  types  of  PG  was  studied. 
The  diffusion  coefficients  of  Nb  in  PG  were  determined  for  two  conditions  of  each 
PG  (as -deposited  and  heat-treated)  and  for  the  different  "crystallographic"  direc¬ 
tions  in  PG  (parallel  and  perpendicular  to  the  deposition  surface).  Autoradiography 
was  employed  to  observe  any  non -homogenous  distribution  of  the  diffusing  atoms. 
The  temperature  range  from  1400  to  1950°C  was  covered,  corresponding  to  a 
range  of  three  decades  in  diffusivity  for  each  direction.  The  structure  of  each 
PG  was  characterized  by  the  following: 

•  Metallographic  structure  at  high  and  low  magnification 

•  Transmission  electron  microscopy 

•  Stacking  error  probability,  p  ,  or  degree  of  graphitization 

•  Crystallite  size,  and 

•  Degree  of  preferred  orientation,  as  measured  by  the  mean  tilt  angle 

•  Orientation  dependence  of  the  structural  parameters  p  and  L 

c 

•  Submicroscopic  void  size  and  distribution,  determined  from  small 
angle  scattering 

•  High  pressure  Hg  porosimetry 

•  Density,  as  determined  by  hydrostatic  weighing 


A  mathematical  analysis  of  diffusion  in  layered  structures  was  made,  and 
correlations  between  the  diffusion  and  structural  parameters  were  noted . 

C .  Choice  of  PG 

Three  different  types  of  PG  were  chosen  to  be  representative  of  the  vari^j 
ous  types  of  PG  available.  The  three  Qrpes  were  all  produced  by  Raytheon^JH 
Company  during  their  development  work  on  PG  and  may  be  described  in  the  | 
Pappis  notation  listed  in  Chapter  II  (74)  as  follows:  I 

R277:  A2  structure,  singularly  nucleated,  medium  cone  size 

R16S:  C3  structure,  combination  singularly  nucleated  and 

regeneratively  nucleated,  fine  cone  size  | 

Rl:  FI  structure,  heavily  regeneratively  nucleated,  medium 

cone  size  | 

Most  of  the  work  was  concentrated  on  the  R16S  material  because  it  had  the 
fewest  delaminations  and  most  uniform  structure.  The  R277  and  Rl  materials 
were  marked  by  a  variation  in  structure  from  one  place  in  the  deposit  to  another 
and  a  tendency  toward  delamination. . 

Since  this  investigation  began,  there  have  been  many  improvements  in  the 
production  of  the  PG  leading  to  thicker  and  more  homogeneous  deposits  with  less 
tendency  toward  delamination.  However,  the  conclusions  reached  by  this  investi¬ 
gation  will  also  apply  to  the  new  material  and  this  will  be  discussed  in  Chapter  IX 


D.  Choice  of  Metal 


||  Initially,  it  was  thought  that  diffusion  of  metals  in  PG  would  occur  with 
appreciable  rate  only  at  very  high  temperatures  due  to  the  high  sublimation 
temperature  of  graphite.  Because  of  this,  metals  with  low  vapor  pressure  at 
high  temperature  were  sought  so  that  the  boundary  conditions  associated  with 
the  diffusion  system  would  not  be  complicated  by  vapor  transport.  Because  the 
vapor  pressure  of  a  metal  is  generally  lower  the  higher  its  melting  tempera¬ 
ture,  only  high  melting  point  metals  were  considered  for  these  studies. 

A  second  criterion  for  choosing  metals  for  the  diffusion  studies  in  PG  was 
whether  a  carbide  forms.  The  description  of  a  diffusion  system  that  includes 
the  formation  and  growth  of  compounds  is  not  as  simple  as  the  description  of  a 
system  without  compound  formation.  Because  of  this,  carbide  formers  should 
be  used  with  caution  in  studies  of  diffusion  in  graphites . 

The  final  criterion  for  choosing  metals  was  the  availability  of  a  convenient 
radioisotope  because  radioactive  tracer  methods  were  to  be  used  to  determine 
compositions . 

With  above  considerations  in  mind,  the  higher  melting  metals  were  divided 
into  three  groups: 

(1)  Mono-carbide  formers;  Ti,  Zr,  Hf 

(2)  Multiple -carbide  formers;  V,  ]^,  T^,  Cr,  I^,  W,  Fe, 

Co,  Th,  U 

(3)  Non-carbide  formers;  Rh,  Ir,  Pd,  Pt 

Those  with  melting  temperatures  in  excess  of  2500°C  are  underlined. 


Re  was  first  selected  based  on  a  lack  of  reported  carbides  (103),  a  low  vapor 
pressure  (54),  and  a  readily  available  radioisotope  (104).  Re-186  has  a  half-life 
of  3.87  days,  decays  principally  by  emission,  and  is  available  as  HReO^  in 
HNO^  from  the  Oak  Ridge  Isotope  Development  Center  (104) with  a  relative  speci¬ 
fic  activity  of  about  1%.  As  discussed  in  die  next  chapter,  Re-186  was  not  satis¬ 
factory  for  these  diffusion  studies,  and  another  metal  with  a  very  high  specific 
activity  had  to  be  selected . 

Nb-95  and  Co-58  were  the  only  readily  available  radioisotopes  with  carrier- 
free  specific  activity.  Nb  was  selected  over  Co  because  Nb  has  a  lower  vapor 
pressure  at  high  temperature  and  is  a  by-product  of  some  nuclear  reactors.  The 
desirability  of  having  a  non-carbide  former  was  considered  to  be  secondary  in 
this  case  because  the  carrier-free  specific  activity  would  permit  ttie  use  of  a 
very  small  quantity  of  metal  in  the  diffusion  couples . 

Nb-95  has  two  isomeric  states.  The  upper  isomer  decays  by  internal  con¬ 
version  with  a  0.23  Mev  electron  and  a  90  hour  half-life.  The  lower  isomer,  on 
which  these  experiments  were  based,  decays  to  Mo-95  by  emitting  a  0.16  Mev 
and  a  0.75  Mev  y-ray.  The  half-life  of  this  isomer  is  35  days  (104). 

E .  Summary 

Re-186  and  Nb-95  were  selected  to  study  the  diffusion  of  metals  in  PG. 
Because  of  difficulties  related  to  specific  activity,  Nb-95  was  the  only  metal  for 
which  diffusion  coefficients  were  determined.  The  diffusion  coefficients  of  Nb-95 
were  measured  as  a  function  of  temperature  and  "crystallographic"  direction  in 


51 


three  different  kinds  of  PG  in  two  different  conditions.  The  structure  of  each 
PG  was  characterized  using  x-ray  and  metallographic  techniques,  and  correla¬ 
tions  between  structure  and  diffusion  data  were  sought  to  determine  the  funda¬ 
mental  mechanisms  of  metal  diffusion  in  PG. 


EXPERIMENTS  WITH  RHENIUM  AND  PYROLYTIC  GRAPHITE 


A .  Introduction 

Initial  experiments  were  concerned  with  high  temperature  phase  relations 
in  the  Re-PG  system,  using  massive  Re-PG  diffusion  couples.  Following  this 


work,  studies  of  the  diffusion  of  radioactive  Re  in  PG  were  started,  but  these 


studies  had  to  be  abandoned  because  of  difficulties  related  to  the  specific  activity 
of  the  Re  isotope.  A  description  of  the  observed  phase  reactions  is  given  below 


as  is  a  summary  of  the  diffusion  studies. 


B.  Phase  Reactions  of  Re  with  PG 

Diffusion  couples  consisting  of  1/8-in.  diameter  rods  of  Re  and  PG  or  other 
carbon  product  were  clamped  in  a  graphite  fixture  and  annealed  in  argon  for  the 
various  times  and  temperatures  given  in  Table  7.  Following  the  diffusion  anneals, 
the  specimens  were  mounted  parallel  to  the  axis  of  the  rods  and  polished  for 
examination. 

Tliree  different  structures  resulted  from  these  couples:  no  adherence, 
adherence  with  no  evidence  of  reaction,  and  adherence  with  strong  evidence  of 
reaction  and  the  formation  of  an  unexpected  phase.  The  couples  numbered  3,  4, 
and  5  did  not  adhere  or  react,  a  failure  attributed  to  insufficient  pressure  applied 
by  the  clamp  during  the  anneal . 
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Table  7 

LIST  OF  DIFFUSION  COUPLES  OF  Re  VERSUS  PG 
AND  OTHER  GRAPHITES 


Coiqple 

Number 

Couole  Comoosition 

Temp. 

rc)_ 

Time 

(hr) 

1 

Re  vs. 

PG<a> 

1800 

3.0 

2 

Re  vs. 

PG<a> 

2100 

2.5 

3 

Re  vs. 

PG<a> 

1960 

2.0 

4 

Re  vs. 

PG<a> 

2070 

3.0 

5 

Re  vs. 

PG<a> 

2150 

3.0 

6 

Re  vs. 

PG<a> 

2100 

3.0 

7 

Re  vs. 

PG<a> 

2180 

3.0 

8 

PG<a>  vs.  Re  vs.  PG<a> 

2200 

4.0 

9 

Re  vs. 

ZTA  Graphite 

2200 

4.0 

10 

PG<c>vs.  Revs.  Carbon 
(non-gri4>hitic)  vs.  Re  vs. 
PG<a> 

2300 

6.0 

In  couples  1  and  2,  the  Re  adhered  to  the  PG,  but  in  couple  1  no  evidence  of 
reaction  was  resolved  at  1500X.  In  couple  2,  several  regions  near  the  interface 
exhibited  a  melted  two-phase  structure  with  graphite  flakes  also  present.  This 
structure  is  shown  in  Fig.  3  where  the  black  phases  are  PG  and  graphite  flake, 
the  white  phase  Re,  and  the  gray,  an  unexpected  new  phase.  Because  this  gray 
phase  was  formed  in  the  presence  of  carbon  and  was  attacked  by  chromic  acid, 
it  was  thought  to  be  a  carbide  (105).  The  line  of  small,  white  particles  across 
the  bottom  of  the  photographs  probably  denotes  the  original  interface  since  the 
carbon  inside  this  line  docs  not  have  the  original  optical  anisotropy  of  PG.  Tliis 
microstructure  will  be  discussed  in  more  detail  in  subsequent  paragraphs. 
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Couples  6  and  7  adhered  but  both  couples  separated  at  the  interface  during 
handling,  leaving  some  reaction  product  on  the  Re.  Couple  7  was  mounted  and 
polished  for  microscopic  examination.  This  examination  of  7  showed  the  same 
kind  of  structure  as  in  couple  2.  The  reaction  product  from  6  was  scraped  off, 
and  a  Debye -Sc herrcr  x-ray  diffraction  pattern  was  made  from  the  powder. 

The  powder  ]rjattem,  obtained  with  Ni -filtered  Cu  radiation  in  a  57.3  mm 
camera,  contained  continuous  lines  of  graphite  and  Re  and  an  additional  set  of 
spotty  or  discontinuous  lines,  lliese  additional  lines  indicated  the  presence  of 
a  small  number  of  grains  of  a  different  phase .  Table  8  includes  an  indexed  list 
of  lines  observed  in  this  diffraction  pattern  together  with  all  the  other  lines  of 
this  phase  observed  in  other  diffusion  couples. 

Couples  8,  9,  and  10  also  adhered  and  showed  strong  evidence  of  reaction, 
thus  demonstrating  that  Re  will  react  with  other  forms  of  carbon  below  2300^C . 
Deby./-Scherrer  patterns  of  scrapings  from  couples  8  and  9  were  weak  and  incon* 
elusive,  but  the  pattern  from  10  was  fairly  strong,  and  gave  all  the  lines  listed 
in  Table  8 . 

To  improve  the  diffraction  patterns,  several  attempts  were  made  to  produce 
larger  quantities  of  the  new  phase  by  heating  pressed  pellets  of  Re  and  C  powder 
of  different  compositions.  In  all  of  the  attempts,  no  lines  or  only  the  first  two 
or  three  lines  of  the  new  phase  were  observed.  Although  there  were  good  indi¬ 
cations  that  the  new  phase  was  a  carbide,  the  phase  may  possibly  have  formed 
from  impurities  in  the  argon.  To  be  certain  that  the  new  phase  was  not  an  oxide 
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Table  8 


COMPILATION  OF  DISCONTINUOUS  DEBYE-SCHERRER 
LINES  OBTAINED  FROM  MATERIAL  TAKEN  FROM 
Re-PG  DIFFUSION  COUPLES 


d 

I 

(h^  +  +  1^) 

a 

2.8 

s(*) 

1 

2.80 

1.98 

S 

2 

2.80 

1.63 

S 

3 

2.83 

1.40 

M<^> 

4 

2.80 

1.26 

M 

5 

2.82 

1.15 

M 

6 

2.82 

Miaaing 

7 

0.995 

vw'®> 

8 

2.82 

0.935 

vw 

9 

2.81 

0. 886 

V'V 

10 

2.78 

0. 852 

VW 

11 

2.82 

0.814 

vw 

12 

2.82 

0. 800 

vw 

13 

2.89 

0. 780 

vw 

14 

2.92 

(a)  Strong 

(b)  Medium 

(c)  Very  weak 
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or  nitride,  Re  powder  was  heated  in  air  at  several  temperatures  up  to  900^C 
at  which  all  the  Re  powder  volatilized  by  oxidation.  X-ray  diffraction  and 
chemical  etct.  ig  techniques  produced  no  correlations  between  these  reaction 
products  and  the  phase  produced  in  the  Re-PG  diffusion  couples. 

Fourteen  lines  on  the  Debye -Scherrer  pattern  of  the  new  phase  were 
observed,  the  analysis  of  which  indicated  that  the  new  phase  has  a  simple 
cubic  lattice  with  a  lattice  parameter  of  2.8-f>  A  (106).  The  fact  that  the  atomic 
diameter  of  Re  in  its  own  lattice  is  2.74  A  (107)  suggests  that  the  unit  cell  might 
contain  Re  atoms  at  the  comers  with  carbon  occupying  some  internal  position(s). 
The  cube  center  would  seem  the  most  likely  carbon  position  (CsCl  structure). 

This  position  is  2 .0  A  in  diameter  along  the  cube  diagonal  and,  since  the  carbon 
atom  is  generally  taken  to  be  1.52  A  in  diameter,  there  is  ample  room  for  the 
carbon  atom.  The  stoichiometry  could  not  be  determined  from  the  relative 
intensities  of  the  diffraction  patterns  because  of  their  discontinuous  nature  and 
the  large  difference  in  atomic  scattering  factors .  Chemical  analyses  were  not 
attempted  because  of  the  jmall  amount  present;  however,  electron  microprobe 
analysis  might  be  appropriate  here . 

There  have  been  some  indications  that  Re  will  form  a  carbide  under  certain 
circumstances.  A  carbide  of  Re  has  been  reported  to  form  between  470  and  600°C 
\/ith  "activated  Re"  and  CO  gas,  but  it  decomposed  above  1600°C  (109).  Other 
workers  have  observed  an  unknown  pha5e(s)  in  the  Re-W-C  system,  but  their 
reported  x-ray  diffraction  lines  do  not  correspond  to  any  of  the  strong- lines  of 
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the  new  phase  (110).  Kempter  and  Nadler  (111)  reported  that  carbides  of  Ru  and 
Os,  isomorphous  with  hexagonal  WC,  can  be  produced  under  certain  pressure 
and  temperature  conditions.  With  the  exception  of  Tc,  about  which  little  is 
known,  the  findings  of  Kempter  and  Nadler  leave  Re  completely  surrounded  by 
carbide  formers  in  the  periodic  table,  another  indication  that  Re  might  form  a 
carbide  -  not  necessarily  a  stable  one  -  under  certain  conditions.  Another 
indication  that  a  carbide  of  Re  would  be  metastable  has  been  given  by  Searcy 
and  Finnic  (112)  who  postulated  that  the  heat  of  formation  of  carbides  of  the  Pt 
group  metals  (Ru,  Os,  Rh,  Ir,  Pd,  Pt)  and  Re  should  be  positive. 

All  of  the  above  evidence  strongly  indicates  that  the  new  phase  is  a  carbide 
of  Re.  If  that  is  the  case,  then  the  melted  structure  shown  in  Fig.  3  must  be 
the  results  of  a  eutectic  reaction.  The  melted  structure  cannot  be  the  result 
of  a  peritectic  reaction  because  the  reaction  occurred  below  the  melting  point 
of  pure  Re .  The  graphite  flakes  within  the  carbide  must  have  been  produced  by 
precipitation  from  the  carbide  on  cooling;  precipitation  of  carbon  flakes  often 
occurs  upon  cooling  carbon-saturated  phases  (113).  The  rhenium-rhenium  car¬ 
bide  eutectic  temperature  was  not  accurately  determined  in  this  study,  but  as 
indicated  by  couples  1  and  2,  it  must  lie  between  1800  and  2100°C,  which  is 
below  the  rhenium -carbon  eutectic  temperature  of  2480°C  (103).  This  low 
eutectic  tempo’  ature  is  further  evidence  that  rhenium  carbide  is  metastable 
because,  as  demonstrated  by  the  metastable  Fe-Fc<jC  phase  diagram  relative 
to  the  stable  Fe-C  phase  diagram  (109),  phase  reactions  with  metastable  com¬ 
pounds  occur  at  lower  temperatures  than  do  reactions  with  stable  compounds. 
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To  summarize,  metallic  rhenium  will  react  with  PG  and  graphite, 
near  2100^C  to  produce  a  previously  unreported  phase.  This  new  phase  is 
gray  in  color,  is  attacked  by  chromic  acid,  and  produces  a  Debye -Schcrrer 
pattern  consistent  with  a  simple  cubic  lattice  with  a  unit  cell  edge  of  2.8+  A. 
There  were  many  indications  that  this  new  phase  is  a  metastable  carbide  of 
rhenium  but  the  stoichiometry  of  the  carbide  could  not  be  determined. 

C.  Diffusion  Studies  with  Re 

As  noted  in  Chapter  IV,  Re  was  initially  selected  for  the  diffusion  studies 
because  of  its  low  vapor  pressure  and  fairly  convenient  radioactive  isotope. 
Techniques  for  preparing,  annealing,  and  sectioning  specimens  were  developed 
for  Re  diffusion,  and  the  details  of  these  techniques  are  described  in  Chapter  VI. 
Some  observations  made  during  the  course  of  the  Re  diffusion  study  are  pre¬ 
sented  here. 

PG  specimens  approximately  one  cm  were  prepared  by  grinding  and  lapping 

as  described  in  Chapter  VI.  Following  preparation,  the  edges  of  the  samples 

were  coated  with  a  vinyl  emulsion  and  the  radioactive  Re  electro -deposited 

-9 

from  a  0.1  N  HNO^  acid  solution  containing  10  molar  HRcO^.  Deposition 
occurred  with  the  PG  as  cither  the  anode  or  the  cathode.  When  PG  was  used  as 
the  anode,  it  was  attacked  by  the  nascent  oxygen  and  the  surface  was  covered  with 
a  thin  layer  of  PG  powder.  Deposition  with  PG  as  the  anode  was  probably  just 
HReO^  in  solution  trapped  in  the  cracks  and  around  the  particles  left  by  the  oxygen 
attack.  For  short  deposition  times,  the  damage  caused  by  the  nascent  oxygen 


was  slight,  and  anodic  deposition  was  used  in  most  of  the  experiments  because 
it  appeared  to  be  more  efficient.  The  method  of  deposition  affected  the  resolu¬ 
tion  of  the  interface  slightly  but  did  not  seem  to  affect  the  diffusion  of  Re  into 
the  PG.  Following  deposition,  the  vinyl  emulsion  was  stripped  from  the  PG  and 
the  specimens  were  given  diffusion  anneals  at  high  temperature  (1100  to  1800°C). 

As  described  in  Chapter  VI,  sections  were  ground  from  the  samples.  In  the 
initial  lapping  experiments,  the  ground  surface  of  each  specimen  was  cleaned 
with  pressure -sensitive  mylar  tape  under  fingertip  pressure.  Later,  it  was  dis¬ 
covered  that  the  uneven  pressure  left  a  residue  of  radioactive  material  from  the 
previous  section  on  the  sample  surface.  Because  the  activity  of  the  Re  remain¬ 
ing  on  the  surface  was  very  high,  this  small  residue  significantly  distorted  the 
measured  concentration  profile.  With  this  distortion,  the  logarithm  of  the 
apparent  concentration  decreased  linearly  with  distance  over  three  to  five  decades 

_3 

and  extended  into  the  specimen  approximately  2  x  10  in.  When  the  cleaning 
pressure  was  applied  uniformly  over  the  surface  with  a  flat  platen,  the  amount 
of  residue  was  decreased  considerably,  leading  to  better  resolution  of  the  inter 
face,  as  shown  in  Fig.  4. 

A  fairly  flat  tail  was  observed  beyond  the  initial  false  profile,  but  the  activity 
was  of  the  order  of  10  c/m.  Because  of  this  low  activity,  the  slope  of  the  line 
could  not  be  resolved  within  reasonable  counting  times.  It  was  assumed  that  this 
flat  portion  represented  true  diffusion  of  Re  into  PG,  possibly  distorted  by  edge 
effects;  but  the  activity  was  deemed  too  low  to  warrant  continuing  the  experiments. 
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PENETRATION  DISTANCE  (lo’^  IN. ) 


Fif.  4  Schematic  Concentration  Profiles  for  Re-186  Diffusion  in  PG. 
Curve  (a)  False  Profile  Caused  by  Non-uniform  Pressure 
Applied  During  Cleaning  Operation. 

Curve  (b)  Resolution  of  Interface  With  Uniform  Pressure 


To  increase  the  counting  rates,  either  the  area  of  the  samples  had  to  be 
increased  by  a  factor  of  100  or  the  specific  activity  of  the  metal  had  to  be  in¬ 
creased  by  a  similar  amount.  It  was  not  practical  to  work  with  large  speci¬ 
mens,  and  Re  with  a  high  specific  activity  could  only  be  produced  at  considerable 
expense.  Therefore,  another  readily  available  metal  with  much  higher  specific 
activity  was  sought.  As  discussed  in  Chapter  IV,  carrier-free  Nb-95  satisfied 
these  requirements.  Furthermore,  the  Nb  produced  counting  rates  of  1000  to 
5000  c/m  for  the  same  conditions  under  which  Re  produced  10  c/m.  As  a  result, 
Nb-95  was  used  to  study  metal  diffusion  in  PG. 
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Chapter  VI 

EXPERIMENTAL  PRCXJEDURES 

A.  Structural  Characterizations 
1 .  Microstructure 

Separate  specimens  from  each  PG,  as-deposited  (AD)  and  heat-treated  (HT), 
were  mounted  for  observation  on  the  plane  of  the  deposit  and  perpendicular  to  it. 
Specimens  polished  perpendicular  to  the  plane  of  the  deposit  first  had  their  depo¬ 
sition  plane  surface  prepared  as  if  they  were  to  be  used  in  a  diffusion  couple. 
Thus,  any  subsurface  effects  of  diffusion  sample  preparation  could  be  observed. 
The  specimens  were  mounted  in  an  epoxy-resin  under  vacuum.  They  were  then 
ground  through  600  grit  paper,  given  alp  diamond  polish,  and  then  given  a 
final  polish  with  Cer-Cre*  and  chromic  acid.  Details  of  the  metallographic 
tecliniquc  have  been  published  elsewhere  (114).  Each  specimen  was  examined 
for  diffusion  sample  preparation  effects,  as  noted  above;  photomicrographs 
were  then  taken  at  10,  100,  and  lOOOX  under  polarized  light,  at  approximately 
15  deg  from  extinction.  The  structure  produced  by  polishing  the  plane  parallel 
to  the  deposition  surface  did  not  clearly  show  the  nodular  structure  of  PG,  and 
tins  teclinique  was  abandoned  in  favor  of  cleaving  the  specimens  parallel  to  the 
deposition  surface  with  a  razor  blade. 

'“‘llurell  Corporation,  2223  Fifth  Avenue,  Pittsburgh,  Pa. 

64 


I 


2 .  Transmission  Electron  Microscopy 

Sjjecimens  for  transmission  electron  microscopy  were  obtained  by  touching 
a  collodion -covered  grid  to  the  PG  powder  produced  when  a  sample  was  ground 
in  the  Goldstein  apparatus*  In  order  to  have  large  flakes  for  observation,  only 
samples  which  had  been  ground  on  the  plane  parallel  to  the  deposition  surface 
were  used.  The  specimens  were  examined  in  transmission  using  an  Hitachi 
HU -11  electron  microscope  operating  at  100  kV. 

3.  X-ray  Diffraction  Studies 

In  this  rection,  the  procedures  for  determining  the  stacking  error  proba¬ 
bility,  unit  cell  parameters,  crystallite  size,  degree  of  preferred  orientation, 
and  orientation  dependence  of  structure  will  be  given . 

As  discussed  in  Chapter  III,  the  stacking  error  probability,  p  ,  in  graphitic 
carbons  can  be  experimentally  related  to  the  observed  interlayer  spacing  (Fig.  2) 
and,  to  soipe  extent,  related  by  theory.  Because  PG  is  also  a  graphitic  carbon, 
this  observed  dependence  will  be  used  to  determine  p  .  The  use  of  Fig.  2 
requires  the  determination  of  a  precise  value  of  dQQQ2'  Accurate  dQQQ2  values 
were  determined  by  extrapolating  the  d-spacing  obtained  from  the  0002  ,  0004, 
0006,  and,  where  possible,  0008  diffraction  peaks  to  29  =  180  deg  using  the 
Taylor-Sinclair  function  to  provide  a  linear  extrapolation  (115).  The  dif¬ 
fraction  profiles  were  obtained  with  a  General  Electric  XRD-5  diffraction  unit 
with  a  Ni-filtered  Cu  radiation,  line  source.  The  peaks  were  scanned  automatic¬ 
ally  with  a  low  time  constant  (0.5  seconds)  and  slow  scanning  speed  (0.1  deg/min) 

*See  Section  B-1  of  this  Chapter. 
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for  the  0002  aiid  0004  profiles,  and  fast  speeds  (1  deg/min)  for  the  0006  and 
0008  profiles.  The  slit  system  utilized  the  0.4  deg  beam  slit,  high-resolution 
Soller  slit,  and  0.1  deg  receiving  slit.  The  specimens  were  1  mm-thick,  flat 
slabs  with  faces  parallel  to  the  deposition  plane,  and  they  were  placed  in  the 
usual  symmetric  Bragg  reflection  geometry .  After  a  precise  value  of  dQQQ2 
had  been  obtained,  the  solid  line  in  Fig.  2  was  used  to  determine  p  . 

The  a  parameter  was  determined  from  the  10  and  bands  in  PG  by  point 
counting  at  0.1  deg  increments  in  the  vicinity  of  the  peak  maximum.  As  noted  in 
Chapter  III,  the  d-spacing  correction  for  the  hk  bands  is  negligible  for  PG  with 
a  crystallite  size  greater  than  100  A. 

There  are  two  crystallite  dimensions  to  be  measured  in  graphites:  the  thick¬ 
ness,  L  ,  and  a  characteristic  diameter,  L  .  The  thickness,  L  ,  was  deter- 
c  a  c 

mined  from  the  0002  ,  0004,  and  0006  profiles.  These  profiles  were  smoothed 
by  hand,  and  the  integrated  intensity,  I(of),  integral  breadth,  B,  and  center  of 
gravity,  20,  were  calculated  using  the  trapezoidal  rule  and  appropriate  equations 


1(0)  =  2  P(20)A20 

B  =  1(0  )/I 

'  '  max 


20  = 


2  20  P(20)  A20 
1(0) 


Here  P(20)  is  the  height  of  the  profile  at  20  and  fhe  peak  height.  The 
integral  breadth  was  corrected  by  use  of  the  Warren  (87)  correction 
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which,  when  applied  to  the  0002  profile,  was  found  to  give  results  similar  to 
those  of  the  more  rigorous  Fourier  analysis  (see  Chapter  III).  Crystallite 
thicknesses  were  then  computed  with  the  Scherrer  equation,  using  the  corrected 
integral  breadth  p 

L  »  - - — 

c  p  cos  6 

The  value  of  b  was  taken  from  a  strain -annealed  PG  specimen  with  the  same 
orientation  and  similar  in  size  to  the  regular  specimens.  When  standards  of 
the  same  absorption  coefficient  and  geometry  are  used  to  correct  integral 
breadths,  separate  absorption  corrections  arc  not  needed.  If  the  absorption 
effects  are  so  large  that  the  ^^^2  not  noticeably  distort  the  stan¬ 

dard  profile,  then  this  correction  need  not  be  made  either.  These  simplifica¬ 
tions  were  applicable  for  the  PG  line  breadths  in  this  part  of  the  study.  With 
the  use  of  equations  developed  b)  Bragg  and  Packer  (93),  the  0002  standard  pro¬ 
file  was  corrected  for  absorption  and  found  t«)  be  only  10%  greater  than  the  dif¬ 
fraction  profile  of  a  polycrystalline  quartz  standard.  This  result  was  taken  as 
sufficient  justification  to  use  the  strain -annealed  PG  as  the  standard.  The  absorp¬ 
tion  corrections  were  programmed  by  Scott  and  Sackman,  and  calculated  on  a 
General  Electric  Computer  System  (116).  . 

The  and  Ij^  peaks  in  PG,  obtained  from  flat  specimens  cut  perpendicular 
to  the  deposition  plane,  were  used  to  determine  the  crystallite  diameter,  L^. 

With  the  specimens  in  the  symmetric  Bragg  reflecfon  geometry,  the  profiles 
were  determined  by  point  counting  at  0.2  deg  increments.  The  various 
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corrections  for  these  profiles  were  applied  in  the  following  order: 

1.  Packgrounds,  determined  at  the  low  IB  side  of  the  1^  profile 

2.  Absorption*,  corrected  with  the  relation  derived  by  Bragg  and 
Packer  (3)  and  programmed  by  Sackman  and  Scott  (116) 

3.  ^1^2  corrected  with  the  Rachinger  equation  (91) 

4.  Preferred  orientation,  corrected  with  the  equation  derived  by  Guentert 
and  Cvickevich  (92) 

5.  Instrumental  broadening,  corrected  with  the  Warren  correction  (90), 

2  2  2 

^  *  B  -  b,  with  the  220  line  of  NaCl  and  331  line  of  Si  as  standards. 
The  doublet  was  also  corrected  in  the  standards. 

The  correction  for  preferred  orientation  requires  that  the  pole  figure  data 
(described  later)  be  represented  as  I  *  Iq  cor.”'  a  .  The  value  of  m  was 
determined  from  a  log-log  plot  of  I  vs  cos  a  . 

Small,  rectangular  specimens,  0.5mm  x  0.5mm  x  20mm,  were  used  for 
studying  preferred  orientation.  The  specimens  were  prepared  by  mounting 
larger  pieces  on  a  slotted  metal  bar  with  hot  wax  and  abrading  them  against  fine 
emery  paper  backed  with  a  flat  platen .  The  specimens  were  cut  so  that  two  of 
the  faces  were  parallel  to  the  deposition  surface.  These  small  rectangular 
specimens  were  placed  in  a  two-circle  goniometer  and,  utilizing  the  rectangular 
symmetry,  were  optically  aligned  in  the  G.E.  XRD-5 X-ray  Diffractometer. 
Further  refinement  in  the  alignment  was  made  until  the  0002  profiles  from 

*The  mass  absorption  coefficient  used  in  these  corrections  was  4.15,  as 
reported  by  Chipman  (117). 
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opposite  sides  of  the  samples  peaked  at  the  same  20,  within  0.01  deg,  and  the 
intensities  were  the  same  within  20%  of  each  other.  With  the  samples  properly 
aligned,  the  0002  line  profiles  were  scanned  at  increments  of  5  deg  in  a  ,  the 
angle  between  the  normals  of  the  deposition  plane  and  those  crystallites  causing 
diffraction.  At  each  increment  in  ot  ,  the  integrated  intensity,  integral  breadth, 
and  center  of  gravity  were  computed  using  the  trapezoidal  rule  as  above.  With 
these  data,  linear  plots  of  log  I(o)  vs  a  were  made  with  the  intensities  nor¬ 
malized  to  1(0)  =  100  and  the  mean  crystallite  lilt  angles, defined  by 


<a 


1(0)  = 


were  calculated  from  the  linear  plots.  Past  experience  (5)  and  data  taken  here 
2  1/2 

showed  that  <a  >  is  insensitive  to  the  specimen  size  over  the  range  0.5 
to  1.0  mm,  indicating  that  with  the  above  arrangement  there  is  no  effect  of 
specimen  geometry  on  the  pole  figure  below  1mm  specimen  thickness. 

The  orientation  dependence  of  structure  as  noted  by  plotting  the  values  of 
p  derived  from  the  d-spacings  and  the  observed  integral  breadths  vs  a  .  The 
integral  breadths  were  not  corrected  for  the  various  effects,  but  any  significant 
changes  with  a  would  show  in  the  uncorrected  values. 


69 


4.  Small  Angle  Scattering 

Small  angle  scattering  of  x-rays  by  PC  was  measured  with  a  special  fixture 
designed  by  Bragg  (118)  for  the  Norelco  X-ray  Diffraction  machine.  The  slit  sys¬ 
tem  limits  the  incident  beam  to  1/12  deg  divergence  and  also  limits  the  aperture 
at  the  receiving  slit  to  1/6  deg.  The  fixture  permits  very  precise  and  repro¬ 
ducible  adjustment  of  the  slit  system  so  data  can  be  taken  very  close  to  20  =0. 
Specimens  5mm  x  30mm  x  1  mm  tiiick  were  prepared  parallel  and  perpendicular 
to  the  plane  of  deposit  and  were  tlie  same  specimens  as  tliose  on  which  the  unit 
cell  dimensions,  crystallite  sizes,  and  densities  were  measured.  Data  were 
taken  by  point  counting  with  tlie  Norelco  Geiger  Counter  and  were  then  corrected 
for  counting  losses  and  background.  The  background  was  determined  by  placing 
two  thicknesses  of  pure,  annealed  A1  foil  (0.075mm  total  thickness)  in  the  speci¬ 
men  position  and  point  counting  in  the  same  range  in  which  the  small  angle 
scattering  data  were  taken.  The  transmitiances  of  the  A1  and  PC  samples  were 
determined  by  measuring  the  attenuation  of  the  diffracted  lx;am  from  the  1 1 1 
planes  of  quartz,  using  Ni -filtered  Cu  radiation.  All  the  transmittances  were 
0.40  ±  0.02  so  relative  transmittance  corrections  did  not  have  to  be  made.  To 
be  certain  tliat  the  A1  foil  did  not  protiuce  any  scattering  of  its  own,  the  back¬ 
ground  was  observed  with  and  without  the  A1  foil  and  for  angles  greater  than 
0.12  deg  20  ;  the  transmittance  of  the  A1  was  the  same  as  that  determined 
with  the  monochromatic  radiation . 


70 


1 


5.  Other  Structural  Parameters 


The  structu'^e  of  each  PG  was  further  characterized  using Hgporosimetry  and 
bulk  density  measurements. 

Small  specimens  (5mm  x  5mm  x  1  mm)  of  each  PG  in  the  as-deposited  and 
heat-treated  conditions  were  separately  placed  in  the  penetrometer  of  an  Aminco 
Hg  porosimeter  and  the  volume  displaced  in  Hg  at  pressures  between  1 .8  and 
5000  psi  pressure  was  measured. 

The  bulk  densities  of  three  specimens  from  each  PG,  us-deposited  and  heat- 
treated,  were  determined  by  drilling  a  hole  in  the  end  of  each  specimen,  sus¬ 
pending  it  from  a  small  wire,  and  weighing  it  in  and  out  of  isopropyl  alcohol. 

The  weighing  was  done  to  the  nearest  0.01  mg  with  a  total  sample  weight 
larger  than  0.20  gm. 

B.  Diffusion  Measurements 

1.  Method  of  Measuring  Diffusion  Coefficients 

A  summary  of  methods  of  measuring  the  diffusion  coefficient,  D  ,  in  solids 
has  recently  been  given  by  Cadek  (119).  Most  of  the  teclmiques  he  described 
involve  determining  the  composition  profile  by  any  of  several  methods,  but 
various  techniques  have  also  been  developed  in  which  D  is  computed  from  the 
decrease  in  surface  activity  or  from  the  total  amount  of  material  that  has  dif¬ 
fused  into  or  out  of  a  phase.  All  the  latter  methods  assume  the  nature  of  the 
composition  profile,  and  the  assumptions  are  exact  only  for  self-diffusion  in 
a  single  crystal.  To  be  certain  of  the  interpretation  of  the  results,  the  compo¬ 
sition  profiles  for  metal  diffusion  in  PG  had  to  be  determined. 
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Profile  determinations  can  be  divided  into  two  general  categories:  those 
in  which  the  sample  is  divided  into  several  sections  perpendicular  to  the  mass 
flux  with  subsequent  analysis  of  each  section,  and  those  in  which  the  specimen 
is  sectioned  only  once  parallel  to  the  flux  with  the  composition  gradient  deter¬ 
mined  by  measuring  some  property  related  to  the  composition.  Properties 
commonly  used  for  this  determination  include  hardness,  color,  autoradiographic 
intensity,  and  fluorescent  x-ray  intensity  excited  by  x-rays  or  an  electron  beam. 
Because  the  very  small  concentration  of  diffusing  species  in  PG  determined  the 
second  of  these  categories,  a  method  of  sectioning  perpendicular  to  the  mass 
flux  was  sought,  to  be  combined  with  tracer  analytical  techniques.  A  propoga- 
tion  of  error  analysis  (120, 121)  showed  that  the  best  techniques  were  those  in¬ 
volving  related  weight  measurements  (determining  the  weight  loss  of  the  sample 
after  each  sectioning)  and  related  thickness  measurements  (determining  the 
sample  thickness  after  each  sectioning).  The  related  thickness  measurement 
technique  was  chosen  because  the  density  of  PG  is  low  and  potentially  variable, 
and  PG  produces  fine  dust  upon  grinding  which  could  cause  difficulties  in  the 
weight  measurements  by  adhering  to  the  sample.  If  sufficient  weight  were 
removed  with  each  section,  this  dust  would  not  be  a  problem;  however,  larger 
samples  and/or  greater  penetration  would  be  required.  In  any  case,  a  thick¬ 
ness  measurement  is  faster  and  more  convenient  than  a  precision  weight  meas¬ 
urement.  The  method  used  employed  a  precision  lapping  device  built  by  Gomez 
(122)  and  similar  to  one  described  by  Goldstein  (123).  The  amount  removed 
with  each  section  and  the  penetration  distances  were  measured  with  a  high 
sensitivity  dial  gage. 
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2.  Sample  Preparation  Prior  to  Diffusion 

Nearly  flat  sheets  of  PG,  10cm  x  10cm  x  1cm  thick  were  cut  into  0.7cm 
square  specimens.  The  original  sheets  had  a  slight  cylindrical  curvature,  and 
the  first  long  cuts  were  made  parallel  to  the  axis  of  the  cylindrical  curve  to 
minimize  its  influence.  Some  specimens  from  each  of  the  PC's  were  then 
given  a  graphitization  heat  treatment  at  2700®Cfor  15 minutes.  These  speci¬ 
mens  were  labeled  IIT  (Heat-Treated)  and  the  untreated  samples,  AD  (As- 
Deposited). 

The  Goldstein  apparatus  (122, 123)  consisted  of  a  rotating  shaft  perpendicu¬ 
lar  to  a  small  table  upon  which  a  circular  glass  flat  was  eccentrically  rotated. 
A  cylindrical  precision -ground  hardened -steel  holder  was  made  to  just  fit  onto 
the  shaft,  the  holder's  top  and  bottom  faces  being  precisely  perpendicular  to 
the  shaft  axis  and  parallel  to  the  glass  flat  surface.  The  holder  was  1/2  inch 
in  diameter,  the  shaft  fitted  into  the  holder  along  3/4  in.  of  its  length,  and 
the  faces  of  the  holder  were  parallel  within  1  x  10  ^  in.  and  perjiendicular  to 
the  shaft  axis  within  ±2  x  10’^  radians. 

For  diffusion  perpendicular  to  the  deposition  plane  (to  be  labeled  the  <c> 
direction)  the  specimen  faces  perpendicular  to  the  intended  flux  were  ground 
flat  and  parallel.  First,  the  substrate  surface  was  attached  to  the  holder  with 
hot  wax  and  the  final  deposition  surface  ground  flat  with  adhesive  240  grit  paper 
which  was  smoothly  affixed  to  the  circular  glass  flat.  The  specimen  was  then 
removed,  remounted  for  grinding  the  substrate  surface,  and  then  remounted 
again  for  final  preparation  of  the  surface  nearest  the  final  deposition  surface. 
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In  this  fashion,  the  resulting  diffusion  flux  was  perpendicular  to  the  substrate 
surface  and,  thereby,  parallel  to  the  c-axis  of  most  of  the  crystallites  (sec 
preferred  orientation  results). 

For  diffusion  parallel  to  the  deposition  plane  (to  be  noted  by  <a>)  the  speci¬ 
mens  were  first  mounted  on  one  edge,  ground,  then  mounted  on  the  other  edge 
and  given  the  final  preparation.  Care  was  taken  to  ensure  that  the  finished  sur¬ 
face  was  perpendicular  to  the  substrate  surface.  Because  the  initial  nucleation 
layer  did  not  have  the  same  structure  as  the  rest  of  the  PC  it  was  removed  after 
the  final  surface  preparation  by  grinding. 

Considerable  care  was  taken  to  ensure  that  the  flat  specimens  were  attached 
parallel  to  the  holder  face  and  glass  flat.  The  holder  was  heated  to  about  200^C 
on  a  hot  plate,  and  a  small  amount  of  Apiezon  wax  melted  on  its  flat  surface. 
With  the  holder  held  with  one  pair  of  tweezers,  a  hot  specimen  was  put  down 
onto  the  melted  wax  and  moved  back  and  forth  until  intimate,  uniform  contact 
with  the  holder  was  established.  The  specimen  tweezers  were  then  carefully 
removed  and  the  holder  set  on  a  copper  plate  to  cool.  With  this  technique, 
specimens  could  be  removed  and  replaced  on  the  holder  with  less  than  7/i  varia¬ 
tion  over  the  entire  surface. 

Final  preparation  consisted  of  grinding  through  4/0  paper  and  giving  the  speci  - 
men  a  metallographic  polish  with  Cer-Cre  and  chromic  acid  on  slipper  satin. 

The  adhesive  paper  was  smoothly  attached  to  the  glass  flat.  The  abrasive  paper 
was  affixed  to  the  flat  and  a  piece  of  ordinary  paper  placed  over  the  abrasive. 
Then,  the  adhesive  paper  was  smoothed  by  pressing  it  against  the  straight  edge 
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of  u  glass  sheet  and  dragging  it  across  that  edge.  The  ordinary  paper  on  top 
of  the  abrasive  {prevented  raking  off  the  abrasive.  Non-adhesive  3/0  and  4/0 
papers  were  affixed  with  cellophane  tape  at  two  edges.  The  glass  flat  for 
metallographic  polishing  was  prepared  by  affixing  a  layered  structure  of  slipper 
satin,  parafilni,  and  double-back  masking  tape  to  it  under  heat  and  pressure. 

The  grinding  action  took  place  with  only  the  weight  of  the  steel  shaft  (1  cm 
diameter  by  15cm  long)  on  the  specimen  but  faster  metallographic  polishing 
occurred  with  an  additional  weight  of  alx)ut  300  grams. 

The  tracer  was  deposited  on  the  metallographically  polished  surface  by 
evaporating  a  drop  of  suitably  diluted  Nb-oxalate  solution  procured  from  the 
Oak  Ridge  Radioisotope  Facility  (104).  Controlled  amounts  were  deposited 
with  a  0.050  cm  pipette,  and  about  50  x  10  c/m  were  deposited  on  each 
specimen . 

3.  Furnace  Description  and  Annealing  Procedure 

The  furnace  used  for  the  heat-treatment  and  diffusion  anneals  consistetl 
of  a  vertical  graphite  resistance  clement,  with  power  supplied  at  water-cooled 
clamps  located  at  one  end  of  the  element.  The  cylindrical  element  was  uniformly 
slotted  in  several  places  to  provide  a  longer  resistance  path,  and  two  sighting 
slits  approximately  1/4  in.  x  2  in.  were  located  in  the  element.  The  furnace 
was  originally  built  for  tensile  testing  at  3000^C,  and  the  sighting  slits  were 
designed  for  observing  the  specimen  gage  length.  The  furnace  was  equipped 
with  a  50  kW  power  supply,  mechanical  and  diffusion  vacuum  pumps.  He  gas 
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supply,  provision  for  protecting  the  quartz  sight  glasses  from  contamination, 
and  a  quenching  chamber  witli  an  air-lock.  Considerable  insulation  was  supplied 
by  several  concentric  shells  of  graphite  with  the  volume  in  Ixjtween  filled  with 
graphite  felt. 

The  furnace  power  was  controlled  by  a  saturable  reactor,  and  an  intcgral- 
plas-proportional  controller,  designed  by  VV .E .  Kppler  of  LMSC,  which  pro¬ 
vided  control.  A  Pyro-Eye  two-color  pyrometer  was  used  to  monitor 

the  temperature.  The  field  of  view  of  the  two-color  pyrometer  was  too  large 
to  inckde  only  the  specimen  holder  surface;  as  a  result,  some  drift  of  holder 
temperature  resulted  from  monitoring  the  combined  heater  plus  sample  tem¬ 
perature.  Tliis  drift  was  corrected  by  hand  and  was  less  than  30^  total 
variation  in  all  but  one  case,  when  a  50°  variation  occurred. 

For  the  diffusion  anneals,  the  specimens  were  placed  face  to  face  in  a 
cylindrical  graphite  holder  wliich  was  attached,  by  pinn  d  connections,  to  a 
0.5  cm  iliameter  x  I5cm  tong  graphite  rod.  This  rod,  in  turn,  was  attached 
to  a  water-cooled  copper  plunger  that  extended  out  of  the  furnace  through  a 
sliding  vacuum  seat. 

I’o  anneal  specimens,  the  holder  was  placed  in  the  hot  zone,  the  furnace 
evacuated,  and  the  temperature  then  increased  manually  to  near  the  desired 
temperature.  At  120vj''C,  the  vacuum  valve  was  closed  and  the  system  brought 
to  a  gage  pressure  of  -10  inches  of  Hg  with  pure  He  gas.  The  heatup  under 
vacuum  served  to  reduce  the  Nb-oxalate  and  to  purify  the  furnace  atmosphere 
while  the  He  gas  served  to  prevent  arcing  at  the  liigher  temperatures,  to  provide 


better  heat  transfer  to  the  holder,  and  to  reduce  the  tendency  of  the  Nb-95  to 
vaporize  during  the  anneal.  Virtually  all  the  Nb-95  remained  on  the  samples; 
no  appreciable  activity  could  be  detected  on  the  holder  or  inside  the  furnace 
after  an  anneal.  The  high  sensitivity  of  the  controller  required  that,  initially, 
the  temperature  be  brought  up  manually.  When  the  temperature  came  within 
100°C  of  the  desired  temperature,  however,  the  controller  was  put  in  the  auto¬ 
matic  mode  with  a  manual  override  and  the  desired  temperature  reached  auto¬ 
matically.  'Die  temperature  was  measured  with  a  Micro -optical  brightness 
pyrometer  focused  on  a  black-body  hole  located  in  the  base  of  the  holder,  just 
below  the  specimens.  Heatup  time  before  the  controller  was  introduced  was 
about  15  minutes,  with  another  10  minutes  being  spent  reaching  the  final  tem¬ 
perature.  No  temperature  overshoot  was  experienced  with  this  heating  method. 
Heatup  time  was  a  small  fraction  of  the  annealing  times . 

During  the  anneal,  the  temperature  of  the  holder  was  monitored  with  the 
manual  pyrometer  and  adjusted  for  the  drift  mentioned  above .  At  the  end  of  the 
anneal,  the  holder  was  raised  out  of  the  hot  zone  and  into  a  water-cooled  copper 
quench  chamber  where  the  holder  cooled  to  below  700®C  in  alxjut  30  seconds. 
When  the  holder  had  cooled  to  about  50°C  (approximately  20  minutes),  the 
couple  was  removed  from  the  quench  chamber  through  the  air-lock.  The  time 
at  temperature  was  measured  to  the  nearest  minute,  with  the  initial  time  esti¬ 
mated  as  the  time  at  which  the  temperature  reached  20°C  below  the  average 
temperature  of  the  anneal . 
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The  manual  Micro -Optical  pyrometer  was  calibrated  with  the  sight  glass 
against  an  NBS  pyrometer  lamp  between  700  and  2200°C  and  then  further  cali¬ 
bration  was  accomplished  by  measuring  the  melting  point  of  Pd  and  Mo  pellets  in 
contact  with  graphite  in  the  furnace.  Nadler  and  Kempter's  cata  (124)  on 
melting  points  of  metals  in  contact  with  carbon  was  used  as  a  reference.  The 
temperature  measurement  errors  are  estimated  as  follows: 

Error  Temperature 

1 .  Due  to  lamp  calibration  ±7°C 

2.  Due  to  temperature  drift  during  run  -^10-30^0 

3.  Estimated  total:  ±15  to  30°C 

between  700  and  2200°C 

The  figures  given  above  do  not  refer  to  the  precision  of  a  temperature  measure¬ 
ment,  which  is  better  than  ±5®C,  but  to  the  accuracy. 

4.  Sectioning  Techniques 

Prior  to  sectioning,  approximately  1-1/2  mm  was  cut  off  the  edges  of  the 
specimens  with  a  jeweler's  saw  and  a  small  steel  miter  box  [WMB  (17)  and  this 
study  showed  that  edge  effects  can  be  appreciable  in  diffusion  studies  in  PGj . 
Following  edge  removal,  each  specimen  was  mounted  on  the  Goldstein  device 
holder  as  described  above,  preparatory  to  lapping. 

Prior  to  each  Upping  operation,  the  combined  thickness  of  sample  and 
holder  was  measured  at  three  to  six  places  on  the  samjile  surface  with  a  dial 
gage  accurate  to  ±1  x  10"^  in.  (Starrett  25-106).  The  specimen  and  holder 
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were  moved  about  on  the  dial  gage  table  to  eliminate  small  dirt  particles  and 
measurements  were  taken  at  the  same  points  on  the  specimen.  With  this  pro¬ 
cedure  the  amount  removed  at  a  particular  point  could  be  measured  to  ±1  x 
10’^  in.  The  specimens  could  be  removed  and  replaced  under  the  dial  gage 
with  the  same  thickness  measurement  obtained.  Thermal  expansion  effects 
were  estimated  to  be  less  than  10'^  in.,  because  the  room  was  air-conditioned 
and  the  sample  was  handled  very  little.  These  factors,  combined  with  a  fairly 
high  thermal  mass  and  low  thermal  conductivity,  made  thermal  expansion  effects 
negligible.  The  surface  of  the  specimen  was  not  perfectly  flat  but  had  a  slight 
cone  shape  due  to  misalignment.  This  cone  shape  was  less  than  10  x  10"^  in. 
high  at  the  center  with  a  base  diameter  the  size  of  the  specimen  (~ 0.375  in.). 

The  first  lap  was  done  with  dry  Ccr-Cre  on  slipper  satin  to  remove  a  maxi¬ 
mum  of  undiffused  material,  if  present,  without  scratching  the  surface  and 
embedding  the  tracer  in  the  scratches.  The  next  lap  was  made  with  4/0  paper 
and  the  next  two  with  3/0  paper.  At  this  point,  approximately  10  —  40  x  10”^  in. 
total  thickness  was  removed  and  this  thickness  is  about  equal  to  the  error  in 
replacing  the  specimen  on  the  holder  after  the  diffusion  anneal .  Following  the 
preliminary  lapping,  the  abrasive  papers  were  selected  according  to  the  amount 
to  be  removed.  The  profiles  were  determined  to  depths  of  3  —  30  x  10  in.  and 
were  analyzed  in  10  to  50  increments.  Slight  variations  in  surface  flatness  pro¬ 
duced  an  average  uncertainty  in  thickness  removed  of  ±2  x  10  ^  in.,  corres¬ 
ponding  to  a  relative  error  of  10%  or  less. 
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After  each  grinding  operation,  the  specimen  was  removed  and  the  surface 
cleaned  with  Scotch-brand  No.  56  yellow  mylar  pressure  tape  by  pressing  the 
sample  face  against  the  adhesive  several  times,  using  a  5  x  8  file  card  on  a  flat 
surface  as  the  backing.  (The  use  of  a  flat  backing  is  very  important  as  uneven 
pressure  allows  material  to  remain  after  each  lap  and.  when  most  of  the  tracer 
has  not  diffused,  this  seriously  distorts  the  composition  profile.)  The  same 
piece  of  tape  was  then  used  to  cover  the  small  square  of  abrasive  paper  which 
contained  the  PG  powder.  Mylar  tape  has  a  peculiar  advantage  in  that  it  acquires 
a  static  charge  which  attracts  the  PG  powder.  This  property  can  be  used  to  con¬ 
trol  the  fine  powder  produced  by  grinding. 

All  the  equipment  used  for  radioactive  counting  was  manufactured  by  Atomic 
Instrument  Co.,  Baird  Associates.  The  abrasive  squares  containing  PG  powder 
and  Nb-95were  placed  inside  a  hollow  Nal  scintillation  counter  (A.I.C.  Model 
810).  The  counting  rate  was  measured  with  a  (3 ~y  spectrometer  composed  of 
A.I.C.  Model  312  power  supply,  215  amplifier,  510  pulse  height  analyzer, 

.134  high-speed  scalar,  and  630  precision  timer.  The  pulse  height  analyzer 
was  set  to  integrate  over  the  0,75  Mev  y  peak  of  Nb  between  0.70  and  0.84 
Mev.  With  this  arrangement,  the  counting  efficiency  was  approximately  25%, 
based  on  the  initial  assay  supplied  by  Oak  Ridge.  Counting  rates  for  each  sec¬ 
tion  were  less  than  10,000  c/m;  generally,  at  least  1000  total  counts  v/ere  taken 
on  each  lap  to  produce  a  standard  deviation  of  approximately  5%  for  a  counting 
rate  as  low  as  300  c/m  and  background  of  70  c/m  (125).  Sections  were  taken 
until  the  residual  activity  fell  Ik*1ow  200  c/m. 
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The  conceiiiruiion  of  Nb-95  was  computed  by  dividing  the  intensity  by  the 
average  thickness  removed;  the  distance  from  the  original  interface  was  com¬ 
puted  from  the  mid-points  of  eacli  lap. 

5.  Auloraihography 

Autoradiographs  were  taken  before  and  during  lapping  by  placing  the  radio¬ 
active  sample  against  Kodak  projector  slide  plates  and  exposing  for  2  to  800 
hours,  depeiuling  upon  the  intensity.  Generally,  the  as-ground  surfaces  were 
used;  on  occassion,  1  iwever,  the  surface  was  repolished  before  an  autoradio¬ 
graph  was  taken.  Repolishing  did  not  alter  the  appearance  of  the  autoradiograph 
other  Ilian  to  improve  its  resolution.  Copies  of  the  autoradiographs  were  made 
by  photographing  the  negatives  at  lOX  using  transmitted  light  and  a  bellows 
camera . 

In  this  fashion,  the  relative  film  darkening  is  preserved,  and  the  darker 
regions  represent  areas  of  higher  Nb  concentration. 
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Chapter  VII 


EXPERIMENTAL  RESULTS  -  STRUCTURE 

A.  Microstructurc 

The  microstructures  of  the  as -deposited  PC's  used  in  this  study  are 
shown  in  Fig.  5  and  the  heat-treated  structures  are  shown  in  Fig.  6.  The 
magnification  used  in  these  photomicrographs  is  lOX.  As  may  be  seen,  these 
PC's  represent  a  variety  of  microstructures.  The  R277  material  is  charac¬ 
terized  by  a  coarse  substrate  nucleated  structure  (structure  A2  in  Chapter  III) 
with  a  moderately  strong  secondary  structure.  The  R165  material  is  also 
substrate  nucleated  but  has  a  fine  primary  structure  and  a  strong  secondary 
structure  (structure  C3).  The  R1  PG  is  a  good  example  of  regeneratively 
nucleated  PG  (FI).  The  as -deposited  R277  PG  shows  a  tendency  to  delaminate, 
as  indicated  by  the  photomicrographs  perpendicular  to  the  deposition  plane. 

As  described  in  Chapter  VI,  the  samples  polished  perpendicular  to  the 
deposition  plane  first  had  their  deposition  plane  .surface  prepared  as  if  to  be 
used  in  a  diffusion  couple.  The  metallographic  preparation  of  these  samples 
was  such  that  the  edges  could  be  examined  within  one  micron  of  the  surface 
prepared  for  diffusion.  No  distortion  or  damage  due  to  the  diffusion  sample 
preparation  was  found . 

The  heat-treatment  did  not  noticeably  alter  the  structure  as  observed  at 
lOX  magnification  (Fig.  6),  but  it  did  produce  some  sort  of  speckling  over  the 
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:  R277  PG.  (cl),  (c2);  Rl  PG 
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(al),  (a2):  R165  PG  Parallel  and  Perpendicular  to  Depositon  Plane, 
(bl).  (b2):  R277  PG.  (cl),  (c2):  Rl  PG 


surface  perj^cndicular  to  the  deposition  plane.  The  source  of  this  speckling 
is  not  know.T,  but  it  does  not  appear  in  the  photomicrographs  at  higher  magni¬ 
fication.  The  speckling  is  probably  due  to  a  shadow  effect  produced  by  slight 
depressions  and  the  oblique  lighting  used  to  take  these  photographs.  The 
tendency  of  R277  material  to  delaminate  is  further  shown  in  the  heat-treated 
samples,  a  similar  effect  being  noted  in  the  R1  PG.  The  R1  material  actually 
consisted  of  a  variety  of  structures,  as  shown  in  Fig.  6  —  C2,  where  both 
substrate  and  regeneratively  nucleated  structures  are  in  evidence. 

Figures  7  through  12  are  photomicrographs  of  each  PG  at  100  and  lOOOX 
magnification.  The  photomicrographs  arc  arranged  so  that  the  as-deposited 
and  heat-treated  specimens  of  each  PG  in  each  orientation  are  together,  and 
consecutive  figures  contain  the  views  parallel  and  perpendicular  to  the  deposi¬ 
tion  surface.  The  views  parallel  to  the  deposition  plane  are  of  cleavage  surfaces, 
and  cleavage  defects  cause  some  regions  of  the  photographs  to  be  out  of  focus. 

The  main  features  to  be  noted  in  this  series  of  photographs  are  the  similarities 
of  the  R165  and  R277  structures,  as  contrasted  with  the  R1  PG,  the  lack  of 
change  in  primary  structure  with  heat-treatment,  and  the  partial  elimination  of 
the  finer  secondary  structure  in  the  heat-treated  material.  The  specimens 
were  prepared  and  photographs  by  Mr.  A.S.  Gleason  of  LMSC. 

B.  Transmission  Electron  Microscopy 

Figure  13  contains  some  transmission  electron  micrographs  of  R277  PG 
taken  with  an  Hitachi  HU-11  electron  microscope  operating  at  100  kV.  The 
powder  u.sed  to  obtain  these  samples  came  from  a  section  out  of  a  Nb  diffusion 
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Fig.  7  Photom*»;rographs  Parallel  to  Deposition  Plane  of  11165  F*G.  (al)  As -deposited,  lOOX. 
(a2)  As -deposited,  lOOOX.  (bl)  Heat-treated,  lOOX.  (b2)  Heat-treated,  lOOOX 
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Fig.  8  Photomicrographs  Perpendicular  to  Deposition  Plane  of  R165  PG.  (al)  As-deposited,  lOOX. 
(a2)  As-deposited,  lOOOX.  (bl)  Heat-treated,  lOOX.  (b2)  Heat-treated,  lOOOX 
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Fig.  9  Hiotomicrographs  Parallel  to  Deposition  Plane  of  R277  PG.  (al)  As -deposited,  lOOX. 
(a2)  As-deposited,  lOOOX.  (bl)  Heat-treated,  lOOX.  (b2)  Heat-treated,  lOOOX 
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Fig.  11  Photomicrographs  Parallel  to  Deposition  Plane  of  Rl  PG.  (al)  As-deposited.  lOOX. 
(a2)  As-deposited,  lOOOX.  (bl)  Heat-treated,  lOOX.  (b2)  Heat-treated,  lOOOX 


Fig.  12  Photomicrographs  Perpendicular  to  Deposition  Plane  of  R1  PG.  al)  As-deposited,  lOOX. 
(a2)  As-depvj»sited.  lOOOX.  (bl)  Heat-treated,  lOOX.  (b2)  Heat-treated,  lOOOX 
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couple  which  was  taken  near  the  beginning  of  tlie  composition  profile.  Tlie 
material  was  examined  to  demonstrate  the  nodular,  wrinkled  sheet  structure 
of  the  PG  and  to  show  the  lack  of  evidence  for  clustering  of  Nb  throughout  the 
material.  At  7O,OO0X  magnification,  150  A  diameter  clusters  would  appear 
as  1  mm  in  diameter  and  could  be  easily  resolved.  The  fact  that  clusters  of 
tills  size  or  even  smaller  were  not  seen  in  PG  flakes  obtained  near  the  Ixjgin- 
ning  of  the  concentration  gradients  is  taken  to  indicate  one  of  the  following 
possibilities: 

1 .  The  Nb  is  uniformly  distrilxited  in  the  PG 

2.  The  Nb  is  distrilxited  in  clusters  in  the  PG  on  a  scale  too  fine 
to  Ik*  resolved. 

3.  Tile  Nb  is  concentrated  in  so  few  regions  that  tiK*  chances  of 
observing  it  are  rather  small. 

The  only  indication  of  high  atomic  weight  clusters  is  shown  in  Fig.  13d.  The 
structure  is  not  well  resolved  and  the  nature  of  these  dark  regions  is  not  known. 
Tlie  transmission  electron  micrographs  were  taken  by  Mr.  P.F.  Lindquistof  LMSC . 

C.  X-ray  Structure 

A  summary  of  the  x-ray  data  is  given  in  Tables  9,  10,  11 .  Tlie  as-deposited 
PG's  are  all  characterized  by  a  90%  stacking  error  probability  (p),  a  crystallite 
thickness  (L  )  of  170A,  and  a  crystallite  diameter  (L^)  of  190  to  260  A.  Tlie 
crystallite  thicknesses  were  computed  from  the  0002  profile,  as  outlined  in 
Chapter  V.  Similar  calculations  on  the  0004  and  0006  profiles  resulted  in 
much  smaller  crystallite  thicknesses,  demonstrating  that  tlie'^e  is  a  spread 
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(a)  H.OOOX 


(b)  53.000X 


(c)  70.  OOOX  (d)  70.  OOOX 

FiK.  13  Transmission  Electron  Micrographs  of  H277  PG 

Taken  Near  Beginning  of  Nb  Concentration  Gradients, 
(a)  14,  OOOX,  (b)  53,  OOOX.  tc)  70,  OOOX.  <d)  70,  OOOX 


in  the  intcrplanar  spacings  (60).  The  crystallite  diameters  were  computed 
from  the  1£  and  profiles  by  the  method  given  in  Chapter  VI.  Of  the  cor¬ 
rections  listed  there,  only  the  background,  the  ^ ^^2  ^o^t)let,  and  the  instru¬ 
mental  broadening  corrections  were  important.  The  absorption  correction 
increased  the  peak  height  and  shifted  it  to  higher  20  but  did  not  increase 
the  breadth.  The  preferred  orientation  correction  was  found  to  be  negligible 
for  ( 20,  -20  \  <  0.7  deg  ,  which  was  the  case  for  all  the  PC's 

in  this  study.  The  preferred  orientation  correction  docs  not  become  important 
until  ^20  -  20^^^  j  ~  1  deg  and,  because  the  doublet  correction  is 

appreciable,  the  best  procedure  is  to  correct  for  the  doublet  first  so  that  the 
preferred  orientation  correction  is  minimized. 

The  graphiiizing  heat  treatment  decreased  p  to  30%  for  the  R165  and 
R277  PC’s  and  40%  for  the  R1  material.  Tlie  R1  PG  varied  somewhat  in  its 
response  to  the  heat  treatment  due  to  the  variable  structure  noted  under  metal¬ 
lography.  This  variation  is  also  reflected  in  the  two  different  values  of  p  given 
in  Table  9  for  Rl-IIT  .specimens.  Ilie  crystallite  thicknesses  increased  two  to 
three  times  upon  heat  treatment,  with  the  R165  material  showing  the  greatest 
change  and  Rl  showing  the  least.  Tlie  different  res|>onses  to  the  heat  treatment 
are  also  ilemonstrated  by  tlie  hk  bands  plotted  in  Fig.  14,  where  the  intensities 
of  the  Rl  profiles  are  always  less  than  those  of  the  other  two  PC's  even  though 
the  gc*ometry  in  each  case  is  identical.  It  was  notetl  in  Cliapter  III  that  there 
are  no  well-ileveloped  line  broadening  theories  for  the  modulated  hl^  bands  of 
partially  grapliitized  PCJ;  liowever,  the  JU[  bands  in  all  three  heat-treated  PC's 
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Table  0 


STACriNG  ERROR  PROBABILITIES  FOR  PG  USED  IN  THIS  STUDY 


PG 

Position 

in 

^*0002 

P 

Avg. 

for 

Deposit 

Centroid 

Peak 

Centroid 

Peak 

each  PG 

R1G5-AD 

Top 

3.423 

3.427 

0.94 

R165-AD 

Middle 

3.427 

3.423 

ns 

0.92 

0.9 

R1G5-AD 

Bottom 

3.423 

3.430 

0.92 

0.96 

R1G5-HT 

Top 

3.3GG 

3.360 

0.  32 

0.18 

R1G5-HT 

Bottom 

3.370 

3.366 

0.40 

0.32 

0.3 

R277-AD 

Top 

3.418 

3.416 

0.  89 

0.81 

R277-AD 

Middle 

3.41G 

3.418 

0.81 

0.89 

0.9 

R277-AD 

Bottom 

3.419 

3.422 

0.  89 

0.91 

R277-HT 

Top 

3.372 

3.367 

0.  43 

0.34 

R277-HT 

Middle 

3.367 

3.366 

0.34 

0.32 

0.3 

R277-HT 

Bottom 

3.377 

3.370 

0.52 

0.40 

Rl-AD 

3.435 

3.427 

0.98 

0.94 

0.9 

Rl-HT 

3.383 

3.381 

0.  58 

0. 56 

0.6 

Rl-HT 

3.369 

0.40 

0.4 

H253-32 

Strain- 

nnnealed 

standard 

3.3560 

3.3560 

0.08 

0.08 

t 
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Table  10 

CRYSTALLITE  SIZES  FOR  THE  PC  USED  IN  THIS  STUDY 


PC 

Position 

in 

Deposit 

Thickness  • 

0002  0004 

-L 

c 

0006 

Average 

Thickness 

0002  0004  0006 

Diameter 

10  1 1  Avg. 

R165-AD 

Top 

176A 

81A 

55A 

R165-AD 

Middle 

179 

92 

55 

175 

85 

65 

246 

250 

250 

R1G5-AD 

Bottom 

176 

82 

81 

R165-HT 

Top 

493 

175 

183 

R165-HT 

Bottom 

484 

175 

143 

485 

175 

155 

595 

595 

R277-AD 

Top 

184 

89 

R277-AD 

Middle 

163 

95 

165 

90 

65 

268 

262 

265 

R277-AD 

Bottom 

164 

89 

64 

R277-HT 

Top 

291 

159 

R277-HT 

Middle 

339 

155 

98 

325 

150 

100 

595 

595 

R277-HT 

Bottom 

338 

141 

Rl-AR 

171 

87 

170 

85 

201 

182 

190 

Rl-HT 

215 

121 

_ 

92 

210 

120 

90 

463 

460 
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Fig.  14  UncorreCied  hk  Bands  in  PG 


arc*  rather  sharp,  and  the  line  width  at  half-height  could  be  used  to  approximate 
the  crystallite  diameter.  The  separation  of  the  1 10  and  111  peaks  in  graphite 
is  about  6  deg  20  with  CuKa  radiation,  the  reason  that  the  modulation  effect  is 
not  very  apparent.  The  crystallite  diameters  given  in  Table  10  for  the  heat- 
treated  PC's  were  computed  from  the  ^  bands,  after  correcting  them  for  the 
O-Ort  doublet,  and  instrumental  broadening.  The  preferred  orientation  cor- 
rection  was  not  applied,  and  the  numerical  factor  in  tlie  Scherrer  equation  was 
taken  to  be  1 .84.  As  n.ay  be  seen,  the  crystallite  diameter  apparently  increased 
somewhat  more  than  two-fold  in  all  cases. 

Tlie  results  from  the  preferred  orientation  studies,  with  the  texture  inter¬ 
preted  as  Gaussian  distributions,  are  shown  in  Figs.  15  and  16  and  Table  11. 

The  value  of  m  in  the  cos^^o  representation  is  also  given  in  Table  11 .  In  all 
cases,  the  R1  material  has  the  broadest  texture,  and,  of  the  as-deposited  PG's, 
R277  has  tiie  sharpest  texture.  Upon  graphitization,  there  is  little  difference 
between  R165  and  R277,  but  R1  still  has  the  broadest  texture. 

The  orientation  dependence  of  the  stacking  error  probability  is  given  in 
Fig.  17.  The  as -deposited  material  exhibits  very  little  orientation  dependence 
V'hile  that  for  tlie  heat-treated  PG's  is  somewhat  greater.  The  change  in  p  is 
somewhat  less  than  previously  observed  (5,  94)  but  not  extraordinarily  so.  In 
the  as-deposited  material,  the  line  breadths  indicated  a  decrease  in  crystallite 
thickness  of  about  15%  at  a  *  45  deg  and  a  25%  decrease  for  the  heat-treated 
material.  Although  the  changes  in  cry'stallite  size  and  stacking  error  proba¬ 
bility  are  definitely  present,  they  arc  not  large  enough  to  v  arrant  consideration 
at  this  time. 
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0002  Pole  Figures  for  R165  and  R1  PG  (a  is  angle  away  from 
deposition  plane  normal) 
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Fig.  16  0002  Pole  Figures  for  R277  PG  {o  is  the  angle  away  from  the 
deposition  plane  normal) 
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Table  11 
SUMMARY  OF 

PREFERRED  ORIENTATION  DATA 


PG 

Position 

in 

Deposit 

Mean 

Tilt  Angle 

2  1/2 
<  a  > 

m 

R165-AD 

23.2* 

5.4 

R165-HT 

— 

15.5* 

11. d 

R277-AD 

Top 

19.4* 

7.8 

R277-AD 

Bottom 

19.5* 

R277-HT 

Top 

16.7* 

9.2 

R277-HT 

Bottom 

15.5* 

Rl-AD 

- 

27.5* 

3.6 

Rl-HT 

— 

20.6* 

7.2 

lOI 


STACKING  ERROR  PROBABIUTY, 


Fig.  17  Orientation  Dependence  of  Stacking  Error  Probability 
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D.  Small  Angle  Scattering 

Small  angle  scattering  of  x-rays  by  solids  or  liquids  is  caused  by  hetero- 
genieties  in  the  electron  density  with  dimensions  from  several  tens  to  several 
hundred  times  the  x-ray  wave  length  (126, 127).  These  heterogenieties  may  be 
regions  of  high  or  low  electron  density;  for  example,  in  solids,  the  hetero¬ 
genieties  could  Ik*  voids.  X-rays  are  scattered  out  of  the  direct  beam  by  less 
than  3  deg  20  and  careful  collimation  is  required  to  allow  data  to  be  taken  as 
close  as  0.1  deg  20  away  from  the  beam.  Small  angle  scattering  data  are 
generally  described  in  terms  of  the  diffraction  vector  h  =  (47r  sin6)/X  . 

This  h  should  not  be  confused  with  the  h  used  in  the  Miller  indices  notation. 

Guinier  (126)  has  shown  that  for  small  values  of  h,  small  angle  scattering 
from  a  collection  of  fixed  identical  particles  having  the  same  orientation, 
behaves  as 

I  =  kNn^  exp  (-h^D^) 

where  I  is  the  intensity  of  the  scatter  x-rays,  k  a  constant,  N  the  total 
number  of  particles,  n  the  number  of  electrons  per  particle,  and  D  the 
inertial  distance  of  the  particles  relative  to  the  plane  of  the  incident  beam. 

This  relationship  is  termed  the  Guinier  approximation  and  is  valid  for  hD  <  1  . 
For  ellipsoids  of  revolution  (126) 

=  b^/5 

wlxjre  b  is  the  particle  dimension  perpendicular  to  the  plane  of  the  beam. 
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For  x-ray  lH*ams  possessing  neither  height  nor  width,  the  small  angle 
scattering  approaches  an  h"'*  dependency  at  large  h  and  a  beam  of  Infinite 

_3 

height  ciisiorts  this  asymptotic  behavior  to  an  h  dependency  (126, 127). 

Bragg  (128)  has  shown  that,  for  freely  moving  mono-disperse  particles,  the 

-4  -3 

onset  of  the  asymptotic  h  or  h  dependency  occurs  at 

hq  »  4.0 

where  q  is  the  smallest  dimension  of  the  particles.  This  approximation, 
which  will  be  referred  to  as  the  Bragg  approximation,  can  be  extended  to  fixed 
particles  having  the  same  orientation  (129)  as 


The  data  at  intermediate  h  generally  do  not  lend  themselves  to  analysis; 
the  equations  given  above,  however,  do  provide  two  independent  estimates  of  the 
particle  size,  one  at  small  h  and  the  other  at  large  h  . 

Small  angle  scattering  of  x-rays  can  bj  complicated  by  the  presence  of 
doubly  diffracted  x-rays.  This  distortion  is  termed  double  Bragg  scattering 
(127, 130).  If  the  diffracting  planes  are  aligned  to  give  strong  double  Bragg 
diffraction,  the  initial  portion  of  the  small  angle  scattering  data  can  be  obscured. 

Using  the  procedure  described  in  Chapter  VI,  the  small  angle  scattering  of 
x-rays  from  as -deposited  and  heat-treated  PC  was  measured  with  the  plane  of 
the  beam  both  parallel  and  perpendicular  to  the  deposition  plane.  The  experi¬ 
mental  configuration  is  shown  schematically  in  Fig.  18,  together  with  the  data 
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Fig.  18  Small  Angle  Scattering  by  R165  PG.  Curves  (al)  and  (bl):  As-Deposited  PG; 
Curves  (a2)  and  (b2):  Heat-treated  PG 
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plotted  as  log  I  vs  log  h.  Because  the  density  of  PC  is  nearly  equal  to  the 
theoretical  x-ray  density,  it  is  assumed  that  the  scattering  is  caused  by  the() 
presence  of  small  voids  having  essentially  zero  electron  density  (95).  The 
asymptotic  h  dependency  is  observed  for  all  the  data,  but  the  data  for  which 
the  x-ray  beam  was  parallel  to  the  deposition  plane  do  not  obey  the  Guinicr 
approximation.  In  this  plot,  the  slope  of  the  data  should  decrease  at  small  h, 
and  curves  for  specimens  of  the  same  condition  but  different  orientation  should 
approach  the  same  value  of  I  as  h  becomes  very  small.  This  deviation  from 
theory  for  specimens  with  the  beam  oriented  parallel  to  the  deposition  plane  is 
the  result  of  double  Bragg  scattering  from  the  0002  planes  of  the  PG.  As  dis¬ 
cussed  by  Bragg  et  al .  (95),  the  layer  planes  in  PG  are  aligned  to  give  strong 
double  Bragg  diffraction  in  this  specimen  configuration.  This  distortion  is  so 
great  here  that  it  does  not  permit  the  application  of  the  Guinier  approximation, 
but  the  void  thickness  perpendicular  to  the  deposition  plane  can  be  estimated 

_3 

from  the  onset  of  the  h  dependency . 

2 

In  Fig.  19,  the  small  angle  scattering  data  are  plotted  as  log  I  vs  h  to 
demonstrate  the  Guinier  approximation.  In  this  plot,  the  data  for  specimens 
of  the  same  condition  but  different  orientation  should  extrapolate  to  the  same 
value  of  I  at  h  =  0  because  the  scattering  is  caused  by  the  same  void  volume 
in  both  cases .  The  distortion  produced  by  the  double  Bragg  diffraction  prevents 
the  curves  noted  by  (a)  from  being  properly  extrapolated. 

Two  dimensions  are  required  to  characterize  the  voids  observed  by' the 
small  angle  scattering:  a  void  thickness  perpendicular  to  the  deposition  surface, 
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X  10^ 

Fig.  19  Guinier  Analyiii  of  Small  Angle  Scattering  Curvea.  Curves  (al)  and  (bl): 
As -Deposited  PG;  Curves  (a2)  and  (b2):  Heat-treated  PG 
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,  and  a  void  diameter  parallel  to  the  deposition  surface,  .  Estimates  of 
these  parameters  as  determined  by  applying  the  Guinier  and  modified  Bragg 
approximations  to  the  data  given  in  Figs.  18  and  19,  are  listed  below: 


Method  of 

As-Deposited 

Heat-Treated 

Estimating 

b 

c 

a  c 

Guinier 

158A 

— 

390A 

Modified  Bragg 

120 

50 

240  60 

The  agreement  between  the  two  estimates  is  fairly  consistent,  showing  that 
during  the  graphitizing  heat -treatment  tne  void  diameter  increased  by  about  a 
factor  of  2  and  the  thickness  remained  essentially  unchanged.  Tbis  is  also  con¬ 
sistent  with  tl  e  increase  in  crystallite  diameter  with  heat-treatment  as  given  in 
Table  10.  The  void  diameter  is  somewhat  smaller  than  the  crystallite  size 
obtained  from  the  line  broadening  experiments;  however,  it  should  be  noted 
that  the  theories  on  which  both  these  estimates  are  based  are  not  accurate  within 
a  factor  of  2,  but  they  do  show  proper  relative  sizes.  Similar  comments  apply 
to  comparisons  between  the  Guinier  and  modified  Bragg  approximations  in  the 
small  angle  scattering  data . 

A  more  detailed  interpretation  of  small  angle  scattering  in  PG  is  being 
developed  by  Bragg  and  co-workers  (129).  As  yet,  this  work  has  not  been 
completed.  Nevertheless,  the  results  reported  here  will  probably  agree  within 
a  factor  of  2  with  the  results  of  the  more  detailed  analysis. 

Mr.  J.C.  Robinson  of  LMSC  assisted  in  obtaining  these  data. 
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E.  Other  Structural  Parameters 


The  mercury  porosimetry  data  showed  that  the  volume  of  accessible  pores 
in  PG  (in  either  the  as-deposited  or  heat-treated  condition)  which  can  be  pene- 
trated  with  mercury  at  pressures  up  to  5000  psi  is  less  than  0.001  cm  /cm  of 
PG.  This  data  represents  the  limit  of  sensitivity  for  the  porosimeter  for  the 
samples  used.  The  results  were  not  unexpected,  when  one  considers  the  high 
density  of  PG  compared  with  its  x-ray  density  and  the  lack  of  porosity  in  the 
microstructure.  The  porosimetry  data  were  obtained  by  Mr.  K.A.  Till  of  LMSC. 

The  bulk  densities  of  the  various  PG's,  as  determined  by  hydrostatic  or 
differential  weighing  in  isopropyl  alcohol,  are  given  in  Table  12,  together 
with  the  calculated  x-ray  densities.  The  density  values  are  about  0.003  gm/ 
cm  too  high  because  the  air  bouyancy  correction  was  not  applied;  however, 

3 

the  accuracy  of  the  determination  is  estimated  to  be  ±0.005  gm/cm  with  the 
major  error  occurring  in  the  surface  tension  effects  on  the  supporting  wire. ' 

As  expected,  there  is  a  small  void  volume  in  the  as-deposited  PG's  (typically 
0.4%),  but  the  graphitizing  heat -treatment  virtually  eliminated  the  porosity 
except  in  the  case  of  the  R1  material,  where  a  slight  porosity  persists. 

The  bulk  density  data  were  obtained  by  Mr .  R .  Nygren  of  LMSC . 

F .  Summary 

The  PG's  used  in  this  investigation  are  representative  of  large  (R277)  and 
small  (R165)  cone -sized  PG  of  the  substrate  nucleated  type  and  of  the  small 
cone-sized,  regeneratively  nucleated  type  (Rl).  The  response  of  each  PG  to 
heat -treatment  is  also  typical,  the  major  features  being  an  increase  in 


109 


Table  12 


COMPARISON  OF  OBSERVED  DENSITIES  AND 
DENSITIES  COMPUTED  FROM  OBSERVED  INTERLAYER  SPACINGS 


PG 

Position 

in 

Deposit 

Observed 

Density 

(gm/cm^) 

Average 

Observed 

Densitv 

(gm/cm^ 

**0002 

A 

X-ray 

Density^*) 

(Px) 

Estimated 

Void 

Volume 

(%) 

R165-AR 

2.  212 

R165-AR 

2.206 

2.214 

3.425 

2.220 

0.3 

R165-AR 

2.223 

R165-HT 

Bottom 

2.  263 

R165-HT 

Bulk 

2.265 

2.264 

3.366 

2.259 

0.0 

R165-HT 

Bulk 

2.  263 

R277-AR 

Top 

2.219 

R277-AR 

Bulk 

2.219 

2.212 

3.417 

2.225 

0.6 

R277-AR 

Bulk 

2.206 

R277-HT 

Middle 

2.260 

R277-HT 

Bulk 

2.254 

2.263 

3.366 

2.259 

0.0 

R277-HT 

Bulk 

2.  265 

Rl-AR 

Bulk 

2.205 

Rl-AR 

Bulk 

2.217 

2.205 

3.431 

2.216 

0.5 

Rl-AR 

Bulk 

2.193 

Rl-HT 

Bulk 

2.251 

2.250 

3.375 

2.253 

0.1 

Rl-HT 

Bulk 

2.249 

(a)Px  =  7.603/d^jg2  gm/cra^ 
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crystallite  diameter  and  thickness  of  about  2.5  times  and  a  decrease  in  stacking 
error  probability  from  90%  to  approximately  35%.  This  was  accompanied  by  a 
decrease  in  mean  tilt  angle  from  somewhat  more  than  20  deg  to  about  15  deg 
and  virtually  complete  elimination  of  void  volume.  Very  little  delamination 
occurred  during  the  graphitizing  heat -treatment. 


Ill 


Chapter  VIII 


EXPERIMENTAL  RESULTS  “  DIFFUSION 

A.  General  Behavior  of  Data 

The  structure  and  many  of  the  properties  of  PG  are  anisotropic  and,  as 
expected,  the  diffusion  of  metals  in  PG  is  also  anisotropic.  For  the  same 
thermal  treatment,  the  diffusion  rate  of  Nb  parallel  to  the  deposition  plane 
is  approximately  30  times  that  of  Nb  perpendicular  to  the  deposition  plane. 

This  factor  is  not  the  same  at  every  temperature  because  the  activation  ener¬ 
gies  for  diffusion  in  the  two  directions  are  not  equal .  The  anisotropy  of  metal 
diffusion  in  PG  gave  rise  to  two  different  experimental  situations. 

The  composition  profiles  fell  into  two  general  categories:  those  in  which 
considerable  activity  remained  on  the  surface  and  in  the  first  three  laps,  and 
those  in  which  there  was  very  little  residual  surface  activity,  ^ecimens  in 
which  diffusion  took  place  perpendicular  to  the  deposition  plane  had  high  residual 
surface  activities  while  those  for  diffusion  parallel  to  the  deposition  plane 
usually  did  not  have  a  high  residual  surface  activity.  Approximately  the  same 
amount  of  Nb  was  placed  on  each  sample,  so  these  results  are  consistent  with 
the  much  higher  D  values  observed  parallel  to  the  deposition  plane .  T\vo 
typical  profiles  for  the  two  different  directions  are  shown  in  Fig.  20.  The  dif¬ 
ference  in  behavior  of  the  two  kinds  of  samples  led  to  the  conclusion  that  two 
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Fig.  20  Typical  Composition  Profiles  Perpendicular  and  Parallel  to 
Deposition  Plane 
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different  solutions  to  the  diffusion  equation  would  be  required  to  treat  the  data 
adequately . 

Tne  concentrations  observed  in  these  experiments  were  very  low .  For  the 
e:q)eriments  approximating  constant  surface  composition  conditions,  the  activity 

3 

concentration  at  the  surface  was  about  0 . 1  mc/cm  .  With  a  counting  efficiency  of 
approximately  25%,  this  observed  activity  corresponds  to  a  true  specific  activity  of 

3 

0.4  mc/cm  PG.  For  Nb-95,  the  carrier-free  specific  activity  without  absorp¬ 
tion,  is  about  4  x  10^  mc/gm  Nb.  From  these  data,  the  maximum  concentration 

-7  3 

of  Nb  in  PG  was  about  10  gm  Nb/cm  PG.  As  will  be  shown  in  Chapter  IX,  this 
concentration  does  not  represent  the  true  solubility  of  Nb  in  PG,  but  rather, 
represents  a  crystallite  edge  decoration. 

B.  Computation  of  Diffusion  Coefficients 

As  demonstrated  in  Fig.  20,  the  data  were  sufficiently  precise  to  distinguish 
the  two  separate  categories.  It  was  felt  that  these  categories  were  the  result  of 
having  two  situations  prevail  during  the  diffusion  anneals:  the  constant  surface 
composition  situation  (CSC)  and  the  instantaneous  source  situation  (IS)  for  semi- 
infinite  media(131, 132).  In  no  case  did  the  data  suggest  the  linear  log  C  vs  x  depen¬ 
dence  predicted  by  various  grain  boundary  diffusion  theories.  This  latter 
observation  will  be  discussed  in  more  detail  in  the  next  chapter. 

For  the  experimental  conditions  of  a  salt  of  a  metal  deposited  on  the 
specimen  surface,  it  is  not  possible  to  predict  the  exact  nature  of  the  boundary 
conditions  that  apply  during  the  diffusion  anneal . 
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If  the  solubility  and  diffusivity  are  extremely  low,  then  even  the  small 
amount  of  tracer  deposited  on  the  surface  will  not  all  diffuse,  and  a  constant 
surface  composition  situation  would  be  approximated .  On  the  other  hand,  higher 
diffusivity  can  result  in  all  the  tracer  diffusing  into  the  PG  and,  if  this  occurs 
near  the  beginning  of  the  diffusion  anneal,  the  instantaneous  source  will  be 
approximated.  An  intermediate  situation  can  arise  if  the  surface  composition 
is  neither  constant  nor  vanishing  near  time  zero,  but  is  a  variable  function  of 
time.  If  an  independent  estimate  of  the  surface  composition  is  available,  then 
it  is  possible  to  account  for  this  last  complication  (133);  however,  this  situa¬ 
tion  did  not  prevail  in  these  PG  diffusion  experiments.  As  a  result,  only  the 
two  extreme  situations  were  examined . 

If  the  CSC  solution  applies,  then  the  composition  relative  to  the  surface 

composition  (C/C  )  should  behave  as 
s 

a 

C  2  r  2  X 

^  =  1  -  “  I  exp  (-a  da)  *  erfe  a  =  erfc 

*  o 

and  if  ..he  IS  solution  applies,  then 

2 

C  =  k  exp  (-a  )  =  k  exp 

will  represent  the  C-x  behavior  (131, 132).  These  two  functions  are  plotted  in 

Fig.  21 .  For  the  IS  equation,  k  was  taken  arbitrarily  to  be  2;  for  the  CSC 

2 

equation,  log  (C/2C  )  is  plotted  against  a  using  1/2  erfc  a  for  the  dependency . 
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Fig.  21  Behavior  of  Inatantaneoua  Source  (a)  and  Constant  Surface  Composition 
(b)  Solutions  to  Pick's  Diffusion  Laws  for  Semi-infinite  Solid 


As  expected,  the  IS  solution  is  linear  over  the  entire  range,  but  the  CSC  solu¬ 
tion  is  approximately  linear  only  beyond  C/2Cg  ^  0.08  and  the  slope  of  the 


line  increases  with  smaller  a  .  Clearly,  if  one  assumed  that  an  IS  solution 
applied  in  an  experiment  in  which  the  CSC  solution  were  more  appropriate, 
and,  if  the  diffusivity  were  computed  from  the  slope  of  the  curve,  some  seri¬ 
ous  error  could  result  for  a  <  l.O  .  Malkovich  (134)  has  treated  this  problem 
analytically.  He  derives  the  following  expression,  relating  the  0  measured  by 
assuming  an  IS  solution  with  the  true  D  appropriate  to  the  CSC  solu- 

(°actual)= 

“meaa  '  <“^)  «rt<=  “l^actual 

meas  ~  ^  ^actual 
“  ’ 

He  also  computed  the  \  factor: 


a 

X 

0.10 

0.16 

0.25 

0.34 

0.50 

0.55 

1.00 

0.76 

2.00 

0.91 

2  2 

Actually,  the  curve  of  log  C/2Cg  vs.  o  or  x  for  the  CSC  solution  is  every¬ 
where  concave  upward;  empirically,  one  can  choose  a  slope  starting  at  a  =  1.0 
that  gives  to  be  0.87  ,  as  depicted  in  Fig.  21 . 
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Figure  21  definitely  shows  that  experimental  conditions  approximating 
the  CSC  solution  can  be  distinguished  from  the  IS  solution  by  plotting  the  data 

9 

as  log  C  vs  X  ,  and  the  data  should  be  treated  accordingly.  In  this  study,  all 
the  data  were  first  plotted  as  log  C  vs  x  .  Those  curves  which  were  linear  for 
C/C  <0.9  were  treated  as  approximating  the  IS  conditions  while  those  which 
curved  upward  for  C/C^  >0.1  were  treated  as  approximating  the  CSC  condi¬ 
tions.  For  the  CSC  samples,  the  first  three  or  four  laps  showed  very  high 
activity  relative  to  the  estimated  value  of  C  ;  this  additional  activity  was 
attributed  to  radioactive  material  remaining  on  the  surface  and  the  lack  of 
perfect  interface  resolution. 

For  the  samples  to  be  treated  as  approximating  the  IS  conditions,  the  dif- 

fusivities  were  computed  from  the  average  of  the  maximum  and  minimum  values 

obtained  from  lines  drawn  with  maximum  and  minimum  slopes  through  the  data 

2 

plotted  as  log  C  vs  x  .  A  standard  deviation  for  each  diffusivity  was  estimated 

by  dividing  the  difference  between  the  maximum  and  minimum  values  by  /2 

These  estimated  standard  deviations  were  between  10  and  20%  of  the  value 

measured  and  occasionally  were  as  high  as  30%. 

To  obtain  a  first  estimate  of  C  ,,  the  data  approximating  the  CSC  conditions 

2  2 

were  extrapolated  back  to  x  =  0  on  the  log  C  vs  x  plots.  Each  composition 
point  was  then  divided  by  2C„  and  the  result  plotted  vs  x^  on  the  lower  half  of 
linear  probability  paper  (Keuffel  and  Esser  Co.,  #359-23)  the  scale  of  which 
is  designed  to  make  1/2  erfe  a  linear  from  0.50  to  0.0001 .  If  the  data  did  not 
extrapolate  back  to  C/2Cg  =  0.50,  then  new  values  of  Cg  were  selected  and 
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1 

* 

each  point  again  adjusted  until  a  straight  line  through  the  data  would  intersect 

the  probability  scale  at  0 . 50  relative  composition .  It  should  be  noted  here  that  • 

^  as  long  as  the  C  value  is  incorrect,  the  data  will  not  lie  on  a  straight  line  but 

s 

will  be  slightly  curved  in  the  direction  of  the  error .  After  a  suitable  value  of 
was  determined,  lines  of  maximum  and  minimum  slope  were  drawn  to  com¬ 
pute  maximum  and  minimum  values  of  D  .  The  diffusivities  were  computed 
from  the  equation 


X 

2^^ 


using  a  single  value  of  C/2C  ,  usually  0.00150.  Average  diffusivities  and 

s 

standard  deviations  were  estimated  as  above. 


C .  Measured  Diffusion  Coefficients 

Table  13  lists  the  diffusivities,  as  determined  above,  together  with  the 
PG  type,  time  and  temperature  of  anneal,  and  groups  the  data  according  to 
diffusion  anneal .  Where  the  word  fail  is  written,  the  data  were  not  used  to 
calculate  a  D  because,  upon  lapping  and  cleaning  the  specimen  with  mylar 
tape,  the  surfaces  became  extremely  uneven  and  produced  very  erratic  pro¬ 
files.  This  damage  to  the  samples  was  prevented  later  by  cleaning  the  heat- 
treated  specimens  with  cellophane  tape  which  is  less  adherent  and  will  not 
tear  the  heat-treated  PG .  The  symbol  (d)  in  the  table  signifies  that  the 
data  are  probably  influenced  by  diffusion  parallel  to  the  deposition  plane . 

The  symbol  <axlenotes  far  greater  penetration  than  expected,  indicating  that 
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Table  13 

DIFFUSION  COEFFICIENTS  (D)  FOR  Nb-95  IN  VARIOUS  TYPES  OF  PG 
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(c)  Md  M>  aolNtioB  iiaed 

(d)  Onta  poMttdy  iafliWBced  by  <a>  dlffysioa 
(c)  <a>  diffiiaioa  doniiuded 


diffusion  parallel  to  the  deposition  plane  definitely  dominated  the  diffusion  process . 
Finally,  the  letter  (c)  denotes  a  D  value  obtained  with  a  Mb  solution  112  days  (3.1 
half-lives)  older  than  the  solutionusedontheother  samples  of  that  annela.  Nb-95 
decays  to  Mo-95  and  this  comparison  was  intended  to  observe  any  effect  the  Mo 
might  have.  As  may  be  seen,  the  Mo  did  not  have  any  discernable  effect. 

Most  of  the  data  were  reduced  by  the  author's  wife,  Betty. 

D.  Temperature  Dependence  of  Diffusion  Coefficients 

Figures  22  and  23  are  Arrhenius  plots  of  the  diffusion  coefficients  grouped 
according  to  the  two  directions  studied.  The  data  for  each  PG,  in  each  condi¬ 
tion  and  direction,  were  fitted  to  the  equation 

D  =  D^expC-Q/RT) 

using  a  leafit  squares  technique  developed  by  Deming  (135)  which  permits  one 
to  include  the  estimated  errors  of  each  data  point  independently.  The  data  were 
fitted  with  the  logarithm  form  of  the  Arrhenius  equation  (the  standard  deviations 
for  D^  computed  in  this  form  are  factors  rather  than  additive  terms).  The 
result  of  these  computations  is  given  in  Table  14 .  Because  the  activation 
energies  for  a  given  direction  were  all  within  their  overall  standard  deviations, 
they  were  averaged  to  produce  one  activation  energy  for  diffusion  parallel  to 
the  deposition  planes,  ,  and  another  perp>endicular  to  the  deposition  planes, 
Qc  •  ^  were  then  recomputed  using  the  average  activation  energies  in 

each  direction .  The  results  of  these  calculations  are  also  given  in  Table  14 . 
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Fig.  23  Arrfaeniiu  Plot!  of  Mb  DUbulon  CoefOcioiit  in  PO 
Perpendicular  to  Depoeition  Plane 
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Table  14 

TEMPERATURE  DEPENDENCE  OF  DIFFUSION  COEFFICIENTS 
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(a)  Standard  deviation  of  D 


Using  the  propagation  of  error  analysis  (120)  it  can  be  shown  (121)  that 


fa  _  ^ 

q2  -  V  Q  ;  p2  ^2 


where  is  the  standard  deviation  of  any  quantity  A  ,  T  the  temperature 
in  °K,  Q  the  activation  energy,  R  the  gas  constant,  and  the  pre- 
exponential  factor  in  the  Arrhenius  equation.  As  determined  by  this  study 
and  that  of  WMB  (17),  the  diffusion  coefficients  of  metals  in  PG  exhibit  very 
high  activation  energies  and  Q/R  is  typically  50  x  10^  °K.  With  the  tem¬ 
perature  range  in  the  vicinity  of  2000°K,  the  factor  Q/RT  is  about  25. 
Reference  to  the  above  equations  shows  that  this  high  value  of  Q/RT  makes 
the  determination  of  Q  rather  insensitive  to  errors  in  D  while  the  p re - 
exponential  factor  is.  extremely  sensitive  to  errors  in  T  .  For  example, 
if  <7£j/D  =  20%  and  Or^l%  =  1.5%,  then 

2 

— ^  -  (.  2of  +  (25)^  (.  015)^ 
o 


and 


0.42 
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The  values  of  and  given  above  are  typical  of  the  conditions 

used  in  these  experiments;  however,  the  resulting  values  of  /D^  from  the 

o 

least  squares  computation  are  much  larger  than  the  value  estimated  above. 

Table  14  indicates  that  the  D^'s  for  a  particular  direction  and  condition  do 

not  differ  appreciably  for  all  the  different  types  of  PG .  The  small  variation  of 

with  microstructure  is  taken  to  indicate  that  the  additional  scatter  in 
o  o 

is  not  due  to  structural  variations  within  a  particular  PG  type.  Instead,  the 
additional  scatter  is  probably  due  to  errors  in  the  temperature  determination. 
Besides  the  sources  of  temperature  error  listed  in  Chapter  VIII,  which  totaled 
less  than  1.5%,  additional  errors  could  arise  due  to  unknown  temperature  gradi¬ 
ents  over  the  stack  of  specimens  and  within  the  holder.  Another  source  of 
error  is  the  possibility  of  temperature  differences  between  the  black  body  hole 
in  the  graphite  holder  and  the  specimens.  The  size  of  these  potential  errors 
is  not  known,  but  they  would  not  have  to  be  very  large  to  appreciably  increase 
the  overall  scatter  in  the  data . 

E .  Autoradiography 

Some  representative  autoradiographs  are  shown  in  Figs.  24,  25,  and  26, 
together  with  photomicrographs  of  the  same  specimens.  These  autoradiographs 
were  taken  at  various  depths  in  the  concentration  profile  and  their  position  in 
the  profile  was  not  a  factor  in  their  appearance.  Figures  24  and  25  contain 
autoradiographs  taken  parallel  to  the  deposition  plane  and  show  how  the  Nb  is 
concentrated  around  the  circles  which  represent  the  major  cone  boundaries  in 
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R165 


R1 

Fig.  24  Direct  Comparison  of  Structure  and  Autoradiographs  Parallel  to 
Deposition  Surface  for  R165  and  R1  As-deposited  PG 


\ 
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As  -deposited 


Heat-treated 

Fig.  25  Direct  Comparison  of  Structure  and  Autoradiographs  Parallel  to 
Deposition  Surface  for  R277  PG,  As-deposited  and  Heat-treated 
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the  PG.  The  R1  autoradiograph  (Fig.  24)  was  heavily  exposed  to  determine  if 
the  centers  of  the  larger  circles  contained  any  appreciable  radioactivity.  As 
is  evident,  the  larger  circles  are  virtually  devoid  of  any  radioactivity  and  the 
areas  containing  many  smaller  cones  show  nearly  complete  film  darkening. 
Similar  effects  are  present  in  the  other  views  of  the  deposition  planes  of  the 
R165  and  R277  PG.  WMB  (17)  obtained  essentially  the  same  autoradiographs. 
They  further  observed  that  autoradiographs  taken  from  samples  in  which  dif¬ 
fusion  occurred  near  the  substrate  surface  of  the  deposit  showed  smaller 
circles  than  did  samples  for  diffusion  near  the  final  deposition  surface.  The 
major  diameter  of  the  circles  increases  from  the  bottom  to  the  top,  as  was 
clearly  delineated  in  the  WMB  autoradiographs.  All  of  these  autoradiographs 
are  proof  that  metal  diffusion  in  PG  perpendicular  to  the  deposition  plane  occurs 
down  primary  cone  boundaries,  with  entry  into  the  cones  taking  place  princi¬ 
pally  from  these  primary  cone  boundaries. 

In  Fig.  25,  it  is  evident  that  the  lieat -treatment  had  very  little  effect  on 
the  manner  in  which  the  Nb-95  is  distributed.  This  is  to  be  expected  because 
the  heat -treatment  did  not  appreciably  alter  the  primary  structure.  The  in¬ 
creased  detail  in  the  R277-HT  autoradiograph  is  the  result  of  repolishing  the 
sample  after  lapping  to  be  certain  that  the  observed  distribution  of  Nb-95  was 
not  an  artifact.  Autoradiographs  and  microstmctures  of  R277  as-deposited 
and  heat-treated  specimens  viewed  on  the  plane  perpendicular  to  the  deposition 
plane  are  shown  in  Fig.  26.  The  major  features  to  be  noted  here  are  the  lack 
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cone  delineation  and  the  presence  of  streaks  parallel  to  the  deposition  plane . 
These  streaks  are  due  to  occasional  delaminations  or  are  reflections  of  the 
PC's  tendency  to  delaminate.  They  demonstrate  the  diffusing  Nb  will  decorate 
places  of  delamination  or  incipient  delamination.  Both  of  these  features  will 
be  discussed  in  more  detail  in  Chapter  IX. 
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Chapter  iX 


DISCUSSION 


A .  Summary  of  Results 

The  diffusion  coefficients  of  Nb  in  several  different  types  of  PG  have  been 
determined  both  parallel  and  perpendicular  to  the  deposition  plane,  and  the 
structure  of  each  PG  has  been  characterized  by  the  application  of  several  tech¬ 
niques  of  metallography.  x-ray  diffraction,  and  density  determinations.  The 
structure  and  diffusion  data,  reported  in  detail  in  Chapters  VII  and  VIII,  are 
summarized  in  Table  IS.  These  dau  will  be  reviewed  briefly  and  then  a  model 
for  the  diffusion  of  metals  in  PG  will  be  developed.  Following  these  discussions, 
the  model  will  be  applied  to  other  forms  of  graphite  and  diffusion  alloying  PG 
will  be  discussed  briefly . 

The  as-deposited  PG's  represented  a  wide  variety  of  microstructures. 

Their  structures  determined  by  x-ray  techniques,  however,  were  quite  similar. 
The  following  structure  parameters  were  derived  from  the  x-ray  diffraction 
studies: 

e  Stacking  error  probability  and  x-ray  density 

•  Crystallite  thickness  and  diameter 

•  Degree  of  preferred  orientation  and  moan  tilt  angle 

a  Void  thickness  and  diameter. 
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The  major  differences  in  the  x-ray  structure  occurred  in  the  degree  of  pre¬ 
ferred  orientation.  Some  of  the  specimens  were  given  a  graphitizing  heat- 
treatment  which  increased  their  density  and  crystallite  size,  sharpened  their 
texture,  decreased  their  stacking  error  probability,  and  altered  their  void 
structure.  The  specific  values  relative  to  these  changes  are  summarized  in 
Table  IS. 

As  noted  in  Chapter  Vlll,  all  the  diffusion  coefficients  in  a  particular 

direction,  regardless  of  PG  type,  could  be  represented  by  a  single  activation 

energy  and  a  variable  pre -exponential  factor.  The  scatter  of  the  data  were  about 

a  factor  of  2.5  and  somewhat  less  for  the  R165  material.  The  graphitizing  heat- 

treatment  increased  D  _  ,  decreased  D  ,  did  not  particularly  alter  the 

O  *  o  c 

spread  in  the  data,  and  diu  not  seem  to  affect  the  activation  energy  for  diffusion 
in  either  direction. 

Autoradiographs  taken  perpendicular  to  the  mass  flux  for  diffusion  perpen¬ 
dicular  to  the  deposition  plane  (D^_^  )  demonstrated  that  diffusion  perpendicu¬ 
lar  to  the  deposition  plane  occurs  in  the  primary  or  major  cone  boundaries . 
Autoradiographs  taken  perpendicular  to  the  deposition  plane  for  diffusion  paral¬ 
lel  to  the  deposition  plane  showed  a  uniform  distribution  of  radioactivity.  In 
this  configuration,  the  cone  structure  was  not  delineated  by  the  diffusing  Nb  but 
some  streaks  occurred  parallel  to  the  deposition  plane  in  the  autoradiographs. 

B.  Correlation  with  Structure 

The  rather  small  variation  of  diffusion  coefficients  for  a  wide  range  of 
microstructures  indicates  that  the  structure  of  the  PG,  as  determined  by  optical 
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microscopy,  has  little  effect  on  the  diffusivity  of  metal  atoms  in  PG .  However, 

the  autoradiographs  demonstrated  that  the  PG  microstructure  determines  the 
distribution  of  the  diffusing  species  and  the  number  of  equivalent  diffusion 
paths  per  unit  area.  Thus,  the  total  mass  flow  of  metals  diffusing  in  PG  will 
be  controlled  by  the  microstructure,  but  the  diffusivity,  defined  from  Pick's 
First  Law  as  the  total  mass  flow  per  unit  concentration  gradient  and  unit 
sample  area,  will  not  be  affected  by  the  microstructurc. 

This  diffusivity  docs  not  depend  upon  the  microstructurc  because  both  the 
total  mass  flow  and  the  concentration  per  unit  sample  volume  are  proportional 
to  the  number  of  equivalent  paths  per  unit  sample  area.  The  statement  can  be 
made  because,  as  will  be  shown  in  this  chapter,  true  volume  of  metals  in  PG 
is  inconsequential,  and  the  diffusing  atoms  arc  restricted  to  the  crystallite 
surfaces.  Thus,  the  situation  is  approximately  equivalent  to  studying  metal 
diffusion  in  a  metal  mesh;  the  measured  diffusion  coefficient  does  not  depend 
upon  the  mesh  size. 

The  somewhat  lower  scatter  in  the  diffusion  data  for  the  R165  material  is 
explained  by  its  more  uniform  microstructurc  because  the  uniformity  of  the 
substructure  is  probably  related  to  the  uniformity  of  the  microstructure.  As 
described  in  Chapter  VIII,  some  of  the  scatter  in  the  diffusion  data  is  definitely 
the  result  of  the  very  high  activation  energies;  even  a  1%  error  in  temperature 
can  appreciably  affect  the  errors  in  the  nre -exponential  factor  in  the  Arrhenius 
equation.  Other  contributions  to  the  scatter  could  arise  from  the  presence  of 
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delaminations,  slight  misorientations  of  the  samples  relative  to  the  deposition 
plane,  and  structural  variations  among  the  samples  of  a  particular  PG. 

Graphitization  did  not  appear  to  change  the  general  modes  of  diffusion 
because  neither  the  autoradiographs  nor  the  activation  energies  for  a  particular 
direction  were  altered  appreciably.  Since  the  activation  energies  were  not 
affected  by  graphitization,  changes  in  the  pre -exponential  factor,  ,  must 
be  responsib'e  for  the  changes  in  the  diffusion  coefficients  with  graphitization . 
The  graphitizing  heat -treatment  increased  the  crystallite  diameter  by  about  a 
factor  of  2.4  and  increased  about  5.2  times .  Because  is  propor¬ 

tional  to  the  distance  squared  (136),  and  because  the  square  of  the  increase  in 
the  crystallite  diameter  is  5.8,  the  jump  distance  parallel  to  the  deposition 
plane  must  be  related  to  the  crystallite  diameter.  This  unusual  correlation 
will  be  discussed  later  in  more  detail.  In  contrast  to  this  behavior,  D 

o-c 

decreased  with  graphitization  even  though  the  crystallite  thickness  increased 
about  two -fold. 

The  observed  5  to  10  deg  decrease  in  mean  tilt  angle  with  heat -treatment 
would  not  appreciably  alter  the  nature  of  length  of  the  diffusion  paths  by  itself. 
This  change,  however,  accompanies  the  many  other  changes  in  structure 
brought  about  by  graphitization,  and  the  diffusion  path  is  thus  indirectly 
related  to  the  texture.  For  example,  a  sharpening  of  the  texture  would  de¬ 
crease  the  width  of  the  crystallite  boundaries  by  decreasing  the  angular  mis- 
orientation  across  them. 
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C .  Possible  Mechanisms  for  Metal  Diffusion  in  PC 


Possible  diffusion  mechanisms  for  metals  in  PG  can  be  separated  into  two 
categories  —  mechanisms  within  the  crystallites  and  mechanisms  external  to 
the  crystallites. 

Within  the  crystallites  the  following  mechanisms  are  possible: 

1)  Direct  substitution  at  a  carbon  site  and  subsequent  diffusion  by 
any  of  the  carbon  set f*dif fusion  mechanisms: 

a)  Layer  plane  'acancy-atom  exchange 

b)  Direct  interchange  of  the  atom  with  a  carbon  atom 

c)  lnterla>er  migration 

d)  Vacancy -atom  exchange  between  layers 

2)  If  the  atoms  are  too  big  or  are  otherwise  not  chemically  suited 
to  substitute  for  a  carbon  atom,  then  diffusion  could  occur  by  the 
following: 

a)  Vacancy  cluster  —  foreign  atom  exchange  within  a  layer 

b)  Vacancy  cluster  —  foreign  atom  exchange  between  layers 

c)  Interlayer  migration 

d)  Migration  along  gross  stacking  errors,  such  as  a  missing 
portion  of  a  layer  plane 

Mechanisms  external  to  the  crystallites  include: 

1)  Migration  between  crystallites  along  the  top  layer  plane,  with 

trapping  at  vacancy  clusters  and  at  the  edge  of  the  top  layer  plane. 
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2)  'Migration  perpendicular  to  the  layers  along  the  crystallite  sides, 
proceeding  layer  by  layer. 

3)  Migration  perpendicular  to  the  layers  along  the  crystallite  sides, 
with  a  jump  distance  comparable  to  the  crystallite  thickness. 

4)  Migration  in  regions  of  gross  mismatch  where  crystallites  of 
differing  orientations  meet. 

5)  Migration  around  the  crystallites,  proceeding  from  trap  to  trap 
along  the  circumference  of  the  crystallites. 

The  mechanisms  listed  above  are  essentially  the  only  ones  possible  for 
metal  diffusion  in  PG.  Gross  pore  diffusion  is  not  a  possibility  because  the 
high  density  of  PG  rules  out  the  presence  of  the  large,  interconnected  pores 
that  occur  in  normal  synthetic  graphites.  As  will  be  discussed  shortly,  the 
small  pores  indicated  by  the  small  angle  scattering  studies  in  PG  are  included 
in  item  1  of  the  mechanisms  external  to  the  crystallites . 

In  the  following  paragraphs,  each  of  the  above  mechanisms  will  be  examined 
to  see  which  of  them  are  probable  mechanisms  for  metal  diffusion  in  PG. 

Considering  the  mechanisms  within  the  crystallites,  the  first  set  requires 
that  the  metal  atom  substitute  for  carbon  in  the  graphite  lattice.  To  substitute 
for  carbon,  the  metal  atom  must  have  the  proper  valence  to  participate  in  the 
graphite  structure  and  must  be  small  e  lough  not  to  seriously  distort  the  lattice. 
The  only  elements  small  enough  not  to  disturb  the  graphite  lattice  are  B,  H,  O, 
and  N.  All  of  these  combine  readily  with  carbon  to  form  compounds  but  none 


have  the  valence  neceaaary  to  participate  in  the  graphite  bonding  structure. 

Of  these  elements,  B  and  N  react  with  each  other  to  form  BN  which  has  struc* 
tures  very  similar  to  the  graphite  and  diamond  structures  (9). 

The  group  IVA  elements  have  the  proper  valence  to  substitute  for  carbon, 
but  even  the  Si  atom  is  50%  larger  than  the  carbon  atom,  which  is  a  serious 
size  difference  in  covalently -bonded  structures.  The  chemical  similarity  of 
Si  and  C  is  well  documented,  and  a-SiC  exists  in  a  variety  of  hexagonal  layer 
structures  (54).  There  is,  however,  no  evidence  for  direct  substitution  of  Si 
in  the  graphite  lattice  and  the  large  size  difference  would  probably  prevent 
direct  substitutior  from  occurring.  From  the  above  discussion,  it  can  be  seen 
that  the  mechanisms  of  metal  diffusion  requiring  direct  substitution  for  the 
carbon  atom  are  not  very  likely  and  the  other  mechanisms  must  be  considered. 

Of  the  second  set  of  mechanisms  for  diffusion  within  the  crystallites, 
those  mechanisms  based  on  the  motion  of  vacancy  clusters  are  not  very  likely, 
considering  the  high  energy  that  would  be  required  to  form  and  move  such  defects 
in  the  covalently -bonded  planes  (see  Chapter  11).  The  interlayer  mechanisms 
are  also  unlikely  except  for  the  very  smallest  atoms  and  hi  situations  in  which 
lamellar  compounds  form.  If  lamellar  compounds  arc  formed,  their  presence 
would  be  noted  bv  the  swelling  at  high  temperature  and  the  increase  in  the  inter¬ 
layer  spacing  that  accompanies  their  formation.  As  noted  in  Chapter  11,  the 
formation  energy  of  an  interlayer  Xe  atom  (4.4  A  diameter)  is  estimated  to  be 
350  kcal/mole  and,  although  the  transition  metals  are  somewhat  smaller  than 
the  Xe  atom,  electronic  interactions  and  the  fact  that  the  "free  space"  between 
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Ityert  decreafei  with  increaiing  temperature  (137)  would  prevent  appreciable 
interlayer  diffuaion.  Gross  stacking  errors  large  enough  to  permit  diffusion 
of  metal  atoms  in  PG  would  greatly  disturb  the  layer  plane  periodicity  and  this 
disturbance  would  be  reflected  in  the  crystallite  size.  Thus,  mechanisms  2d 
under  diffusion  within  the  crystallites  is  equivalent  to  mechanisms  1  under 
diffusion  external  to  the  crystallites.  From  this  discussion,  it  can  be  seen 
that  diffusion  in  PG  must  occur  external  to  the  crystallites  for  ail  but  the  smallest 
atoms . 

Metal  diffusion  in  PG  is  anisotropic.  This  discussion  will  now  be  separated 
into  mechanisms  for  diffusion  parallel  to  the  deposition  plane  and  mechanisms 
perpendicular  to  the  deposition  plane. 

Of  the  mechanisms  postulated  for  diffusion  external  to  the  crystallites, 
only  1,  4,  and  S  apply  to  diffusion  parallel  to  the  deposition  plane.  Mechanism 
1,  migration  between  crystallites  along  their  top  and  bottom  planes,  requires  a 
consideration  of  metal  atom  sizes  just  as  in  the  case  of  the  interlayer  mechan¬ 
ism  discussed  directly  above.  Metal  atoms  are  typically  3  A  in  diameter,  the 
"free  space"  between  carbon  atoms  within  a  crystallite  being  about  1.5  A  the 
2000^0  (137).  Therefore,  migration  between  crystallites  along  their  top  and 
bottom  planes  requires  that  the  crystallites  be  separated  sufficiently  to  provide 
room  for  the  diffusing  atoms.  The  voids  indicated  by  the  small  angle  scatter¬ 
ing  studies  will  provide  just  this  separation.  The  void  diameters  are  slightly 
less  than  the  crystallite  diameters,  and  the  thicknesses  are  approximately  50  A. 
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The  thickness  values  are  biased  to  the  larger  values  (95),  but  smaller  void 
thicknesses  are  adequate  to  allow  for  metal  migration  between  the  crystallites. 
In  mechanisms  1  (external  to  the  crystallites),  trapping  of  the  diffusing  atoms 
will  occur  at  the  unsatisfied  bonds  along  the  edges  of  the  top  and  bottom  layers 
of  the  crystallites  and  at  vacancy  clusters  in  the  tqi  and  bottom  planes . 

Because  there  was  no  evidence  of  preferential  boundary  penetration  for 
diffusion  parallel  to  the  deposition  plane,  mechanism  4  (external  to  the  crystal¬ 
lites)  does  not  apply.  Mechanism  S,  migration  around  the  circumference  of  the 
crystallites,  although  applicable  to  diffusion  parallel  to  the  deposition  plane, 
is  unimportant.  The  Jump  distance  for  mechanism  5  is  a  few  A  while  that  for 
mechanism  1  is  a  few  hundred  A.  Therefore,  the  diffusion  rate  based  on 
mechanism  5  would  be  inconsequential  compared  to  that  for  mechanism  1 . 

The  mechanisms  external  to  the  crystallite  labelled  2,  3,  snd  4  apply  to 
metal  diffusion  perpendicular  to  the  deposition  plane.  Migration  rates  perpen¬ 
dicular  to  the  layers  will  depend  upon  how  well  the  crystallites  meet  along  their 
common  boundaries  and  upon  the  nature  and  extent  of  contamination  by  other 
foreign  atoms  at  the  layer  plane  edges.  Furthermore,  the  unsatisfied  bonds 
at  each  layer  plane  edge  will  provide  many  trapping  sites,  and  the  average  jump 
distance  will  be  of  the  order  of  the  layer  plane  separation  distance.  This  situa¬ 
tion  eliminates  mechanism  3,  in  which  the  jump  distance  is  comparable  to  the 
crystallite  thickness. 

In  PG  there  are  regions  of  differing  misorientation;  within  a  minor  cone 
the  misorientation  between  crystallites  is  rather  small,  while  at  minor  cone 
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boundaries  it  is  somewhat  greater,  and  at  major  cone  boundaries  it  is  the 
greatest  of  all.  Because  the  major  cone  boundaries  contain  the  greatest  degree 
of  mismatch,  they  will  offer  the  least  resistance  of  all  the  boundaries  to  the 
diffusion  of  metals.  The  unsatisfied  bonds  at  each  layer  and  the  regions  of 
differing  misorientation  are  evidence  for  a  diffusion  mechanisms  perpendicu¬ 
lar  to  the  deposition  plane  that  includes  mechanisms  2  and  4. 

To  summarize: 

a  Metal  diffusion  in  PG  can  only  occur  external  to  the  crystallites 
a  Traps  exist  along  the  layer  plane  edges  and  at  vacancy  clusters 
in  the  top  and  bottom  layer  planes  of  the  crystallites 
a  Atoms  diffusing  on  the  top  layer  planes  of  the  crystallites  are  rela¬ 
tively  free  to  move,  provided  there  is  sufficient  room  between 
between  the  crystallites 

a  Primary  cone  boundaries  contain  regions  of  greater  misorienta¬ 
tion  than  do  secondary  cone  boundaries 
a  Diffusion  perpendicular  to  the  layer  planes  occurs  principally  in 
the  preimary  cone  boundaries,  the  regions  of  greatest  layer  plane 
misorientation 

a  Atoms  diffusing  perpendicular  to  the  layer  planes  encounter  frequent 
trapping  sites,  probably  at  each  layer  plane  edge 
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D.  Model  for  the  Diffusion  of  Metals  in  PG 


Two  separate  processes  for  metal  diffusion  in  PG  are  recognized;  one 
parallel  to  the  layer  planes  and  another  perpendicular.  The  high  degree  of 
preferred  orientation  of  the  layer  planes  in  PG  makes  it  possible  to  assume 
a  close  proportionality  between  the  diffusion  coefficient  measured  parallel  to 
the  deposition  plane  ^tid  the  diffusion  coefficient  parallel  to  the  layer 

planes.  The  nearly  straight  cone  boundaries  and  an  analysis  given  in  the  Appen¬ 
dix  make  it  possible  to  relate  the  diffusion  coefficient  measured  perpendicu¬ 
lar  to  the  deposition  plane  )  with  the  diffusion  coefficient  perpendicular 

to  the  layer  planes . 

The  expression  relating  the  diffusion  coefficient  for  non -substitutional 
diffusion  to  the  jump  distance  (d)  and  jump  frequency  (v)  will  be  given  here 
because  it  will  be  needed  later  in  the  development  of  the  model  for  metal  dif¬ 
fusion  in  PG  (136).  This  expression  is 

D  =  k  V 

where  k  is  a  numerical  factor  concerning  the  number  of  different  jumps 
possible  (136).  The  jump  frequency  is  simply  given  by  the  vibration  frequency, 

,  and  the  Boltzman  factor,  therefore 

D  =  k  d^  exp  (-AG*/RT) 

where  AG*  is  the  standard  Gibbs  free  energy  of  activation  per  mole  required 
to  surmount  the  energy  barrier  for  the  jump  process. 
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MeUl  diffusion  in  PG  parallel  to  the  layer  planes  takes  place  between  the 
crystallites  and  along  delaminations  in  the  material,  as  is  evidenced  by  auto¬ 
radiography  and  by  deductions  concerning  the  lack  of  diffusion  within  the  crys¬ 
tallites.  Because  the  top  and  bottom  surfaces  of  the  crystallites  are  relatively 
free  of  traps  and  the  edges  of  the  crystallites  have  the  highest  trap  density,  the 
important  Jump  distance  parallel  to  the  layer  planes  is  the  crystallite  diameter. 
One  can  determine  that  the  jump  distance  is  equal  to  the  crystallite  diameter  by 
considering  that  the  ratios  of  the  's  for  the  as -deposited  and  partially 

graphitized  PC's  are  approximately  equal  to  the  ratio  of  the  squares  of  the  crys¬ 
tallite  diameters  in  the  two  conditions.  The  basis  of  this  correlation  is  given 
in  the  above  equations  relating  the  diffusivity  to  the  square  of  the  Jump  distance. 
Vacancy-cluster  trapping  can  also  occur,  but  it  v'ill  not  be  important.  The 
energy  to  form  even  a  divacancy  is  estimated  to  be  150  kcal/mole  (sec  Chapter 
U),  and,  unless  an  active  source  of  vacancies  exists  the  few  available  clusters 
will  soon  saturate.  After  saturation,  vacancies  or  vacancy -clusters  will  not 
appreciably  affect  the  overall  migration  rate.  As  noted  above,  metal  diffusion 
across  the  surfaces  of  the  crystallites  requires  a  slight  void  region  between 
the  crystallites  to  allow  passage  of  the  diffusing  atoms;  the  small  angle  scatter¬ 
ing  from  PG  is  the  evidence  for  these  voids. 

The  microstructure  of  the  PG  is  not  delineated  by  diffusion  parallel  to  the 
layer  planes  because  there  are  about  as  many  trapping  sites  per  unit  area  near 
the  major  cone  boundaries  as  there  are  away  from  them.  Diffusion  parallel  to 
the  deposition  plane  in  the  major  cone  boundaries  should  occur  at  about  the  same 
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rate  as  cone  boundary  diffusion  perpendicular  to  the  deposition  plane .  The 
data  showed  that  the  diffusion  rate  measured  parallel  to  the  deposition  surface 
is  greater  than  that  measured  perpendicular  to  it  at  any  temperature;  there¬ 
fore,  diffusion  parallel  to  the  deposition  surface  is  dominated  by  diffusion 
across  the  crystallite  surfaces  and  not  b>  cone  boundary  diffusion. 

Diffusion  parallel  to  the  deposition  plane  occurs  on  approximately  flat 
planes  containing  the  crystallite  surfaces;  the  relatively  slower  diffusion  rate 
perpendicular  to  the  deposition  plane  prevents  interaction  between  the  parallel 
planes  of  diffusion.  If  this  interaction  were  appreciable,  it  would  complicate 
the  measurement  of  •  Diffusion  parallel  to  the  deposition  plane  is  slightly 

complicated  by  diffusion  around  the  circumference  of  the  crystallites  but,  as 
described  earlier,  mass  transport  by  this  circumferential  path  is  negligible 
compared  to  diffusion  across  the  cryt  tallite  surfaces.  The  delaminations  men¬ 
tioned  above  can  slightly  disturb  the  concentration  profile,  but  there  is  little 
difference  between  surface  diffusion  along  a  delamination  and  surface  diffusion 
across  a  crystallite  if  traps  exist  along  the  delamination  with  about  the  same 
spacing  as  the  crystallite  diameter.  The  extent  to  which  the  traps  in  the  delami¬ 
nations  do  not  approximate  the  traps  within  the  material  will  determine  the  effect 
of  the  delaminations  on  measuring  the  diffusion  coefficient  parallel  to  the  layer 
planes . 

In  apparent  contradiction  of  the  results  and  model  presented  above,  Wolfe, 
McKenzie  and  Borg  observed  that  extreme  graphitization  in  PC  (annealed  at 
3000°C  for  15  min)  decreased  the  diffusion  coefficient  of  Ni-63  parallel  to  the 
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layer  planes  by  about  three  orders  of  magnitude.  PC  this  well  graphitized 
would  be  characterized  by  a  very  large  crystallite  size,  an  extremely  sharp 
layer  plane  texture  (5),  and.  probably,  a  complete  lack  of  small  angle  scatter¬ 
ing.  The  lack  of  misorientation  between  the  crystallites  would  permit  consider¬ 
able  inter-cry  SI  Hite  bonding  over  adjacent  layer  plane  faces,  and  this  bonding 
would  tend  to  eliminate  the  voids  between  the  crystallites.  In  turn,  elimination 
of  the  voids  would  reduce  the  jump  distance  and/or  increase  the  activation  energy, 
and  both  of  these  changes  would  decrease  the  observed  diffusion  coefficient.  The 
heat-treated  PC's  used  in  this  study  were  not  nearly  as  well  graphitized  as  heat- 
treated  PG  studied  by  WMB . 

The  activation  energy  for  diffusion  across  the  crystallite  surfaces  consists 
of  two  terms:  that  for  "untrapping"  and  that  for  motion  across  the  crystallite 
surface.  The  untrapping  energy  should  be  called  the  energy  to  form  a  crystal¬ 
lite  surface  atom  **  should  be  related  to  the  bond  energy  of 

the  metal  for  carbon  and  the  size  of  the  metal  atom.  Conversely,  the  energy 
for  motion  across  the  crystallite  will  primarily  depend  upon  the 

crystallite  separation  and  the  metal  atom  size.  From  all  the  above  discussion, 
the  diffusion  coefficient  parallel  to  the  deposition  plane  can  be  expressed  as 

°<a>  '  ‘‘a  ^a  "o'-P  <-'^®*f<a>  ' 

As  noted  above,  k„  is  a  numerical  factor  related  to  the  number  of  different 

A 

jumps  possible,  is  the  crystallite  diameter  and  jump  distance  for  the  dif¬ 
fusion,  and  is  the  vibration  frequency. 
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The  mechanism  for  diffusion  perpendicular  to  the  layer  planes  will  now  be 
discussed.  Diffusion  in  this  direction  must  occur  down  the  crystallite  edges  for 
the  reasons  given  above .  Between  the  crystallites  within  the  minor  cones  and 
along  the  minor  cone  boundaries,  the  relative  misorientation  is  slight,  which 
probably  results  in  a  low  degree  of  disorder.  Large  lattice  strains  would  be 
required  for  metal  diffusion  along  these  boundaries,  making  the  activation 
energy  for  diffusion  much  larger  here  than  at  the  major  cone  boundaries  where 
the  misorientation  is  greatest.  Because  the  pre-exponential  factor  for  all  these 
paths  will  be  about  the  same,  the  diffusion  rate  at  a  given  temperature  will  be 
greatest  for  the  path  of  least  activation  energy,  which  is  the  primary  cone 
boundary.  Primary  cone  boundary  diffusion  is  demonstrated  in  the  autoradio¬ 
graphs  of  diffusion  perpendicular  to  the  layer  planes  and  in  the  observation  that 
partial  graphitization  decreases  the  diffusion  coefficient  measured  in  that 
direction.  Further  evidence  for  this  portion  of  the  model  is  given  by  the  effect 
of  graphitization  on  .  Partial  graphitization  results  in  a  decrease  in  the 

misorientation  angle  between  crystallites  and  a  decrease  in  .  These 

observations  are  consistent  with  the  arguments  just  given  concerning  diffusion 
rates  in  regions  of  differing  misorientations . 

Diffusion  perpendicular  to  the  layer  planes  is  complicated  by  a  high  lateral 
diffusion  rate  parallel  to  them;  however,  ir  ':an  be  shown  (see  Appendix)  that 


<c> 
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where  I  is  the  major  cone  boundary  thickness,  the  boundary  diffusion 
coefficient,  and  the  value  obtained  by  treating  the  observed  profile  as 

being  the  result  of  diffusion  into  a  semi-infinite  solid  under  constant  surface 
composition  conditions .  This  proportionality  between  Dj^  and 
important:  it  shows  that  the  activation  energy  computed  from  the  D^j.^data 
is  the  same  as  that  for  Dj^  .  It  also  states  that,  although  the  scatter  in  the 
diffusivity  data  is  dominated  by  errors  in  temperature,  some  contribution  to 
this  scatter  can  come  from  having  experimental  co..ditions  that  do  not  meet  the 
assumptions  used  in  the  derivation  given  in  the  Appendix.  These  assumptions 
are  that  the  observed  profile  can  be  represented  by 


and  that  constant  surface  composition  conditions  existed  throughout  the  diffusion 
anneal.  If  the  tracer  reservoir  on  the  surface  is  depleted  during  the  experiment, 
one  of  the  boundary  conditions  used  in  the  derivation  will  be  violated. 

It  is  reasonable  to  assume  that  the  jump  distance  for  diffusion  perpendicular 
to  the  layer  planes  is  equal  to  the  layer  plane  separation.  The  edges  of  each 
layer  plane  have  a  large  number  of  unsatisfied  bonds;  traps  can  thus  be  expected 
to  exist  at  each  layer  plane.  Such  trapping  yields  a  jump  distance  of  3.4  A, 
which  is  virtually  independent  of  the  crystallite  thickness.  If  the  jump  distance 
is  independent  of  crystallite  thickness,  the  activation  energy  for  diffusion  per¬ 
pendicular  to  the  layer  planes  would  apparently  have  to  increase  in  the  graphi- 
tized  material  to  account  for  the  observed  decrease  in  D^^^^with  graphitization. 
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There  was  no  significant  increase  in  activation  energy;  therefore,  the  cause 
of  the  decrease  in  have  to  be  found  elsewhere.  In  the  Appendix,  it 

is  shown  that 


It  is  reasonable  to  expect  that  the  width  of  the  boundary  containing  the  less 
organized  carbon  would  decrease  slightly  upon  graphitization  due  to  the  de¬ 
crease  in  misorientation.  A  one-half  decrease  in  i  would  account  for  the  slight 
decrease  in  for  the  partially  graphitized  PG.  A  gradual  increase  in 

the  activation  energy  for  diffusion  perpendicular  to  the  layer  planes  with  increas¬ 
ing  graphitization  is  to  be  expected,  however,  because  the  major  cone  boundaries 
will  become  regions  of  less  disorder  as  the  misorientation  decreases  and  these 
boundaries  will  approach  the  state  of  the  minor  boundaries.  This  state  was  not 
reached  in  the  graphitizing  heat -treatment  used  in  this  study. 

From  the  above  discussion,  the  observed  diffusion  coefficient  perpendicular 
to  the  layer  planes  can  be  written 

=  k  d^..„  i'  exp  (AG*  /RT  «  ID. 

where  i  is  the  boundary  thickness,  a  numerical  factor  related  to  the  num¬ 
ber  of  different  jump  directions,  dQQQ2  the  interlayer  spacing,  the 

standard  Gibbs  free  energy  for  motion  perpendicular  to  the  layer  planes,  and 
the  vibration  frequency . 
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The  ratio  of  the  pre -exponential  factors  for  diffusion  parallel  and  perpen¬ 
dicular  to  the  layer  planes  should  be  equal  to  the  ratio  of  the  numerical  factors 
times  the  square  of  the  inferred  Jump  distances: 


o-a 


-A 


*^o-c  K  **0002 


Assuming  a  ,  and  taking  data  for  as -deposited  PG  from  Table  IS, 


o-a  _  3.2  >^10 


o-c  1.2x10 


r2 


2.7  xlO 


0002 


Similarly,  for  heat-treated  PG 


mj 


0.  54  X  10 


o-a 


1.8  X  10' 


o-c  6. 5  X  10 


-3 


2.5  xlO' 


a 


0002 


(e) 


0.31  X  10 


For  the  as -deposited  PG's,  the  ratios  differ  by  a  factor  of  5  and  for  the  heat- 
treated  PG's  the  ratios  differ  by  a  factor  of  8  .  If  the  values  of  ^  for  the 
as-deposited  and  heat-treated  PG's  were  normalized  to  the  same  primary  wonc 
boundary  thickness,  the  agreement  would  be  improved.  In  the  model  discussed 
above,  it  was  postulated  that  the  primary  cone  boundary  thickness  was  decreased 


150 


* 


I 


by  one-half  during  the  graphitizing  heat-treatment;  therefore,  normalization 
of  the  Dq_c  terms  would  require  multiplying  (HT)  by  2.  With  this 

factor  of  2  ,  the  ratios  of  the  D^'s  and  the  ratios  of  the  squares  of  the  Jump 
distances  would  differ  from  each  other  by  a  factor  of  5  .  The  agreement  could 
be  improved  further  if  were  about  5  .  'Fhese  numerical  factors  arc 

proportional  to  the  reciprocal  of  the  number  of  different  Jump  directions  from 
a  particular  site  (136),  and  estimates  of  these  factors  arc  given  below. 

For  diffusion  perpendicular  to  the  deposition  plane,  an  atom  diffusing  in 
the  primary  cone  boundary  can  Jump  in  four  different  directions:  to  the  next 
layer  plane  above  its  position,  to  the  plane  below  its  position,  to  the  left  of  its 
position,  or  to  the  right  of  its  position.  All  of  those  jumps  arc  of  equivalent 
size.  Conversely,  for  diffusion  parallel  to  the  layer  planes,  the  number  of  dif¬ 
ferent,  equivalent  Jumps  is  restricted  to  two:  forward  or  backward  one  crystal 
lite  diameter.  Jumps  to  either  side  arc  possible,  but  their  disunce  is  much 
smaller  than  the  crystallite  diameter  Jump  disunce,  and  therefore,  they  are 
not  of  consequence  in  the  diffusion  process.  These  arguments  yield  k.  «  2  k. 
Another  factor  of  two  might  be  obtained  from  arguments  concerning  the  vibra¬ 
tion  frequencies;  however,  the  D^'s  are  not  known  with  sufficient  accuracy  to 
warrant  such  refinements.  The  standard  deviation  of  is  about  6%,  a  slight 
increase  in  Q  for  Nb  would  also  reduce  the  pre -exponential  factor  ratio  and 
tend  to  improve  the  agreement. 
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B.  Discussion  of  Grain  Boundary  Diffusion  Theories 

Grain  boundary  diffusion  in  metals  has  been  studied  rather  extensively 
(136,138  through  141)  and  LcClairc  (142)  has  recently  revic\^'ed  the  subject. 

All  the  grain  boundary  theories  presented  to  date  have  been  developed  for  con¬ 
current  boundary  and  volume  diffusion  in  an  otherwise  isotropic  solid  and  all 
predict  that  the  logarithm  of  the  concentration  should  decrease  approximately 
as  the  first  power  of  the  distance.  When  grain  boundary  diffusion  is  important 
this  dependency  is  observed  for  metals  diffusing  in  polycrystalline  metals  (143). 
In  contrast  with  these  results,  the  data  of  this  investigation  and  those  of  Wolfe, 
McKenzie  and  Borg  definitely  show  that  this  dependency  is  not  observed  for 
metals  diffusing  in  PG . 

The  reason  for  this  difference  between  grain  boundary  diffusion  in  metals 
and  in  PG  lies  in  the  lack  of  volume  diffusion  and  the  peculiar,  anisotropic  defect 
structure  of  PG .  Diffusion  in  PG  parallel  to  the  layer  planes  is  equivalent  to 
diffusion  on  many  parallel  planes,  each  with  a  line  source.  With  modifications 
of  the  concentration  terms,  this  situation  is  equivalent  to  linear  diffusion  from 
a  plane  into  a  solid  (144).  The  same  could  be  said  for  diffusion  within  the  pri¬ 
mary  cone  boundaries  perpendicular  to  the  layer  planes  except  for  the  lateral 
diffusion  parallel  to  the  layer  planes .  These  sinks  extending  away  from  the  cone 
boundaries  are  supplied  at  the  rate  of  diffusion  down  the  boundary  and,  because 
only  lateral  and  cone  boundary  diffusion  can  occur,  the  concentration  profile 
measured  by  sectioning  techniques  is  simply  proportional  to  the  boundary  con¬ 
centration  profile  (see  Fig.  27).  Therefore,  the  measured  diffusion  coefficient 


is  approximately  proportional  to  the  boundary  diffusion  coefficient.  The 
mathematics  of  this  discussion  arc  given  in  the  Appendix. 

F.  Discussion  of  the  Model  and  Comparison  with  the  Data  of  Wolfe, 

McKenzie  and  Borg 

The  data  of  Wolfe,  McKenzie  and  Borg  arc  in  general  agreement  with  the 
data  of  this  investigation.  While  WMB  studied  the  diffusion  of  several  different 
metals  in  PG,  the  present  study  concentrated  on  the  diffusion  of  one  metal  in 
several  different  PC's.  A  comparison  of  all  the  data  shows  that  the  diffusion 
of  metals  in  PG  is  dependent  upon  the  kind  of  metal,  the  direction  of  diffusion, 
and  the  crystallite  diameter,  but  not  upon  the  microstructure.  The  model  pro¬ 
posed  by  WMB  was  based  primarily  upon  entropy  arguments  concerning  the  large 
values  of  D^_^  .  These  arguments  led  them  to  propose  that  the  Jump  distance 
was  approximately  equal  to  the  crystallite  dimension  in  the  direction  of  diffusion. 
Their  model  is  similar  to  that  proposed  here,  except  for  the  Jump  distance  per¬ 
pendicular  to  the  layer  planes.  The  WMB  model  assumes  this  Jump  distance  to 
bo  the  crystallite  thickness,  whereas  in  the  model  proposed  in  the  present  in¬ 
vestigation,  the  Jump  distance  {perpendicular  to  the  layer  planes  is  the  inter¬ 
layer  spacing. 

Some  of  the  results  for  WMB's  work  and  this  investigation  arc  summarized 
in  Table  16,  together  with  other  information.  First,  it  may  be  seen  that 

4 

Do-.  >  Do-c  in  all  cases,  with  the  ratio  being  about  10  for  Nb  diffusing 
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Table  16 

COMPARISON  OF  DATA  FROM  THIS  INVESTIGATION  WITH  THOSE  OF  WOLFE,  MCKENZIE  AND  BORG 
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(b)  Reference  10*» 


in  as -deposited  PC  and  somewhat  less  for  WMB’s  data  on  Th,  U,  and  Ni .  The 
model  proposed  in  this  study  requires  that,  to  a  first  approximation,  this  ratio 
must  be  constant  for  a  particular  PC,  regardless  of  the  diffusing  metal.  This 
ratio  is  not  constant  for  WMB's  data;  however,  it  should  be  noted  that  is 
very  sensitive  to  the  choice  of  Q  ,  and  it  is  possible  for  the  ratio  of 

to  vary  by  a  factor  of  10  through  the  choice  of  and  .  This  ratio 
is  as  low  as  20  for  the  case  of  U  diffusing  in  PG,  as  determined  by  WMB.  Such 
a  low  value  cannot  be  explained  by  a  smaller  cyrstallite  diameter  because  WMB 
estimated  that  the  crystallite  diameter  of  the  PG  used  in  their  investigation  was 
about  150  A.  This  diameter  is  somewhat  smaller  than  the  values  obtained  for 
the  PC's  used  in  this  study,  but  the  difference  would  only  account  for  a  factor  of 
2  in  the  ratio  of  the  pre -exponential  factors. 

The  relative  values  of  the  activation  energies  parallel  and  pemcndicular  to 
the  layer  planes  also  differ  for  the  different  metals  studied  in  this  investigation 
and  that  of  WMB.  As  shown  in  Table  16,  for  Nb  and  Th,  and  for 

U  and  Ni,  Q  <  Q  .  In  the  model  proposed  here. 


Q  =  AH,  +  AH 
^a  f<a>  m<a> 


and 


Q  »  AH 
c  m<c> 


To  a  first  approximation,  one  might  expect  “ 

represent  untrapping.  This  approximation  would  make 


H 

m<c> 

Qn  >  Qr 


because  both 
for  all  metals. 
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which  is  not  the  ease  for  U  and  N».  In  hgiit  of  this,  the  approximation  concern¬ 
ing  the  activation  energies  does  not  appear  to  be  valid. 

It  is  also  possible  that  AH,  ^  and  AH  ^  are  related  to  the  metal - 
*  f<a>  m<c> 

carbon  lx)nd  strength.  Table  16  lists  the  heats  of  formation  for  the  carbides  of 
Nb,  Til,  and  U  at  298°K.  As  may  be  seen,  the  activation  energies  for  diffusion 
are  two  to  three  times  the  heat  of  formation  of  the  carbides,  indicating  that 
either  the  metal  atoms  are  more  tightly  bound  in  the  PC  structure  than  in  the 
carbide  or  that  the  contribution  of  strain  energy  to  the  activation  energies  for 
diffusion  is  large.  Both  effects  probably  contribute  to  the  difference.  The  high 
value  of  for  Tli  in  PG  is  very  likely  related  to  its  large  size.  In  the  model 
developed  by  this  investigation,  diffusion  parallel  to  the  deposition  plane  occurs 
by  atoms  moving  between  the  crystallites.  Thus,  the  much  larger  Tli  atom 
should  have  a  higher  motion  energy  for  that  process  than  do  the  smaller  atoms. 

As  listed  in  Table  16,  the  activation  energies  for  diffusion  of  Ni  and  Ag  in 
PG  are  much  lower  than  those  of  U,  Th,  and  Nb.  Ni  and  Ag  do  not  form  carbides 
and  their  activation  energies  for  diffusion  are  aliout  equal  to  those  for  the  carbide 
formers  after  the  absolute  value  of  the  heat  of  formation  of  the  carbide  has  been 
subtracted.  This  observation  is  in  agreement  with  the  model,  in  that  the  activa¬ 
tion  energy  for  diffusion  of  a  carbide  former  in  PG  should  lx*  larger  than  that  of 
a  non -carbide  former  by  the  heat  of  formation  of  the  carbide  Ixcause  the 
"untrapping"  energy'  for  a  carbide  former  should  include  the  heat  of  formation 
of  the  carbide.  An  apparent  exception  to  this  argument  is  the  value  of  99  kcal/ 
mole  for  diffusion  of  Ra  parallel  to  the  deposition  plane  in  PG.  This  value  is 
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larger  than  that  observed  for  Nb  but  the  atom  sizes  differ  appreciably.  The 
atom  size  of  Ra  is  not  well  known;  however,  the  atom  radius  of  Ra  may  bo 
estimated  from  the  atom  sizes  of  Kr  and  Xe,  using  the  values  given  in  Cullity(106) 
as  2.5  A.  An  atom  size  this  large  would  result  in  a  high  activation  energy  for 
diffusion  parallel  to  the  deposition  plane  foi  the  same  reasons  given  above  for 
the  Th  diffusion. 

G.  Relevance  of  Model  to  Diffusion  in  Other  Graphites  and  Pyrolytic  Graphites 

The  model  for  diffusion  in  PG  developed  in  this  investigation  should  have 
application  to  metal  diffusion  in  other  forms  of  graphite  and  in  the  recently 
improved  forms  of  PG.  Some  of  the  model's  implications  for  diffusion  in  the 
other  graphites  will  be  given  here. 

As  discussed  in  Chapter  III,  conventional  synthetic  graphites  consist  of 
well  graphitized  regions  surrounded  by  somewhat  less  graphitized  regions  that 
are  derived  from  the  carbonization  and  graphitization  of  the  petroleum -coke 
binder.  Conventional  graphites  contain  a  spectrum  of  crystallite  sizes  from 
over  1000  A  to  less  than  50  A,  with  the  well  graphitized  regions  having  the 
larger  crystallites  and  the  binder  regions  the  smaller.  The  average  diame¬ 
ter  of  crystallites  is  about  500  A  (27).  Conventional  graphites  also  possess  a 
spectrum  of  pore  sizes,  some  of  which  arc  accessible  to  high  pressure  Hg 
porosimetry.  Because  the  easily  sheared  layer  planes  tend  to  lie  parallel  to 
the  flow  lines,  some  preferred  orientation  exists  if  the  graphite  body  has  been 
subjected  to  hot  working. 
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Because  the  activation  energies  parallel  and  perpendicular  to  the  layer 
planes  are  comparable  and  the  pre -exponential  factor  for  diffusion  paral¬ 
lel  to  the  layer  planes  is  large,  the  diffusion  rate  parallel  to  the  layer  planes 
will  dominate  in  conventional  graphite.  This  effect  would  be  accentuated  in 
conventional  grapliite  because  of  large  average  crystallite  size.  Another  com¬ 
plication  would  arise  due  to  the  wide  spectrum  of  crystallite  sizes  which  would 
result  in  an  even  greater  spread  in  the  pre-exponential  factor.  Furthermore, 
there  would  definitely  be  an  anisotropy  of  diffusion  related  to  the  layer  plane 
texture  of  the  graphite.  However,  the  ratio  of  the  diffusivities  perpendicular 
and  parallel  to  the  texture  would  be  considerably  less  than  in  PG  because  con¬ 
ventional  graphites  do  not  have  the  high  degree  of  layer  plane  texture  possessed 
by  PG. 

If  diffusion  parallel  to  the  layer  planes  is  the  dominate  mechanism  for  metal 
diffusion  in  conventional  graphites,  then  one  might  expect  to  observe  the  same 
activation  energy  as  that  for  diffusion  in  PG  parallel  to  the  deposition  plane. 
Tliere  are,  however,  additional  complications .  It  is  not  likely  that  conventional 
graphites  would  possess  the  uniform  layer  plane  voids  found  in  PG  because  of  the 
differences  in  manufacture.  When  the  volatile  elements  are  driven  off  during 
the  manufacture  of  conventional  graphite,  the  carbon -foreign  element  bonds 
are  broken  and  new  carbon -carbon  bonds  are  formed.  Many  of  these  form  what 
are  called  cross-links  (57)  which  are  bonds  that  do  not  form  part  of  the  layer 
plane  network  but  extend  away  from  it  in  different  directions.  These  bonds 
interfere  with  the  graphitization  process  and  also  would  interfere  with  the  motion 
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of  foreign  atoms  diffusing  over  the  surfaces  of  the  crystallites.  In  addition, 
the  purity  of  conventional  graphites  is  generally  less  than  that  of  PC,  with  the 
foreign  atoms  tending  to  congregate  at  the  edges  of  the  crystallites  for  the 
reasons  discussed  earlier.  Because  these  foreign  atoms  occupy  the  same 
p’aces  that  the  diffusing  atoms  should  occupy,  some  interference  is  inev.  able. 
For  these  reasons,  tlie  activation  energy  for  metal  diffusion  in  conventional 
graphites  should  be  somewhat  greater  than  that  for  diffusion  parallel  to  the 
deposition  surface  in  PG. 

Examination  of  WMB' s  data  (T;.!jle  3)  indicates  that  the  diffusion  of  L  -232 
in  both  National  Carbon  ZTA  graphite  (high  density)  and  Speers  S700  graphite 
(low  density)  follows  the  above  description.  The  p  re -exponential  factors  are 
two  orders  of  magnitude  greater  than  the  pre -exponential  factor  parallel  to 
the  PG  deposition  plane  and  the  activation  energies  are  slightly  higher  than  the 
activation  energy  parallel  to  the  deposition  plane  in  PG.  The  anisotropy  of  Th 
diffusion  in  ZTA  also  fits  the  above  discussion. 

None  of  the  other  data  for  metal  diffusion  in  graphites  listed  in  Chapter  II 
shows  activation  energies  as  high  as  those  found  in  this  investigation  or  by  WMB 
for  the  same  metals  studied.  It  should  be  noted,  however,  that  these  data 
generally  represent  very  questionable  results  derived  from  two  or  three  points 
on  the  concentration  profiles  with  data  at  three  or  so  temperatures  used  to  deter¬ 
mine  activation  energies.  As  discussed  earlier,  tiie  studies  of  the  migration  of 
fission  fragments  are  additionally  complicated  by  radiation  damage  and  the  pres¬ 
ence  of  other  elements;  however,  the  generally  lower  activation  energies  are 
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probably  indicative  of  mass  transport  within  or  along  the  pores  rather  than 
between  the  individual  crystallites,  as  in  PG.  All  of  the  elements  listed  in 
Table  2  have  rather  low  melting  temperatures  and  high  vapor  pressures, 
whic  1  would  tend  to  enhance  the  pore  migration  mechanism. 

Recent  developments  in  PG  technology  have  led  to  the  production  of 
thicker  deposits  with  a  greater  uniformity  of  structure  and  less  tendency 
toward  delamination  than  were  previously  available.  The  differences  between 
these  newer  PC's  and  the  ones  used  in  this  study  are  differences  only  in  quality 
and  uniformity .  Because  it  was  shown  that  metal  diffusion  in  PG  does  not  depend 
upon  microstrucure  but  rather  upon  crystallite  size,  metal  diffusion  coefficients 
in  the  new  FG's  would  not  differ  appreciably  from  those  found  in  this  study  or  that 
of  WMB.  Any  variation  should  be  predictable  on  the  basis  of  crystallite  varia¬ 
tions,  as  described  above. 

H.  Commems  Concerning  Diffusion-Alloying  in  PG 

One  of  the  purposes  of  this  investigation  was  to  consider  the  possibility  of 
diffusion-alloying  PG  with  metals.  From  the  diffusion  data  obtained  in  this 
study  and  that  of  WMB,  it  is  evident  that  the  diffusion  coefficients  of  many  metals 

in  PG  are  high  enough  to  allow  appreciable  pcietration  within  reasonable  times. 

-7  1 

However,  the  concentrations  involved  are  very  small  (~10  gm/cm  )  and  the 
effect  of  the  alloy  element  might  be  small.  It  was  shown  that  the  diffusing  metal 
atoms  are  restricted  to  the  crystallite  surface  and,  since  the  mechanical  (5)  and 
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oxidation  (15)  behavior  of  PG  is  largely  controlled  by  the  interaction  of  the 
crystallites  and  their  surfaces,  these  few  atoms  concentrated  on  the  crystal¬ 
lite  surfaces  could  appreciably  alter  some  of  the  properties  of  PG.  The  si7.e 
of  the  effect  would  have  to  be  determined  by  experiment. 
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Chapter  X 


SUMMARY 

The  purposes  of  this  investigation  were  to  determine  the  fundamental 
mechanisms  of  metal  diffusion  in  PC,  to  increase  the  understanding  of  the 
structure  of  PC  through  these  studies,  and  to  examine  tlie  possihlity  of  diffusion¬ 
alloying  PC.  A  literature  survey  revealed  that,  with  the  exception  of  recent  data 
of  Wolfe,  McKenzie,  and  Borg  (17)  and  possibly  that  of  Hennig  (18),  there  were 
virtually  no  reliable  datr.  for  the  diffusion  of  metals  in  graphites.  This  lack  of 
data  was  one  of  tlie  underlying  reasons  for  tlK*  investigation. 

In  the  initial  studies,  diffusion  reactions  of  Re  with  PC  w’ere  studied  and 
a  metastable  carbide  of  Re  was  discovered.  Tliis  carbide  could  not  be  pro¬ 
duced  in  quantity,  but  information  regarding  its  phase  reactions  with  Re  metal 
and  its  cr>'stal  structure  was  obtaincxl.  Tlic  crystal  structure  of  tlx;  carbide 
was  found  to  Ik-  simple  cubic  with  a  unit  cell  edge  of  2.8-f  A.  However,  because 
the  atom  positions  within  tlK-  unit  cell  could  not  be  determined,  the  stoichiometry 
of  the  carbide  is  not  known.  Tlie  largest  quantities  of  the  carbide  were  pro¬ 
duced  in  diffusion  couples  by  a  eutectic  reaction  with  Re.  Tlie  eutectic  tempera¬ 
ture  betw’een  Re  and  Re  carbide  was  estimated  to  be  ~20(X)°C  (±200*^0),  com¬ 
pared  with  2480®C  for  tlw  eutectic  temperature  beuvt-en  Re  and  C  (103). 

Tlie  techniques  for  studying  metal  diffusion  in  PC  were  developed  with  a 
radioisotope  of  Re,  but  hlglier  specific  activity  was  required  than  was  readily 
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availal)I(,*  wiili  Kc  and  a  carrier-free  isotope  of  Nh  was  chosen  to  continue  the 
study  of  metal  diffusion  in  PG. 

'Hie  diffusion  coefficients  of  Nb  in  PG  were  measuretl  l)etween  14  50  ard 
1950*^C  in  iliree  different  kinds  of  PG.  llie  PG's  were  cliosen  to  l)e  representa¬ 
tive  of  sulisirate  nucleated,  combination  substrate  and  re^eneralively  nucleated, 
and  regenera  lively  nucleated  microstruciures.  The  diffusion  coefficients  were 
determined  for  two  conditions  of  each  PG  (as-deposited  and  partially  graphili/ed) 
aiul  111  llie  two  different  "crystallographic"  directions  of  PG  (parallel  and  perpen¬ 
dicular  to  the  deposition  surface).  Autoradiography  was  employed  to  ohsi*rve  tlie 
distrilHition  of  the  diffusing  atoms. 

Tile  structure  of  each  PG  was  characterized  by  the  application  of  many  tech- 
ni(|ues  of  metallography,  x-ray  diffraction,  and  density  determinations .  Wiili 
these  tccliniques,  llie  structure  was  de.scrilied  in  terms  of; 

•  Metallographic  structure  at  high  and  low  magnification 

•  Structure  observed  by  transmission  electron  microscopy 

•  .Stacking  error  probability  or  degree  of  graphitization,  p 

•  Crv.stalliie  thickness  anil  iliameler,  L  and  L 

c  a 

•  IX'gree  of  preferred  orientation 

•  Orientation  dependence  of  structural  parameters,  p  anil  L, 

u 

•  Submicroscopic  void  size,  determined  from  small  angle  scattering 
of  x-ravs 

•  Oisplaceil  volume  density  anil  x-ray  density 
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A  mathematical  analysis  of  diffusion  in  layered  structures  was  made,  and  corre¬ 
lations  lx;tween  the  diffusion  and  structural  parameters  were  noted. 

Tile  diffusion  coefficients  for  Nb  in  PC  in  a  particular  direction  could  lx.* 
represented  by  a  single  activation  energy  ami  variable  pre-exponential  factors. 
Tlie  activation  energy  was  independent  of  the  PG  microstructures  ami  states  of 
graphitization  used  in  this  study.  Tlie  activation  energ\’  for  Nb  diffusion  in  PG 
parallel  to  tiK*  de{X)sitlon  surface  was  96.2  kcal/mole  while  that  for  Nh  diffu¬ 
sion  perpendicular  to  the  deposition  surface  was  73.1  kcal/mole 

Tlie  pre -exponential  factors  (D^)  for  **  given  direction  in  the  as -deposited 
PG  did  not  differ  appreciably,  regardless  of  microstructure.  Tie  re  was,  how¬ 
ever,  a  significant  change  in  the  Dq’s  brought  alxnit  by  the  graphiiizing  heat- 

treatment.  Parallel  to  the  deposition  plane  the  pre-ex|>onential  factors 

2  2 

(Do.g  )  were  all  approximately  3.2  x  10  cm  /sec  for  the  as-deposited  PG’s 
3  2 

and  1 .8  10  cm  /sec  for  the  heat-treated  materials.  For  diffusion  peipen- 

dicular  to  the  deposition  plane,  these  factors  (l)„  „  )  were  approximately 
1 .2  ^  10*^  cm^/sec  for  the  as-deposiietl  PG's  and  6.5  x  10  cm^/sec  for  the 
iK'at-treated  PG's.  There  were  slight  differences  in  the  pre-exponential  factors 
in  each  direction  among  the  various  IXj’s  Ihji  the  differences  were  of  the  oriler 
of  the  uncertainty  in  the  D^’s 

Autoradiographs  taken  perpemlicular  to  the  mass  flux  for  diffusion  per¬ 
pendicular  to  the  dejxisition  plane  demonstrated  that  diffusion  peipendicular  to 
the  deposition  plane  occurs  principally  along  the  primary  cone  Ixiundaries  which 
arc  the  regions  of  greatest  layer  plare  misorieiitation .  Autoradiographs  taken 


» 
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pcrpcncllcuhir  to  the  deposition  plane  for  diffusion  parallel  to  the  deposition 
plane  showed  a  uniform  distrilnition  of  radioactivity,  except  for  occasional 
streaks  caused  hy  delamination. 

A  consideration  of  the  possible  mechanisms  for  diffusion  within  tlie  crys¬ 
tallites  showed  that  such  mechanisms  are  very  unlikely  except  for  the  smallest 
atoms:  U,  H,  O,  andN.  I'rom  valence'considerations,  none  of  these  will 
substitute  for  carlxjn  in  the  layer  planes,  Ixit  an  interlayer  mij;ration  mechan¬ 
ism  is  possible.  All  other  atoms  are  too  larp.*  to  diffuse  within  tlie  crystallites 
except  in  situations  in  which  lamellar  compounds  form.  Hvidence  was  given 
that  lamellar  compounds  did  not  form  during  this  investigation. 

Considerations  of  tho  possible  mechanisms  external  to  the  cr>'stallites  led 
to  the  following  conclusions: 

•  Traps  exist  along  the  layer  plane  edges  and  at  vacancy  clusters  in 
the  top  and  Ixittom  layer  planes  of  the  crystallites 

•  Primary  cone  Ixiundaries  contain  regions  of  greater  misorientation 
than  ilo  secondary  cone  Ixiundaries 

•  Atoms  diffusing  on  the  top  layer  plane  of  a  crystallite  are  relative!) 
free  to  move,  provided  there  is  sufficient  rt>om  lK*tween  the 
crystallites 

•  Diffusion  per|KMuhcular  to  the  layer  planes  occurs  principally  in  the 
primary  cone  Ixiundaries,  along  the  paths  of  greatest  layer  plane 
misorientation 


165 


•  Atoms  diffusing  perpendicular  to  ilie  layer  planes  encounter  frequent 
trapping  sites,  prolwhly  at  each  layer  plane  edge. 

Tl>e  alxjve  conclusions  and  an  analysis  of  the  ratios  of  the  pre-exponential 
factors  for  the  different  directions  and  conditions  revealed  tliat  the  luniinmrmnl 

jump  dtstancr  parallel  to  the  deposition  plane  is  the  crystallite  diarm  ter  u  hile  the 
jump  distance  perpendicular  to  the  layer  planes  is  the  layer  plane  st  paranon  distance. 

The  following  equations  resulted  when  these  conclusions  were  comhinetl 
W’ith  elementary  interstitial  diffusion  theory: 


D 

D 


<a> 

<c> 


k  exp  (-AG*  ) 
n  A  Q 

‘‘c  '*M02 


where  is  the  crystallite  diameter,  interlayer  spacing, 

G*  ^  and  G*  the  standard  free  energies  of  activation  for  the  pro- 
<a>  <c>  ^  ' 

cesses  parallel  and  perpendicular  to  the  deposition  plane,  respectively,  and 

k  and  k  numerical  factors  related  to  the  numlxrr  of  different  jump  directions 

for  a  given  position.  A  consideration  of  the  different  possible  jumps  led  to  the 

relation  *  2  k  . 

a  c 

An  analysis  of  metal  diffusion  down  the  primary  cone  ixiundaries  of  PC  with 
lateral  diffusion  occurring  only  parallel  to  the  deposition  surface  yielded 


D  ^  a  ID, 
<c>  b 
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f 


where*  /  IS  the  primary  cone  tx)un(lar>  thickness,  1)^  the  Ixjuiulary  diffusion 
coefficient,  and  I)  llie  measured  diffusion  coefficient .  Tlie  analysis  re - 
quired  lliat  1)^^, ^I>e  obtained  l)y  treating  the  composition  profiles  as  if  they 
were  the  result  of  diffusion  into  a  semi -infinite  solid  undi*r  constant  surface 
composition  comlitions .  Tliis  relationship  explained  how  was  dec rea sell 

with  partial  graphitization  without  any  change  in  This  relationship, 

together  with  the  *  2  relationship,  also  provided  quantitative  confirma¬ 
tion  of  the  equations  given  above  relating  the  iliffusivities  in  the  different  ilirec- 
tions  to  their  respective  jump  distances. 

Sime  apparent  anomalies  for  Tli  and  Ra  diffusion  in  PC,  reported  by  Wolfe, 
McKenzie  and  Borg  (17),  were  explained  by  the  model  develojxrd  in  this  investiga¬ 
tion,  as  were  some  of  the  difficulties  involved  in  studying  metal  diffusion  in  con¬ 
ventional  graphites.  The  diffusion  coefficients  of  .several  metals  in  PG  were 
show-n  to  lx;  sufficiently  large  to  permit  appreciable  diffusion-alloying  of  PG 
within  reasonable  times;  however,  the  efficacy  of  this  methoil  of  alloying 
remains  to  lx*  determined  by  experiment. 
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Appendix 

MATHEMATICS  OF  DIFFUSION  IN  LAYEIIED  STRUCTURES 


Using  various  arguments,  it  was  shown  in  Chapter  IX  that  metal  diffusion 
in  PG  occurs  only  over  the  surfaces  and  down  the  edges  of  the  PG  crystallites. 
Because  most  of  the  crystallites  in  PG  lie  parallel  to  the  deposition  surface, 
the  diffusion  coefficient  measured  parallel  to  the  deposition  surface, 
corresponds  to  that  parallel  to  the  layer  planes.  Conversely,  the  diffusion 
coefficient  measured  perpendicular  to  the  deposition  surface,  corres¬ 

ponds  in  some  fashion  to  a  boundary  diffusion  coefficient,  D^.  It  is  the  purpose 
of  this  appendix  to  examine  the  relationship  of  and  Dj^. 

Figure  27  shows  a  somewhat  idealized  model  of  metal  diffusion  in  PG  per¬ 
pendicular  to  the  deposition  surface.  From  the  arguments  prck'cntcd  in  Chapter 
IX,  the  volume  diffusion  coefficient  is  essentially  zero  so  the  only  paths  into  the 
PG  from  the  surface  are  the  major  cone  boundaries.  Approximately  perpendicu¬ 
lar  to  the  major  cone  boundaries  there  exist  high  diffusivity  paths  which  arc  the 
surfaces  of  the  crystallites,  and  they  have  a  spacing  equal  to  the  crystallite 
thickness,  .  Since  is  very  small  relative  to  the  concentration  profile, 
the  problem  can  be  treated  as  if  the  PG  were  an  anisotropic  solid  with  a  slab  of 
different  diffusivity  extending  into  it. 

Figure  28  shows  the  mathematically  idealized  structure  which  will  be  used 
to  treat  cone  boundary  diffusion  in  PG.  The  model  consists  of  a  single  boundary 
perpendicular  to  the  deposition  surface,  2f  units  wide,  and  e.xtending  in  the 
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MAJOR  CONE  BOUNDARY 


DEPOSITION  SLT^FACE 
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SLTtFACE 


-  £ 


L 
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+  £ 


Y 


OUTSIDE  BOUNDARY 
SLAB 


Fig.  28  Mathematical  Model  for  Diffusion  Perpendicular  to  Layer  Planes 
in  Pyrolytic  Graphite 
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x-direction  into  a  semi-infinite  solid.  The  distance  to  the  next  cone  boundary 
is  assumed  to  be  so  large  that  the  lateral  diffusion  can  be  treated  as  if  the  solid 
were  semi-infinite  in  the  y  and  z  directions  also. 

If  the  subscript  ^  refers  to  processes  outside  the  boundary  and  the  sub¬ 
script  1^  refers  to  processes  within  the  boundary,  then,  in  the  quadrant  con¬ 
taining  positive  X  and  positive  y  ,  we  have  two  differential  equations  to 
describe  the  diffusion  processes  outlined  above. 


ac 

a^c 

DE  1 

_ V 

=  D  — J 

at 

ay 

ac 

a^c 

DE  2 

i  — ^ 

=  ID.  - 7^ 

at 

•>  8x2 

.  dC 

V  ay 


y  =  £ 


(1) 


where  i  the  diffusion  coefficient  measured  parallel  to  the  deposi- 

tion  plane.  Note  that  a  term  involving  d  C^/dx  is  absent  because  the  volume 
diffusivity  in  that  direction  is  assumed  to  be  zero.  Assuming  that  there  are  no 
concentration  gradients  in  the  y-direction  within  the  boundary, 

Cy  =  (X.y.t) 

The  initial  conditions  are 

IC  1  (x.y.o)  =  0 

IC  2  Cj^  (x,o)  =  0 
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and  the  boundary  conditions  for  constant  surface  comF>osition  conditions  and 
semi-infinite  media  are 


BC  1 

Cjj  (o,t)  = 

C 

8 

ac 

1 

BC  2 

(x,«,t)  = 

V 

ay 

tl 

8 

II 

o 

BC  3 

Cj,  («,t)  = 

0 

The  coupling  condition  is 


CC 


(X.i.t)  =  Cjj  (x,t) 


Sectioning  experiments  involve  the  use  of  an  average  composition  taken 
over  a  slice  parallel  to  the  yz  plane,  so  we  can  define  an  average  composition 
by 


C  (x,t)  =  fCj^(x,t) 


Cy  (x,y,t)  dy 


(3a) 


and  this  can  be  differentiated  with  respect  to  t  yielding 

ao 


(3b) 
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t 


f 


From  DE  i  and  2,  and  BC  2 


a  f 

SrJ  Q,(x,y,t)dy  -  -  -r- 


at  /  'V 


V  ay 


a^c 


y  ■  I 


b 

»  3*2  ^  “v  ay 


y  m  I 


Subitltuting  into  Eq.  (3b)  wa  obtain 


8C  .  ID  ^ 
at 


Integration  with  respect  to  x  once  and  applying  BC  2  and  3  yields 


-ID, 


r  “d, 


b  ax  ^  at 


Integrating  again, 


•0  ao 

// 


^^^dn  dt 


Interchanging  the  order  of  integration  and  integrating  once,  we  obtain 


IDjjCjj(x.t) 


•  V 

// 
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(4) 


“  /  fir) 

X 


Applying  BC  1  to  Eq.  (4)  yields 


(5) 


At  this  point,  an  approximation  will  be  made.  It  was  observed  that  the 
concentration  profiles  for  Nb  diffusing  into  PG  perpendicular  to  the  deposition 
surface  could  be  represented  very  well  by  the  error  function  complement. 
Because  of  this,  we  shall  assume  a  solution  to  the  above  equation;  namely 

00 

C(x,t)  s  C*(x,t)  =  j  exp(-s^)ds  =  erfc  ^  (6) 

X  ^ 

^t 


♦ 

where  C*  is  the  observed  concentration,  is  the  surface  composition 
obtained  by  extrapolating  the  data  to  x  =  0,  and  D  is  the  diffusion  coef- 
ficient  measured  perpendicular  to  the  deposition  plane.  The  brackets  <> 
around  c  have  been  omitted  here  for  convenience.  Differentiating  Eq.  (6)  with 
respect  to  t  yields 


ac*  ^  S  X 
df  ■  7?"  rTirr  ®*p 

c 


2  \ 
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t 


and  from  Eq.  (5) 


'V.  ^  I 


«0 


8 _ I 

7r2t  J 


o 


(7) 


ID.  C  a  C*  D 
b  s  sc 

C* 

'I’b  “  r  “c 


(8) 


which  demonstrates  that,  approximately,  the  boundary  diffusion  coefficient,  D|^  , 

is  proportional  to  the  measured  diffusion  coefficient. 

Equation  (8)  is  only  approximate  because  C*  ,  the  measured  surface 

s 

composition  is  a  function  of  time.  Evaluating  Eq.  (3a)  at  x  =  0,  we  obtain 

oo 

C*(o,t)  =  fCg  +  j  C^(o,y,t)dy 

£ 

C*(o,t)  =  fCg  +  f(t) 

But,  the  approximation  given  by  Eq.  (6)  assumes 

C*(o,t)  =  C* 

which  is  not  a  function  of  time.  Thus,  the  approximation  given  by  Lq.  (6)  is 
not  valid  near  x  =  0. 
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C*(o,t)  can  be  evaluated  at  x  *  0  because 


Cjj(o,t)  *  «  C^(o,l,t) 


Cy<o.y.t)  * 


therefore 


•0 

C*(o.t)  =  1C,  +  2VT5;;i  I  C,  erfc  ^ 


2V  D  t 

C*(o,t)  =  ICg  +  Cg  jexp  (-IV4D^t)  -  I  ei 


i 


“=  2/I);t ) 


The  size  of  I  is  very  small,  being  of  the  order  of  a  few  A  ,  while  2VD^t 

g 

is  of  the  order  of  10  A.  Therefore 


CJ  =  C(o,t)  a  IC^ 


(9) 


Applying  Eq.  (9)  to  Eq.  (8) 
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But,  i  «  / D^t  ,  therefore 


(10) 


This  result  is  still  only  approximate  because  it  is  based  upon  the  assumption 

that  C*  is  not  a  function  of  t ,  and  Eq.  (9)  shows  that  C*  is  a  function  of  time. 

fl  s 

In  any  case,  all  of  the  above  derivations  demonstrate  that 

“c  “  '“b 

The  physical  significance  of  this  analysis  is  that  the  average  composition 
near  x  »  0  continues  to  increase  as  material  diffuses  out  the  first  lateral 
paths,  but  beyond  x  =  0  somewhere,  the  approximation  given  by  Eq.  (6) 
adequately  represents  the  data.  An  exact  solution  to  the  differential  equations 
describing  the  situatioii  was  not  found  and  the  above  approximation  was  made 
because  Eq.  (6)  represented  the  data  very  well  over  virtually  all  of  the  concen¬ 
tration  profile.  As  noted  in  the  body  of  the  text,  the  initial  concentration 

points  were  greater  than  the  extrapolated  value  of  C  and  this  was  attributed 

8 

to  a  lack  of  perfect  resolution  of  the  interface  with  the  lapping  technique.  It  is 
clear  fi'om  this  analysis  that  part  of  the  extra  concentration  at  the  very  beginning 
of  the  profile  is  due  to  the  results  shown  in  Eq.  (9). 
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The  differential  equationa  and  boundary  conditiona  given  above  can  alao  be 
approached  by  the  method  of  the  Laplace  tranaform.  The  following  diacueaion, 
due  to  Profeaaor  J.  C.  Shyne,  ueea  Laplace  tranaform  a,  but.  unfortunately  the 
inveree  of  the  laat  tranaform  ia  not  available.  Thia  approach  ia  preaented  here 
becauae  the  reault  ia  capable  of  being  numerically  evaluated;  however,  thia 
numerical  evaluation  haa  not  been  carried  out. 

The  equationa  and  conditiona  are  aa  before. 

Tranaform  ing  fir  at  with  reepect  to  t: 

t  ■  a 

Tranaform  ed 


Boundary  Condition 

Boundary  Condition 

11] 

c^(x,y,t) 

- 

C^(x,y,i) 

(2] 

Cjj(x,t) 

- 

Cjj(X,B) 

13] 

Cjj(«,t)  =  0 

- 

Cjj(«,e)  -  0 

14] 

Cjj(x,t)  *  C^(x,f,t) 

- 

Cj^(x,i)  =  C^(x,l,e) 

15] 

C^(x,«,t)  »  0 

- 

*  0 

(6] 

Cjj(x,o)  -  0 

- 

Cjj(x,o)  -  0 

17] 

C^(x,y,o)  *  0 

- 

C^(x,y,o)  *  0 

[8] 

Cjj(o,t)  = 

- 

C^^(o,e)  «  C^/S 
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Transforming  the  differential  equations,  DE  1  and  DE  2,  with  respect 


to  t 


•Cy(x,y,a)  -  C^(x,y,o)  « 


a^C^(x.y.s) 

dy^ 


eCjj(x,s)  -  C^(x,o)  «  Djj 


e‘Cjj(x,s)  aC^(x.y,s) 


dx 


2  i 


dy 


y  *  I 


Applying  conditions  [6]  and  [7]  above 


sC^(x,y,s) 


(11) 


•Cjj(x,s) 


D  dC 

_y _ V 

1  8y 


y  ■  £ 


(12) 


Now,  transforming  with  respect  to  x  : 


it  X  =  r 


Transformed 


Boundary  Conditions 

Boundary  Condition 

11] 

Cjj(x,s)  - 

^^(r. ») 

12] 

Cy(x,y,a)  ^ 

5y(r,y,a) 

[3] 

Cy(x,«,s)  =  0  - 

5y(r,«,s)  »  0 

[4] 

C^(x,s)  =  C^(x,£,s)  - 

5jj(r,s)  =  5y(r,l,s) 
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Transforming;  Kcjs.  (1  1)  and  (12)  with  !  C‘S|U‘cl  to  \ 


sC^(r,N,s)  - 


(»“  L\  (r,  v,s) 


(JV 


(13) 


sCj^(r,s)  =  D|^ 


2^ 


d  C|^((),  s) 


r  Cyr.s)  -  cyo.s)  -  — - 


1)^  c)C^(r,£,s) 
1  Oy 


(14) 


Now 


0^(0.  s)  =  — 


and 


9Cjj(o,8) 


9x 


=  </>(s) 


where  <p{8)  must  be  determined  later. 


Substituting  for  Cj^^o.s)  and  9Cj^(o,s)/9x  into  Eqs.  (13)  and  (14): 


rC 


D,  9  5  (r,£,s) 


d2  5 


iC  =  D 


V  V  .  2 
dy 
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Hcarranj?inn: 


c)C^,(r,r,s) 

<>y 


± 

D 


(15) 


(16) 


A  solution  for  Eq  (16)  is 


=  A  exp 


y 


+  B  exp 


Since  C^(r,oo,s)  =  0  ,  A  =  0  . 


and  since  C^(r,£,  s)  -  C.  =  B  exp  y  £ 

V 


B  =  C.  exp  yT  ^ 


and 


(17) 
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Now,  differentiating  Eq.  (17): 


Substituting  into  Eq.  (15) 


and  the  inverse  transform  of  with  respect  to  r  is: 

C 

C|jj(x,s)  =  cosh  yx  +  sinh  yx  (18) 
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The  function  (;p(8)  can  now  be  evaluated. 


BC.  C 

yx  +  <p(8)  cosh  yx 


3C^ 

3x 


=  0  at  X  =  0 


s  2 

0  =  —  y  tanh  y«  +  ^(s)  y 

B 


Thtrefore 


ip(§)  = 


y 


Substituting  this  relationship  into  Eq.  (18)  yields 

C 

Cjj(x,s)  =  —  exp(-yx)  (19) 


Transforming  Eq,  (17)  yields 


C^(x,y,8) 


(20) 
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The  problem  could  be  solved  if  the  Inverse  transforms  of  Eqs.  (19)  and 


(20)  were  available.  The  resulting  functions  C^(x,y,t)  and  Cj^(x,t)  would 
then  be  substituted  into  Eq.  (3a)  in  order  to  functionally  determine  the  observed 
compositicm-distance  data.  Unfortunately,  the  inverse  transforms  of  Eqs.  (19) 
and  (20)  are  not  available.  An  alternative  iq)proach  is  to  express  Eqs.  (19)  and 
(20)  in  the  complex  plane  and  evaluate  them  numerically;  however,  it  was 
decided  that  this  effort  would  not  be  worthwhile  at  present,  considering  that  the 
approximation  given  in  Eq.  (6)  fit  the  data  rather  well. 
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